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Beryllium-Uranium System 
R. W. Buzzard 


The phase diagram of the system beryllium-uranium was constructed from data obtained 


by thermal analysis, metallographic examination and X-ray diffraction 
characterized by one intermetallic compound, UBejs, melting at about 2,000° C, 


The system is 
a eutectic 


having a composition of approximately 0.12 weight percent of beryllium, occurring between 


uranium and UBe,;, at 1,090° C 
percent of beryllium 
bility of uranium in beryllium was detected 


are not affected by additions of beryllium. 


1. Introduction 


(he literature contains very little data on the 

bervllium-uranium system. This lack of information 
an be ascribed to the technical difficulties resulting 
from the high vapor pressure of beryllium and the 
large difference in the densities of the two compo- 
nents. A tentative phase diagram of the beryllium- 
iranium system was presented in the literature re- 
view by Buzzard and Cleaves. The diagram shows 
one compound, approximating UBe, in composition 
and two eutectics, one at 1,060° C between gamma 
iranium and UBeg, and the other at 1,200° to 1,250° 
( between beryllium and U Be, 

In order to have available a more accurate diagram 
of the bervllium-uranium system, work was initiated 
at the National Bureau of Standards at the request of 
the Atomic Energy Commission. A series of beryl- 
lium-uranium alloys was prepared with compositions 
ranging from 0- to 100-percent beryllium and analyzed 
metallographically, by X-ray diffraction, and by 
thermal analysis. 


2. Preparation of Alloys 


2.1. Stock materials, apparatus, and procedures 


The stock materials consisting of extruded beryl- 
lium and uranium were obtained from the Atomic 
Energy Commission. The alloys were prepared in 
a specially constructed, high-frequency induction 
furnace, in which the charges could > melted and the 
ingots poured in an atmosphere of purified helium. 
This was accomplished by enclosing the beryllia 
crucible and mold assemblies in a silica tube and by 
using bottom-pour crucibles. Details of the con- 
struction are shown schematically in figure 1. It is 
important that the beryllia be fired at approximately 
|,900° C and be of high purity. The pouring temper- 
atures were determined by sighting through the 
“sight”’ tube on the bottom of the stopper with an 
optical pyrometer. The stopper was essentially a 
closed end beryllia tube. 

Investigation sponsored by the Atomic Energy Commission at the National 
Bureau of Standards. 


hk. W. Buzzard and H. E. Cleaves; Issue No. 1, p. 26-27, Journal, Metallurgy 
Ceramics (TID 65), July 1948. 





, and a monotectic occurring at 1,550° C and 60-weight- 
The solid solubility of beryllium in uranium is slight; no solid solu- 
The transformation temperatures of uranium 
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2.2. Segregation 


The density of the beryllium-uranium compound 
was considerably greater than the beryllium-rich 
liquid phase, thereby creating a tendency for the 
compound to collect at the bottom of the crucible or 
mold as it precipitated from the cooling liquid. 
This settling-out tendency was not so “cag eager in 
the uranium-rich phase. It was possible to prepare 
beryllium-rich alloys reasonably free of segregation 
by pouring the ingots with the power on full into a 
cold mold of small cross section. The large supply 
of power, immediately before pouring, caused a 
vigorous stirring action, which kept the compound 
uniformly distributed. 

The distribution of uranium in the ingots in the ‘‘as 
cast condition” was checked both by chemical 
analysis (table 1) and microradiographic examination 
(figs. 2 and 3). Chemical analysis of the ingots 


Figure 1, 





showed satisfactory uniformity of composition from 
top to bottom, with the exception that the ingot in 
the 50-weight-percent range had a deviation of 5 
percent over its length. The variation from top to 
bottom in the other ingots was within | percent, and 
the ingots having compositions in the high-uranium 
and high-beryvllium fields showed little deviation in 
composition from top to bottom. The microradio- 
graphs showed good distribution of the compound in 
the bervllium-rich ingots, with the exception of the 
99-weight-percent beryllium ingot. The swirl dis- 
tribution of the compound shown in the photomicro- 
radiograph (fig. 2, A) is caused by the pouring action 
and is apparent in this ingot because of the low 
concentration of the compound. 





Ficure 2. 
beryllium-rich ingots in the 
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TABLE 
’ 
( aq] 
Alloy des hare 
mation inal ysis 
ignha of i 
‘ 
22B 0.2 
21B 0.5 
1B 
2B 1.2 
‘B 14.58 
1i6bB 190 
5B w.0 
4B 0.0 
7B m0 
aR “5.0 
17B 73 
9B 75. ¢ 
ISB 77.0 
15 85.0 
10B 95.0 
11B 1.0 
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“4B wo 05 
13B wo. 6 
* Wet chemical analysis for | 
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Fiacure 3. 


A, 82.4-weight-percent and B, 91.7 


3. Thermal-Analysis Apparatus and 
Procedure 


rhe solid-phase arrests were determined by time- 
mperature analysis, using a modification of the 
ipparatus described by Teitel and Cohen.® The 
ecessory apparatus included a Brown Air-o-line pro- 
vram controller to maintain uniform heating and 
ooling rates. A furnace atmosphere of purified 
argon was used in order to minimize the evapora- 
tion of beryllium. Arrests with this equipment 
could be determined with an accuracy of +2 deg C 
For regions involving a liquid phase, the high- 
frequency induction furnace used to prepare the cast 
alloys was modified for thermal analysis by substi- 
tuting a tapered crucible for the bottom- -pour-casting 
wible. A beryllia plug was inserted in the stop- 
wr tube to a height such that temperature observa- 
tions were made at the center of the melt, using a high 
precision optical pyrometer. The cooling rates were 
controlled by manually reducing the power of the 
convertor at a uniform rate in a manner that cooling 
rates of the order 0.5 to 2.0 deg/min were obtained 
through the thermal arrests. The solidus and solid- 
phase arrests of the beryllium-rich alloys were deter- 
mined in this furnace, using platinum-rhodium ther- 
mocouples enclosed in high-fired-alumina protection 
tubes and inserted through the sight tubes. The tip 
of the couple was seated on the top of the beryllia 
plug. The convertor control was modified so that 
the power could be reduced mechanically with an 
electric motor operating through a gear box. The 
\ime-temperature curves were recorded on a Brown 
electronic recorder. 


4. Results 


4.1. Beryllium-uranium compound 


Metallographic examination showed all of the cast 


oys (table 1) to consist of two phases. One phase 
creased in amount from either the uranium or 
ryllium end of the system until at approximately 
)-weight-percent beryllium, the alloy appeared to 
composed almost entirely of this phase (fig. 4, BE, 
Allovs in this range were hard and very brittle. 
rrans. Am. Inst. Mining Met. Engrs. 185 


Cohet 


el and Morris 


Photomicroradiographs of vertical sections of beryllium-rich ingots in the 


as cast condition” 


weight-percent beryllium X6 magnification 


It was found that this phase could be readily 
| isolated from a uranium matrix by electrolytically 
dissolving the uranium in a 10-percent chromic- 
acid solution. The separated phase was identified 
as the intermetallic compound UBes having a cubic 
structure, as reported by Baenziger and Rundle.‘ 
Precision determination of the lattice constant gave 
a value of a=2 5.1849. +0.0001 A at 27° C. 

In order to check the steichiomatele formula, 
UBe,s, a density determination and chemical analysis 
were made. The values obtained were 4.420+0.002 
g/em® at 26° C and 32.15 +0.05-weight-percent 
beryllium as compared to the computed values of 
4.373 g/em® and 33.0-weight-percent beryllium. The 
actual and computed values are in sufficient agree- 
ment to conclude that the formula is UBe,s and not 
UBe, as previously reported. The thermal-analvsis 
results indicate that the melting point of UBey 


| was approximately 2,000° C, 





| 


4.2. Beryllium-rich alloys 





Data (table 2) from the time-temperature cooling 
curves obtained for the beryllium-rich alloys show a 
monotectic at 60-weight-percent beryllium and 
1,550° C (fig. 5). The liquidus and liquid solubility 
temperatures above 1,600° C were influenced by the 
vapor phase of the beryllium, and the accuracy is 
not greater than +25 deg C. The maximum tem- 
perature indicated for the UBe,; composition is a 
selected value based on an evaluation of the thermal 
data, casting temperatures, and extrapolation of the 
liquidus curve. The solidus (1,280° C +3°) coin- 
cided with the melting point of beryllium at 1,283° 
C. A pronounced arrest occurred below the solidus 
at 1,250° C +10°, which can be ascribed either to 
impurities or an allotropic transformation in the 
beryllium phase, although the strong arrests observed 
| at this poimt indicate a possibility of the latter. 
Previous investigators * ° have observed similar 
double points in the melting-point determinations of 
beryllium and have concluded impurities to be the 
cause, but lacking conclusive proof as to the cause of 
the observed double arrest, the line is dotted on the 
diagram. 

‘N.C. Baenziger and R. E. Rundle, U1 
mission, AE C D-2506 April 8 1v4y 

G. Ocsterhe ld, Z. anorg. 0 allgem. Chem 


L. Losano, Alluminio 8 67 (1939 
H. A. Sloman, J. Inst. Metals 49 365 
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FIGURE 
Etchant 5 percent HF, A, B, C, D, E, G, and H at X100, F at x500. 
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9 8. . ; 
"ypical microstructures of beryllium-uranium alloys in the ‘as cast condition”’. 


A, 9.49-weight-percent beryllium, B, 98.86-weight-percent beryllium; C, 80.15-weight-per- 





cent beryllium; D, 60.7-weight-percent beryllium, dark constituent, U Bez, in light matrix of beryllium; E and F, 35.1-weight-percent beryllium, dark areas at low 
magnification are shown at high magnification to be beryllium (pale grey) and the mottled background is U Bew; G, 0.17-weight-percent beryllium, small angular 
areas, sharply outlined, are U Bex particles; H, 0.00-weight-percent beryllium. 


TABLE 2. 


wt % 
100 

w. 5 
w.0 
wm 5 
«0 
06.5 
96.0 


@.5 
92.5 
91.0 
9.0 
87.5 
85.0 
82.5 
80.0 
77.5 
75.0 
72.5 
70.0 
67.5 
65.0 
4.0 
61.0 
0.0 
57.5 
55.0 
%.0 
45.0 
40.0 
35.0 


* All points not otherwise indicated were determined in high-frequency induc- 
tion furnace, using an optical pyrometer. 
> Points determined in resistance furnace. 


Miscibility 
of La+L 


Liquidus 


Thermal arrests * 


¢ Points determined in high-frequency furnace with’thermocouple. 


Summary of thermal data obtained for beryllium-rich 
alloys 


arb—-Laté Solidus e+6-+n+6 
Lf c Cc 
> 1, 263 
> 1,24 
1,233 
1, 22 1, 223 
1,233 1,234 
f 1, 281 © 1,241 
{ 1,233 1, 228 
1, 2 
1, 271 
f °1, 280 © 1, 238 
\ 1,24 1, 239 
1, 236 © 1,240 
1, 20 
°1,239 
>1,m4 > 1, 267 
1. 47 
1, 550 © 1,282 °1,253 
1, 48 
1, 556 
1, 54 
1, 552 ° 1,270 1,237 


The uranium-U Bey eutectic is shown in an uranium matrix. 
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Figure 5. 
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Beryllium-uranium phase diagram. 
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letallographic examination showed that the cast 
‘tures of the alloys in the composition range 
to 100-weight-percent beryllium consisted of 
icles of UBe,,; embedded in a beryllium matrix. 
nogenizing treatments up to 1,225° C produced 
change in the structures. The compound 
eared in an angular form (figs. 4, A, B, D) over 
entire composition range, with the exception of 
arrow range between 92- to 80-weight-percent 
vllium, in which the compound was observed to 
e a globular form (fig. 4, C). All evidence indi- 
ed the solubility of uranium in beryllium was 
vy limited and was not detected in this investiga- 
I 
4.3. Uranium-rich alloys 


The cooling curves (table 3) obtained for the 

nium-rich alloys indicate a eutectic at approxi- 
itely 0.12-weight-percent beryllium and 1,090° C. 
Arrests were obtained at 770° and 660° C that 
respond to the two allotropic transformation 
temperatures, indicating little or no solid solution of 
beryllium in alpha or beta uranium (fig. 6). 


LE 3. Summary of thermal data obtained for the uranium- 
rich alloys 


lhermal arrests * 


Liquidu Solidus 


* All points not otherwise indicated were determined in high-frequency induc 
furnace, using an optical pyrometer 
Points determined in resistance furnace 


Metallographic examination showed that alloys 
n the vicinity of the uranium-UBe,; eutectic, as 
indicated by thermal analysis, contained a definite 
eutectic structure (figs. 4, G, H). Homogenization 
near the solidus spherodized the U Be. X-ray 
analysis of alloys quenched from 1,075° C in the 
range 0.027 to 0.079-weight-percent beryllium, placed 
the solid solubility of beryllium in gamma uranium 
between 0.042 to 0.079-weight-percent at 1,075° C. 
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Uranium-rich section of beryllium-uranium phase 
diagram. 


Figure 6 


5. Summary 


The phase diagram of beryllium and uranium was 
constructed from data obtained by thermal analysis, 
metallographic examination, and X-ray diffraction 
(figs. 5 and 6). The beryllium-uranium system is 
characterized by one intermetallic compound, UBe,,:, 
having a melting point of approximately 2,000° C. 
A eutectic, having a composition of approximately 
0.12-weight-percent beryllium, occurs between ura- 
nium and UBe,, at 1,090° C. A monotectic occurs 
at 1,550° C and 60-weight-percent beryllium. 

A slight solubility of beryllium in gamma-uranium 
is indicated, but no solubility of uranium in beryl- 
lium was detected. The transformation tempera- 
tures of uranium were unaffected by additions of 
beryllium. 


The author expresses his appreciation to John H. 
Darr, J. T. Sterling, and S. J. Rosenberg for their 
assistance in developing the thermal data; to R. A. 
Lindberg for the preparation of the beryllia refrac- 
tories; to H. C. Vacher, Ellen A. Buzzard, and Robert 
Liss for the X-ray and metallographic analysis; 
to Martha S. Richmond for the chemical analysis; 
and to Martha M. Darr for the spectrographic 
analysis. 


WasHINGTON, October 1, 1952. 
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Pulsed Light System for Multiple-Cell 
Ultracentrifuge Rotor 


Febri 


M. Wales, P. G. Sulzer, and L. C. Williams 


A system for photographing any spot on a moving rotor is described. 
triggered by a signal obtained from 2 smell alternator on the rotor shaft, is used. 
the phase of the light, any desired location on the rotor may be kept in view. 


A pulsed light, 
By shifting 
This svstem 


is used to record data from an ultracentrifuge having four symmetrically placed cells on the 


rotor. 


1. Introduction 


In certain types of ultracentrifuge measurements 

is desirable to investigate more than one sample 
at once. Beams ' has devised a method for making 

any sedimentation velocity measurements simul- 
taneously, in order to carry out the Gofman test 
for lipoprotein on a large number of samples. Sedi- 

entation equilibrium measurements often require 
. long waiting period for the attainment of equi- 
librium.? It would, therefore, be advantageous to 
carry out measurements on as many solutions as 
possible at the same time, especially when an extra- 
polation to infinite dilution must be made. The 
present paper gives an account of an arrangement 
that enables four symmetrically placed cells, all 
containing solutions, to be photographed separately, 
using the Svedberg * equilibrium machine. 


2. Principle of Operation 


The rotor of late models of this ultracentrifuge 
has four holes to accommodate three solution cells 
and a reference cell as shown in figure 1. The 
three solution cells are unequally spaced, radially, 
in order to separate the scale images photographed 
through each one. This arrangement has two dis- 
advantages: (1) the column of solution is so short 
that the regions of optical uncertainty near the bot- 
tom of the cell and the meniscus comprise a rather 
large fraction of the cell, and (2) great care must be 
used in filling the cells when maximum column length 
is desired, in order to obtain visible menisci. Fur- 
thermore, it has been suspected for some time that 
the use of a separate reference cell, where part of 
the scale is photographed through air to give refer- 
ence lines,‘ is unnecessary. Reference lines photo- 
graphed through the quartz windows in the upper 
unfilled) part of a solution cell gave reproducible 
reference readings. Therefore, four longer, sym- 
metrically placed solution cells were constructed,’ 
with the cell inserts ‘ of type 316 stainless steel and 

J. W. Beams, private communication 

M. Wales, J. Applied Phys. 22, 735 (1951 
'M W ales, F. T. Adler, and K. E. Van Holde, J 
‘ Tr. Svedberg and K 

dge, England, 1940 


Constructed by R. E. Ward and associates, NBS Shops Division, according 
1 design originated by E. Hanson, University of Wisconsin 


Phys. Colloid Chem. 55, 


©. Pedersen, The ultracentrifuge (Oxford Press, Cam- 
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4 sealer-counter circuit is included for speed measurement. 


the barrels of dural. Dural inserts had given cor- 
rosion troubles in aqueous solution and could not 
be permanently plated. Other types of stainless 
steel corroded in aqueous salt solutions. 

The four cells could be photographed separately by 
interrupting the light source in such a way as to 
illuminate the same cell at each flash, while all other 
cells would be dark. This could be accomplished, in 
principle, in several ways: (1) mechanically, by 
gearing a rotating disk with a window in it to the 
drive-shaft of the motor, and providing suitable 
means for phase change, (2) electrically, either by 
driving a rotating disk with a motor synchronized to 
the centrifuge shaft, or by pulsing a gaseous discharge 
tube in synchronism. The first method was not 
tried because of the difficulty of obtaining access to 
the centrifuge shaft without extensive changes in the 
design of the machine. However, electrical syn- 
chronization appeared promising, because the manu- 
facturer had installed a small coil and magnet alter- 
nator on the end of the shaft. This was intended to 
be used in connection with a pulse-shaping circuit, 
oscilloscope, and oscillator to measure the speed of 
rotation. It was found to be simpler and less ex- 
pensive to feed the output into a sealer circuit with a 
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iltracentrifuge rotor showing path of 
light. 
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Fiaure 2. 


scale of 64, driving a register, and to count the speed | 
directly. This arrangement was checked for in- 
ternal consistency by means of an applied 60-cycle 
pulse (a-c line), the time interval being measured by 
a synchronous automatic timer, also driven from the 
line. The circuit was also checked against standard 
frequencies available at the National Bureau of | 
Standards. 

From the 1:8 and 1:16 stages of the scaler, pro- 
vision was made to tap off the output through a 
cathode-follower circuit and amplifier, to drive a 
small synchronous motor. Unfortunately, the wave 
form obtained was so irregular that the motor would 
not run in synchronism except at a 1:1 frequency 
ratio, with the scaler out of the circuit. This made it 
inadvisable to use the motor because the speed range | 
of the centrifuge is 30 to 300 rps, and the motor could | 
be run only at 30 to 80 rps. 

The final arrangement, which proved to be success- | 
ful, is shown in figure 2. A Stroboluxpulsed light- | 
source assembly was driven with the amplified signal 

| 


from the centrifuge shaft, via the scaler circuit at a 
frequency ratio of 1:8 (or 1:16 if desired). Pro- 
vision for shifting the phase of the signal by + 180° 


was made in order to illuminate any desired cell. 





Circuit diagram, coupling of scaler output to Stroboluz. 


Cells were distinguished by means ofjidentifyin, 
dots painted on the far ends of the windows, awa, 
from the column of solution. The cell image was 
centered in the field of view in the camera by means 
of the potentiometers R, and R,, which were brought 
back to the vicinity of the camera. The images 
were remarkably stable and remained well centered 
for as long as 15 min without adjustment. 

Light from the Strobolux was concentrated on the 
entrance window of the centrifuge optical system by 
means of a truncated cone of kraft paper covered on 
the inside with smooth aluminum foil. This mate- 
rial is commonly used for insulation and is readily 
available. A thin frosted-glass diffuser was neces- 
sary for uniform illumination of the scale. A filter, 
designed originally to isolate the mercury green line, 
was also used in order to obtain nearly monochro- 
matic light. Under these conditions exposures of 
3 to 4 min were required to obtain good scale photo- 
graphs on CTC Kodaline plates. 

It is expected that this arrangement will materially 
increase the utility of the equilibrium ultracentrifuge, 


Wasaincton, December 1, 1952. 
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Damping of Elastically Supported Element 
in a Vacuum Tube 





No. 2, February 1953 Research Paper 2391] 





S. Levy, E. V. Hobbs, W. D. Kroll, and L. Mordfin 


it is shown that it is possible to damp the motion of an element in a vacuum tube by 


making the tube part of a two-degree-of-freedom system. 
The optimum values of the parameters are determined on SEAC, 


system are presented. 


The equations of motion of this 


the National Bureau of Standards electronic automatic computer, for a particular set of 


design requirements. 


In carrying out this computation, SEAC is used not only to compute 


the displacements of the system for given values of the parameters, but also to choose that 


set of parameters for which the performance is optimum. 


A description is given of a tube 


mount built in accordance with the computations and of the degree of damping achieved. 


l. Introduction 


In the development of the Ramberg vacuum-tube 
wecelerometer,' the problem arose of eliminating 
excessive electrical output in the neighborhood of the 

sonant frequency. One solution of this problem 
vas to include an electrical filter in the output of 
the accelerometer to electrically remove the effects 
of resonance of the mechanical system. As long as 
the resonant amplitude was in the linear range, 
such a filter was adequate. In some applications, 
however, large resonant amplitudes, beyond the 
range of linearity between mechanical displacement 
and electrical output, might be present, and for 
these it was considered desirable to damp the me- 
chanical vibration directly. 

The damping system described herein was espe- 
cially designed to give good damping for the Ramberg 
vacuum-tube accelerometer. The design require- 
ments were those for an accelerometer, namely, 
constant displacement amplitude and phase shift 
changing linearly with frequency, for frequencies 
from zero up to the natural frequency. The methods 
presented for achieving these design requirements are 
general and can be used where other types of damping 
are required 


2. Mechanical System 


A diagrammatic sketch of the mechanical system 
used to damp response of an elastically supported 
element in a vacuum tube is shown in figure | (a). 

A difficulty in damping the motion of elements of 
a vacuum tube by dampers attached to the tube 
base is that the elements are usually so light that 
their motion causes negligibly small motion of the 
tube base. To increase the motion of the tube base, 
and thus the absorption of energy from a particular 
vibrating element, the auxiliary mass m is used in 
the design shown in figure 1 (a). This mass is sup- 
ported on a spring k so chosen that the natural 
frequency of the mass, with the tube base clamped, is 
the same as the natural frequency of the element to 





*Walter Ramberg, Vacuum tube acceleration pickup, J. Research N BS 37, 391 
December 1946) RP1754. The measurement of acceleration with a vacuum tube, 
rans. Am. Inst. Elec. Engrs. 66 (1947). Samuel Levy, Recent developments 
xn 153 vacuum tube accelerometer, Proc. Soc. Exp. Stress Analysis, 
TX, No. 1 (1951). 
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Figure 1. a, Schematic diagram of mechanical arrangement 
for damping response of an elastically supported element in a 
vacuum tube; b, equivalent two-degree-of-freedom mechanical 
system. 


1, Auxiliary external mass elastically supported on cantilever springs; 2, tube 
envelope ; 3, elastically supported element; 4, case; 5, cantilever springs mounting 
tube in case. 


be damped. Under these circumstances, the re- 
sponse of the tube element and the mass, m, to 
motion of the tube base is the same. The mass m 
is large enough to force appreciable motion of the 
tube base, with a resultant damping of the motion of 
both m and the tube element by damper ec. 

The springs K support the tube from the outside 
case. 


3. Equations of Motion 


For purposes of analysis, the mechanical system of 
figure 1 (a) is idealized to that of figure 1 (b). The 
displacement (z—y) of mass m (and therefore of the 
tube element) to motions z of the case will new be 
determined. 

Writing the equations of equilibrium for masses m 
and M gives 


m(d*z/dt®) +k(z—y) =0. (1) 
M(d*y/dt®)+ K(y—z2)+k(y—2z) + e(dy/dt) —e(dr/dt) =0. 


These equations are solved on the assumption that 
the disturbance at the case is of constant amplitude a 
and radian frequency 


r=a sin wt. (3) 















The following substitutions are made 
I Mi er AK M (: ‘ V/ > uw M p 4 


After solving eq (1) and (2), it is found that the 
acceleration of the mass m is 


h y) m=—F(z—y aw*(A sin w+ B cos wt), (5 
where 

{ Or ( w'— wf wl — w* pf + 0 wy" (w 4 
wy \w 0 w* wf wD w pr A )° 

RB wyFQ*(w* ) + yw w* wt wth w pt - (re) 

; hdl # 
wy \w i i w wt w°th w* pt + O° )° 

The ratio of the maximum acceleration &(- Y) max 


of the mass m to the maximum acceleration aw’ of 
the case is 


Response ratio, R=, A’*+B (6) 


The phase shift of the amplitude (2—y) with respect 
to the motion « of the base is seen from eq (5) to be 


@=tan'(B/A),. (7) 


4. Application of SEAC 


It is desired to determine the values of p=m/M, 
Q=AK/M, and y=c/M for which the response ratio, 
eq (6), has a constant value near unity for values of 
the driving frequency w between 0 and 6. This was 
done by preparing a code for SEAC,’ which required 
it to carry out automatically the following successive 
steps: 

(1) For values of #/@=0.1, 2@=0.1, y/@=1.0, 
and p=0.6 compute the values of A and /?, eq (5), 
and from them the response ratio, 2, eq (6). Store 
R temporarily. 

(2) Increase the value of w/@ successively to 0.3, 
0.5, 0.7, and 0.9, repeating step (1) for each new 
value. 

(3) Determine the maximum, R,,,, of Roa, Ro.s, 

Ro.g and store temporarily. 

(4) Compute a criterion ( of flat response from 
C= Ry. +0.3Ro.9+0.1Ro.5+ 03Rinax— 1.43 and store 
temporarily. 

(5) Increase Q/@ in steps of 0.1 from 0.1 to 1.6, 
7/@ in steps of 0.1 from 1.0 to 2.0, and p in steps of 
0.2 from 0.6 to 2.0, and repeat steps (1) to (4) 
saving the results for the best five cases (those 
where the absolute value of ( is smallest). 

(6) Print the results for the five best cases. 

The best values were 


2/6=0.1, p=1.8, and y/@=1.9. (8) 


The curve of response ratio versus frequency ratio 
corresponding to these values is shown in figure 2. 


4 S3EAC, the National Bureau of Standards Eastern Automatic Computer, 
NBS Technical News Bulletin %, p. 9, (Sept. 1950 
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Figure 3. Phase shift for optimum values of design variables 
(y/@=1.9; p=1.8; 2=—0.1 


The response ratio is within | percent of 1.00 for 
values of w/@ less than 0.67. It is within 6 percent 
of 1.00 for w/@ less than 1.34. The phase shift for 
this case was computed from eq (7) and is plotted in 
figure 3. The departure from a straight line is 
negligible for w/@ less than 0.75. It should be noted 
that the additional points computed by hand for 
figure 2 fit in perfectly with those computed on 
SEAC. 

An important operation of SEAC for this problem 
was that it not only computed the responses but 
also, with the aid of the criterion C, chose the 
parameters for which the response curve was 
flattest. In this way, the plotting by hand and the 
visual comparison which might otherwise have been 
needed were avoided. 

In the coding of this problem for SEAC, no 
effort was made to include the linearity of the phase 
shift as part of the criterion C. This could have 
been done, with only minor changes in the coding, 
if it had been found necessary. 
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RE 4 Ramberg vacuum tube accelerometer as part of a 
»-deqree-of-freedom systen having optimum parameters 
lamping resonant amplitude and giving the accelerometer 


freq 


cency-response cur 





5. Damper Design 


\ damper designed as indicated in eq (8) was con- 
weted and is shown in figure 4 with the outer case 
(} removed. The masses m and .\V/, the springs & and 
and the damper ¢ are labeled on the figure. The 
ment of the Ramberg vacuum-tube accelerometer 
which damping was desired had a natural fre- 
wney @=2r(16l)rad/sec. The tube, together 
th the necessary base attachments, had a mass 
\/ 0.000285 Ib-see?/in. Using the optimum param- 
ters given in eq (8), the elements in figure 1 were 
lesigned to have the following values: 


} o\J—1.8(0.000285 0.000513 Ib-see?/in. 


kb m® = (0.000513)427(161)?=526 Ib/in ( 


+ \J—1.90\f = 1.912 7) (161) (0.000285) 
0.548 lb-see in 


K = MQ? = \1(0.180)* = (0.000285)427(16.1 
2.92 |b in. 


6. Electrical Output at Natural Frequency! 


To determine whether the damper was successfu 
in limiting electrical output at the natural frequency- 
the tube was subjected to a sudden change in accel 
eration. Such a change in acceleration should excite 
vibration at the natural frequency. Since the elec- 
trical output of the Ramberg vacuum-tube accelerom- 
eter is a direct measure of the displacement of the 
elastically supported element, the effect of the case 
on the vibration of the tube element could be con- 
veniently determined by recording the accelerometer 
output. The desired sudden change in acceleration 
vas obtained by mounting the tube at the end 

' a cantilever beam A, figure 5, whose tip was 
lefleeted sideways by a force, W, of several times 
the weight of the mounted tube, 7. Cutting the 

ipport wire, S, caused the desired sudden change 

| acceleration. 

The results of this test are shown in figure 6. It is 
ident from comparison of figures 6 (a), with the 
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Figure 5 Test set-up for determining the degree of damping 
of natural frequency achieved by making tube part of a two- 


de aree -of-freedom s system 
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Figure 6. Records of acceleration obtained in tests to determine 
effectiveness of damper in reducing response at 161 c/s natural 
frequency of tube. 









damper locked to prevent relative motion between 
the tube and the cantilever tip, and 6 (b) with the 
damper unlocked that a marked reduction in natural 
frequency electrical output results from use of the 


damper. In figure 6 (c) it is seen that the acceler- 
ation-time history obtained by fairing a curve 
through the natural frequency response in figures 
6 (a) and 6 (b) gives nearly identical results. This is 
expected since none of the output in figure 6 (a) 
exceeds the linear range of about 6 g for the 161 ¢/s 
model of the Ramberg vacuum-tube accelerometer. 


7. Accelerometer Performance 


Preliminary tests have been completed to de- 


(2) The National Bureau of Standards electron 
automatic computing machine, SEAC, is an excelle: 
means for computing the displacements of dynam 
systems, particularly where it is desired to comput 
the displacements for many different values of th 
mechanical parameters. 

(3) SEAC can be conveniently used as a substitut 
for hand plotting of results and visual inspection ; 
curves in selecting optimum parameters of a syste: 
by having SEAC compute, in addition to a solutio: 
of the problem, a criterion of quality and having i: 
select the solution for which the criterion is bes‘ 

(4) Further work is indicated to investigate the 


| cause of and to eliminate the large amplitude electrica| 


termine the extent to which the performance of the | 


accelerometer is in agreement with the computed 
curves, figures 2 and 3. In these tests, the accelerom- 
eter was driven by an electromechanical shaker at 
frequencies from 10 to 200 c/s, and the output of the 
accelerometer was viewed on an oscilloscope. It was 
observed that large amplitude electrical output at 
the natural frequency was present. This has been 
tentatively ascribed to nonlinearity. Other damping 
fluids and dash-pot designs will be considered in an 
effort to eliminate t is difficulty. 


8. Conclusions 


The following conclusions can be drawn from this 
investigation: 

(1) The vibration of an element in a vacuum tube 
due to external disturbances can be damped by 
making the tube a part of a two-degree-of-freedom 
mechanical system having suitable constants. 


output of the accelerometer when driven at its 


natural frequency. 
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Yglorimetric Determination of the Half-Life of Polonium 
D. C. Ginnings, Anne F. Ball, and D. T. Vier ' 


The heats of radioactivity of four samples of polonium have been measured with a 


Bunsen ice calorimeter over a period of about seven months. 


With samples ranging in 


initial powers from 0.17 to 1.4 watts, the half-life values calculated from these measurements 
were found to agree within 0.1 percent, or the equivalent of 0.0003 watt, whichever was the 


larger. 


of 138.3 days | 
ric method, 


l. Introduction 


Of all published determinations of the half-life of 
wlonium, probably the most accurate is that by 
Beamer and Easton [1],2 who used a calorimetric 
method and observed a sample for as long as 97 days. 
heir method consisted in measuring the temperature 
difference between a container with the sample of 
lonium and an identical container that was empty. 
Both containers were surrounded by an isothermal 
iacket, and the relation between temperature differ- 
ence and power was determined by means of electric 
calibration experimenis. The value of the half-life 
was found to be 138.3 days (+0.1%). A Bunsen 
ice calorimeter is also suitable for measuring radio- 
active power (and half-life) of radioactive materials. 
Because this method is entirely different from that 
ised by Beamer and Easton, and because much larger 
samples of polonium were available, it seemed 
worth while to measure the half-life of several sam- 
ples that were suitable for measurements with the 
ice calorimeter. 

The Bunsen ice calorimeter seems ideally suited for 
several reasons to the measurement of the heat 
evolved by curie or multicurie quantities of alpha, 
weak beta, and with modification [2], other radio- 
active materials. First, the heat leak is small (per- 
haps 0.0002 w). Second, the calorimeter requires 
very little attention during the measurements, 
usually a few minutes every hour or two for ice-bath 
replenishment. Third, no electrical or temperature 
measuring instruments are required because the 
measurement of heat requires only weighing of mer- 
cury. Fourth, the calibration factor of the ice 
calorimeter is a fundamental physical constant, 
which has been determined to about 0.01 percent by 
electrical calibration experiments [3]. 

On the other hand, the ice calorimeter in its present 
application has an inherent variation, or error, which 
is comparable with its heat leak. While this does 
not limit the precision of experiments involving about 
|-w power, it would limit the precision of measure- 
ments of much smaller powers. Other ice calorim- 
eters designed specifically for small radioactive 
powers [2] have proved sensitive to as little as 
0.00003 w. An improved ice calorimeter has been 


Los Alamos Scientific Laboratory, University of California 
Figures in brackets indicate the literature references at the end of this paper 
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The results with the sample with the largest power gave a half-life value of 138.39 
days, with an uncertainty of 0.1 percent (0.14 day). 


This is in agreement with the value 


+0.1 percent) reported by Beamer and Easton, who used a different calorimet- 


in use at the National Bureau of Standards for some 
time in the measurement of heat capacities at high 
temperatures [3]. Although this calorimeter was not 
intended for measurements of radioactive power, it 


| seemed suitable for measurements on certain samples 


of polonium (alpha emitter) furnished by the Los 
Alamos Scientific Laboratory, which initially de- 
veloped powers in the range of 0.17 to 1.4 w. By 
measuring the decay of radioactive power of a sample 
of polonium over time intervals comparable with its 
half-life, values were obtained of its half-life. 


2. Experimental Details 


The ice calorimeter used in these experiments is 
shown in figure 1. It is the same calorimeter pre- 
viously described [3]. An ice mantle (1) was frozen 
around the central calorimeter well (A) on a system 
of copper vanes (F) designed to increase the area of 
contact of the central well with the ice. The sample 
was suspended in the bottom of this well by a small 
wire that was pushed over against one side of the 
well by the gate G. The heat developed in the 
sample melted some of the ice, thereby decreasing 
the volume of the ice-water system and causing mer- 
cury to be forced into the calorimeter from beaker 
B. A small flow upward of dry helium in the calori- 
meter well was maintained to prevent condensation 
of water vapor from the room in the calorimeter and 
also to increase thermal contact between the sample 
and the calorimeter. 

The polonium samples were sealed in glass or 
metal containers, which were enclosed in brass outer 
containers made to fit closely to the central well of 
the ice calorimeter. Because the radiation from 
polonium is essentially all of the alpha type, the 
containers were completely effective in converting 
all the radioactive energy into heat within the con- 
tainer. A sample was lowered into the calorimeter 
(with the two platinum shields above it, as de- 
scribed elsewhere [4]), and about 10 or 15 min was 
allowed for it to come to essentially a steady tem- 
perature, the valve V being open. Then the mer- 
cury meniscus was adjusted to the upper part of 
the scale, C, by partly evacuating reservoir R. 
Valve V was immediately closed, thereby causing 
the meniscus at C to start to move downward due 
to the heat input to the calorimeter. The time wa. 
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Fieure | Schematic diagram of ice calorimeter. 


4, calorimeter well; B, beaker containing mercury; C, glass capillary; E, ice 
hath; F, copper vanes; G, gate; I, ice mantle; M, mercury; P, Pyrex containers; 
RK, mercury reservoir; T, mercury “‘tempering™ coil; V, needle valve; W, water 

1, 2, 3, and 4, thermocouple junctions 


recorded when the meniscus passed an arbitrary zero 
on this seale, the beaker of mercury (B) was re- 
placed by a weighed beaker of mercury, and valve 
V opened again. At the end of an experiment, the 
valve was closed again and the time recorded when 
the meniscus again passed its zero. The difference 
in the weights of the beaker gave the weight of mer- 
cury sucked into the calorimeter in a measured time. 
A small correction (perhaps | joule) was made for 
heat leak. The corrected mass of mercury was con- 
verted to energy by use of the calibration factor of 
the ice calorimeter as previously determined [3], 
270.46 abs j/g of mercury for the “ideal’’ calorimeter 
or 270.47 abs j/g for the actual calorimeter used. 
From the duration of the experiment (difference in 


times at which the mercury meniscus passed the | 


arbitrary zero), the energy was converted into 
average power. 

The values of power calculated in the manner de- 
scribed would be correct if the sample and all other 
parts in thermal contact with the calorimeter did 
not change temperature during the experiment. 
However, as heat is put into the calorimeter, the 
thickness of the water layer between the ice and the 
central well with its vanes increases, and as there is 
flow of heat across this water layer, the temperature 
drop across this water layer must increase. As a 


result, the temperatures of the central well and | 


vanes, and the sample with its container, all in- 
crease slightly during an experiment, so that part of 
the heat produced by the sample does not melt ice. 
This amount of heat depends on the power of the 
sample, the heat capacities of the various parts, and 
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on the thickness of the water layer at the beginnin 


of the experiment. The analysis of the correctio 
for this heat is given as follows. 

Consider a system (such as the radioactive samp 
plus metal calorimeter well with vanes) separat: 
from ice by a laver of water. Consider also that t] 
thermal contact between sample and calorimet 
does not change during an experiment. Let 

‘)= the total amount of heat (joules) that has be: 

introduced into the ice calorimeter at an 
particular time, 

P=the power (j/sec) developed by the sample, 

H=the heat capacity (j/deg C) of the system 

(sample plus calorimeter well with vanes 

A=the effective area of contact of the system 

with ice, em’, 

K=the thermal conductivity of water at 0° C 

0.0052 w em™'deg™, 
m=the average thickness of the water laver, em 
‘=the heat of fusion of ice, 333 j g. 
Then the grams of ice melted =) 0/333. 

The volume of water between the ice and the vanes 
is the same as the volume of ice melted because any 
void created by the difference in density of ice and 
water is filled almost instantaneously by water out- 
side the mantle passing through and around the 
cracks in the ice mantle. Hence, the volume of 
water layer formed = (specific volume of ice) (Q/333 

(1.09) (Q/333)=0.0033Q0; m=0.0033(Q/A); and 
the temperature difference across the water film 
Pm/(0.0052).A = (0.0033/0.0052)(PQ/A?) = 0.63(P0 
A’). Thus the heat stored in the system after ( 
joules have been introduced (assuming no tempera- 
ture gradient in the metal system) = (0.6377) (PQ/A®), 
and the heat stored in the water after Y joules have 
been introduced (assume linear temperature gradient 
across the film) = (4.18) (0.0033Q) (0.63/2) (PQ/A*) 
0.0043 (P(?/A*), and the total heat stored=(P/A’ 
(0.63H0+0.0043("). If we start an experiment 
with (, joules of heat already put into the calorim- 
eter and a steady state of heat flow from the sample 
to the ice mantle, and end the experiment with @, 
joules of heat in the calorimeter, and the same stead) 
state of heat flow, then the difference in heat stored 


(P/A*)(0.63H(Q.— Q,) +-0.0043 (G3 — Q7)I. 


The relative error, EZ, is 


— P 
 A(QG,—Q 


E= 7 (0.631 + 0.0043 (Q. + Q,)]. 


E , (0.631 (Q.—Q;) +0.0043(QG—QD], or 


If E is plotted against (+ 4,), the slope of the 
resulting straight line is 
> 


{* 


’ 


Slope = 0.0043 


and the intercept of this line with the F axis is 


HP 


Intercept = 0.63 Ar 
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is analysis shows that for any given power /’ of 
sample, the slope of this line is dependent only 
the effective area of contact (A) of the metal vanes 
h the ice. This effective area depends, of course, 
the length of the sample. The intercept of this 

e with the relative error axis depends not only 

on the area A but also upon the heat capacity, /7/ 

the sample-plus-vane system. 

For any given sample with a known heat capacity 

d length, values of A and //7 can be calculated, 

ing the dimensions and heat capacity of the cor- 

sponding portion of the calorimeter well and its 
stem of vanes. It is believed that a better evalu- 
ion of A and 7 can be obtained by electrical cali- 
ation experiments, putting in known quantities of 
ectric heat into the calorimeter over a length 
juivalent to the radioactive samples to be measured 
lf the heat is distributed in the electrical experiments 
he same as in the experiments with the radioactive 
umples, the value of A obtained from the electrical 
yperiments should also be the same. However, the 
ilue of 77 obtained from the electrical experiments 
liffers from that for the radioactive experiments by 
the difference in heat capacities of the electric heater 
ind the radioactive samples. 

The results of the electric calibration experiments 
ure indicated by the circles in figure 2, which shows 
the relative error, FE, plotted against ((,+,). 
The best straight line through these circles is labeled 
‘Electric Experiments” and was determined by the 
method of least squares to be /=0.000633 + 2.992 > 
1O-*(,+Q2), where 9, and © are expressed in 


joules. The electric calibration experiments were 


made with a power of about 1.4 w to correspond 
approximately to the initial radioactive power of 
sample 4, which was the largest of the four samples 
and therefore had the largest correction for the water 


laver The slope of the solid line (electric experi- 
ments) indicated that the effective area of contact 
A) between the ice and the source of heat was about 
120 em*. This is roughly equivalent to a 9-cm 
length of the calorimeter well, as compared to an 
actual length of the electric heater of about 7 em. 
While the value of A obtained from the electric 
calibration experiments agreed as well as expected 
with that estimated from the dimensions, the value 
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Figure 2. Corrections for water layer in calorimeter. 


Electrical experiments: &, 10-11-49; @, 10-12-49; ©, 10-14-49 


of heat capacity 77 (after correction for the heat 
capacity of the heater corresponded to about a 
17-em length of the calorimeter well, which is con- 
siderably more than expected. Some of this differ- 
ence may be due to the assumption in the theory 
that the temperature gradient in the copper vanes ts 
negligible. If the position of the line had been cal- 
culated entirely on the basis of a 9-cm length of the 
calorimeter well, it would have the same slope as the 
line given by the electric experiments but would be 
shifted to decrease by almost 0.0003. Because it is 
believed that the electric calibration experiments give 
more reliable values of the effective A and //, the 
solid line in figure 2 was used as basis for the correc- 
tions for the experiments with the radioactive 
samples. The dashed lines shown in the figure for 
sample 4 were calculated on this basis, as well as 
corresponding lines (not shown) for the other samples 
Values of apparent calibration factor (which is 
equivalent to 270.47 (1+/)) were read from the 
dashed lines for the corresponding values of (Q,4+-Q.). 


3. Results 


The results of the measurements with the four 
samples are given in table 1. The mean time of 
experiment is given, so that all experiments in one 
group can be corrected for the radioactive decay to 
bring them all to the same time (noon, EST) on the 
reference date. The mass of mercury was corrected 
for the effeet of heat leak between the calorimeter 
and its surroundings (this effect is usually only a 
few milligrams of mercury). The quantities Q, and 
(J. are as previously described. Using the quantity 
@-+- Qs, the apparent calibration factor of the ice 
calorimeter was obtained from figure 2, from the 
curve for the appropriate sample. The product of 
this factor and the mass of mercury, divided by the 
duration of the experiment, gives the power at mean 
time of experiment. Using an approximate value of 
the half-life of polonium, this power was converted 
to power at noon (EST) on the reference date. The 
deviations from mean power are the deviations of the 
results of the individual experiments from the mean 
power of the group of experiments. The values of 
half-life given in the last column are based on the 
mean powers and reference dates, and the values 
indicated in parentheses are the authors’ estimated 
uncertainties of the half-life based on both accidental 
and systematic errors. 

It will be noted that more measurements were 
made with sample 4 than with the other three sam- 
ples. This was because sample 4 had the highest 
power of the four samples, and it was believed that 
the highest relative calorimetric accuracy could be 
obtained with this sample, resulting in the best de- 
termination of the half-life of polonium. The power 
of sample 1 was so low that calorimetric errors could 
account for the low value of half-life calculated. 
The uncertainty in heat leak of the calorimeter was 
believed to be about 0.0002 w. This uncertainty, 
when compared with the power measured for sample 
1 on May 14, 1950 (0.05698 w), could explain most 
of the discrepancy with the results with the other 









































Tasie l. Results of experiments 
| 
| Power at Power at Devia- 
Apparent + 
Mean time of | Duration ’ ~~ mean noon (EST) iis tions : 
Date of experi experiment of experi M —_ of an U calorimeter | time of on Reference from Mean Half-life 
ment Par mercury factor date power 
EST ment M0471+E) experi reference mean 
— : ment date power 
Sample 1 
sec g gHe og He ig abs*u abs w percent abs u Days 
10-349 f 32pm $ 557.8 2. 3048 0 2 270. 49 0. 17523 «0. 17451 +0. 16 
aM | 6236p. m 12 367.2 7. ORE 2 10 270. 49 17474 17411 1G~ 4-49 | 0.07 1 6 a7493 
10-4-49 f 1232 p.m 9 928.0 6. 39% 10 17 270. 0 17417 17419 ‘ : | —0. 02 — 
: \ 3:20 p. m 1 126.5 6. 5205 17 23 270. 3 17418 17430 +4). 04 
137.7 
53%) L:14p. m D 34.0 4. 5260 a4 Ps] 270. 48 OO016 O54 | +0). 07 &’ (41.0 
4) 1:19 p. m 19 440.0 4.2913 2 34 270. 48 05071 05680, { 0). 32 | on 
eh f 216p.m 21 909. 5 4. 6221 0 ‘ 270. 48 O5706 05709 14-50 +0. 19 0. 05608 
it) | 74pm 13 948.0 2. 9396 5 8 270. 48 05700 05709 | +0.19 
Sample 2 
| il: Wa. m 6 168.5 @ 1630 23 32 270. 4% 0. 40190 0. 40183 | | +0. 04 | 
M 4y 12:45 p. m 5 224.0 7. 7486 32 0 270. 57 4133 40139 W- 549 —0. 06 0. 40165 
| 2:15 p. m 5 478.0 Ss. 1204 0) 4s 270. 58 wis W173 | | +0. 02 | | 138. 52 
> (+0.4 
+-27-%) 1:48 p. m 16 514.9 8.8470 4 33 270. ® 14491 14425 | on_* f+o.18 | , 
4+-28-™) 1.00 p. m 17 900.8 9. 5623 43 43 270. 50 14370 14373 j #2350 | lo ig =f 0. 14300) 
Sample 3 
{ L24p.m 5 134.4 12. 25 0 12 270. 58 0. 643° 6 0. 64323 | | +0.02 | 
10 &-49 2:45 p. m 5 207.5 12. 3623 12 24 270. 61 4241 (4278 10- &49 0. 05 » 0. 64312 
| 406 p. m 4 10, 8695 24 35 270. 63 64280 64335 | l+o.o4 | | ainied 
138. 52 
1-19-50 1:49 p. m 12 692.3 417 0 270. 51 24471 2358 | 0.01 | "GAS 
+-20-5 f 11:35 p.m 10 313.6 9. 2808 il 20 270. 52 24367 24365 20-5) +0). 02 0. 24361 
oe \ 231 p.m 10 808. 6 9.7279 2 30 270. 52 24347 236) «=| 0.00 | 
Sample 4 
12:55 p. m 2 634.7 14. 0020 a 42 270. 84 1. 43987 1. 43307 —0. 01 
| 1:35 p. m 2 76.9 14. 7006 42 57 270. 91 1. 43004 1. 43382 +0. 04 
i449 1:23 p.m 2 622.8 13. 8785 0 “4 270. 72 1. 43251 1. 43206 —0). 02 
| 2:07 p. m 2 676.2 14.1607 | 148 270. 78 1. 43275 1. 43345 +0. 01 
2:57 p. m 2 733.6 14. 4582 P-2) 42 270. 84 1 1. 43538 10- 7-49 +0. 01 1. 4332 
3:43 p. m 2 717.2 14. 3735 42 57 270. 91 1 1. 43423 +0. 07 en 
10-749 | 4:28 p.m 2 736.5 14. 4643 57 71 270. 97 1 1. 43361 +0. 02 
P 4:14 p. m 2 680.3 14. 1840 71 sth 271.03 1 1. 43105 —0. 15 138. 39 
6:00 p. m 2 879.7 5. 2113 Ne) 101 271. 0 ? 1. 43382 \ +0. 04 > (+0.14 
11:00 a. m 3 170.1 10. 2320 10 2 270. 64 0. 87335 +0. 05 
11:4 a. m 3 345.9 10. 7884 2 31 270. 67 87272 —0. 03 
12:50 p. m $ 350.8 10. 8281 31 42 270. 70 7 s 0.04 x 
1-14-50 144 p. m 3 067.7 9.8837 | 42 52 270. 73 1-14-50 0.04 | 0: 87204 
2:35 p. m $115.1 10. 0448 52 62 270. 74 348 +0). 06 
3:23 p.m $178.3 10. 2384 fi2 72 270. 78 87202 +0. 00 138. 40 
6 (+0. 18 
10:4 a. m. 8 282.3 14. 8070 3 18 270. 4 48370 48350 +0). 00 
12:25 p. m 8 305.9 15. 0148 18 33 270. 58 . ASSAY ASSUS +0. 10 
2:46 p. m 8 54.5 15. 2198 33 45 270. 61 48257 48284 R108 —0. 12 —_ 
12-50 5:00 p.m 7 430.2 13.2613 48s 270. 63 48243 48203 12-530 |) 6.10 || 0: 48342 
6:55 p. m 6 0.2 11.1775 61 72 270. 64 48322 48391 +0. 09 
8:43 p. m 6 047.3 12. 3852 72 5 270. 66 48252 48339 0.00 
* Experiment weighted one-third because of short duration of experiment 
» Figure given for uncertainty based upon authors’ estimate of both accidental and systematic errors. 
samples having higher power because it would re- | probably comparable with all other errors. In addi- 


quire an increase of only 0.6003 w in the value ob- 
served on May 14 to agree with the other samples. 

It is interesting to note the corresponding effects 
of an uncertainty of 0.0003 w on the other samples | 
having higher powers. With samples 2 and 3, this 
uncertainty in the May values would result in half- | 
life uncertainties of about 0.21 percent (0.29 day) 
and 0.12 percent (0.17 day), respectively. With 
sample 4, nell the same power uncertainty in 
the values on January 14 results in only 0.04 percent 
(0.06 day) uncertainty in the over-all half-life. | 
The uncertainty in heat leak is probably the largest | 
error in the measurements with samples 1, 2, and 3. | 
With sample 4, however, the heat leak uncertainty is | 
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tion, sample 4 was observed at an intermediate date 
(0.7 half-life) in an effort to detect any change in the 
decay constant with time. 

The measurements with the four samples agree 
within +0.1 percent, or 0.0003 w, whichever is the 
larger. The slightly different half-life values on 
sample 4 for the two different periods are not 
significant because the difference is much less than the 
experimental error. There is no significant evidence, 
therefore, of a change in the decay constant, such as 
would occur either if the polonium was contaminated 
with radioactive impurities having values of half-life 
different from polonium, or if some secondary chemi- 
cal or nuclear reaction produced heat that was not 





ectly proportional to the radioactive energy of the 
onium sample. In addition, the polonium sam- 
s used were of a purity that should have precluded 
: possiblity of appreciable contamination by other 
lioactive elements. The polonium had _ been 
rified both by distillation and by electrodeposi- 
yn on platinum foils from dilute nitrie-acid solution 
polonium nitrate. Of course, the stable lead 
stope formed from the decay of polonium was pres- 
t but could not affect the results of this investiga- 
mn. 
The best value of the half-life of polonium was 
timated from the two values for sample 4, giving 
ore weight to the earlier value, where the sample 
id a larger power. On this basis, the authors be- 
ve that the best value of the half-life derived from 
these measurements is 138.39 days. It is believed 
that the calorimetric uncertainties (based upon both 


accidental and systematic errors) in the experiments 
may result in an error in this value of 
(0.14 day). 

The results of these measurements are in agree- 
ment with the value reported by Beamer and Easton 
[1] of 138.3 days +0.1 percent, determined calori- 
metrically by another method. 


t+0.1 percent 
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Means of Alkaline 


Determination of Molecular Weights of Dextrans by 


Copper Reagents’ 


H. S. Isbell, C. F. Snyder, N. B. Holt, and M. R. Dryden 


Reducing values were determined for samples of clinical dextrans, 
isomaltose with several alkaline copper reagents 


gentiobiose, and 
The results for the dextrans were con- 


verted to number-average molecular weights on the assumption that equimolecular quanti- 
ties of dextran, isomaltose, and gentiobiose have approximately the same reducing power 


The molecular weights, thus obtained, are in 


recently developed radioactive cyanide method 


fair agreement with values obtained by the 
The results show that copper-reducing 


methods are suitable for the determination of the number-average molecular weights of 


clinical dextrans 


1. Principles of the Method 


There is widespread interest in degraded dextrans 
prepared from material formed by the action ol 
Leuconostoc mesenteroides and other bacterial agents 
on sucrose. Various fractions of these carbohy- 
drates are under investigation as blood-plasma ex- 
tenders. Therefore, simple and reliable methods of 
analysis of these clinical products are desirable. 
Dextran is composed of glucose units in 1,6-glycosidic 
linkage, with occasional branching at carbons 3 or 4 

1, 2, 3, 4].2 The terminal reducing unit consists of 
a glucose molecule substituted at carbon 6 with an 
alpha glucosyl group. Ordinarily, each molecule of 
clinical dextran has one, and only one, reducing end 
group. This condition arises from the manufactur- 
ing process, which ordinarily involves hydrolytic 
cleavage of a high-molecular-weight polysaccharide 
5). Each cleavage results in two molecules. One 
contains the reducing end group originally present; 
the other contains the reducing end group formed by 
cleavage of a glycosidic link. Thus there is one 
reducing group per molecule, and the number- 
average molecular weight * can be calculated from 
the number of reducing groups per unit weight. 

The reducing end groups are largely responsible 
for the capacity of the polysaccharide to reduce 
alkaline copper and other reagents. For determina- 
tion of reducing groups in polysaccharides of rela- 
tively low molecular weight, alkaline copper reagents, 
iron salts, and dinitrosalicylic acid have been most 
widely used. These reagents do not react stoichio- 
metrically, and factors must be derived for convert- 
ing the reducing values to molecular weight. For 
the determination of the molecular weights of 
starches, Meyer, Noelting, and Bernfeld [6] applied 
Sumner’s dinitrosalicylic acid method and standard- 
ized the results with maltose. The reducing end 
eroup of dextran differs from that of starch in that 
it has a 1,6-a glucosyl link in place of a 1,4-a-glucosyl 
link. Isbell, Pigman, and Frush [7] have shown 
that the position and the configuration of the gluco- 


rhis work was supported, in part, by the Medical Research and Develop- 
ent Board, Office of The Surgeon General, Department of the Army. It 
was carried out under the Plasma Volume Expanders Program of the Sub- 
mmittee on Shock, National Research Council. The project was sponsored 
y the Office of the Surgeon General, U. 8. Department of the Army 


' Figures in brackets indicate the literature references at the end of this paper. 
The number-averege molecular weight, Myg=2Z(niM,)/Zng, where ng is the 
mber of molecules of molecular weight Mg 
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sidic link in a disaccharide influences the molecular 
reducing power, and that a parallelism exists in the 
molecular reducing power of substances of the same 
or similar structure. The reducing power per mole 
cule of a disaccharide having a 1,6-link is substan- 
tially higher than that of the corresponding 
monosaccharide. Presumably the higher reduction 
arises from cleavage of the disaccharide link in the 
partially oxidized product during the analysis with 
alkaline copper reagents. Glycosidic links in poly- 
saccharides near the reducing groups would be broken 
more readily than links one or more hexose units 
removed from the reducing end group. For this 
reason the hexose units more than one unit removed 
from the reducing end group would not contribute 
appreciably to the reducing power, and a molecule 
of dextran should have the same reducing power as 
a molecule of isomaltose. Thus the molecular weight 
of a clinical dextran might be calculated from the 
reducing power of isomaltose (6-a-p-glucopyranosyl- 
p-glucose). Isoma!tose has not been crystallized, 
and the amorphous material available is of doubtful 
purity. Therefore, for most of the comparisons 
given in this paper and those given in [8] and [9], 
the closely related substance, gentiobiose, (6-8-p- 
glucopyranosyl-p-glucose) was selected as a reference 
standard, even though it has the beta linkage rather 
than the alpha linkage common to dextran and 
isomaltose. 

In the present study copper-reducing values were 
measured for samples of clinical dextrans and genti- 
obiose by several methods, and the number-average 
molecular weights of the dextrans were calculated on 
the postulate that a molecule of dextran has the 
same reducing power as a molecule of gentiobiose 
Table 1 gives the molecular weights thus obtained, 
together with those estimated by the recently de- 
veloped radioactive cyanide method [10]. The 
results obtained by different copper-reducing methods 
are in fair agreement with one another and with 
those obtained by the evanide method. The results 
show the practicability of determining number- 


average molecular weights by copper-reducing 
methods, but further refinements in the method 


seemed desirable. 
Following the suggestion of one of us [8] that 
equimolecular quantities of dextran and isomaltose 












TABLe 1. 


Somogyi 


Sample PO, 


, 100 
5, 300 
. 300 
5, 500 
34, 700 


A verage 


* Percentage deviation from arithmetic mean 


have like reducing powers, workers in other labora- 
tories developed analytical procedures based on an 
isomaltose standard. A group at the Commercial 
Solvents Corp. applied Somogyi’s method [11] to 
dextran and isomaltose. With a 20-min heating 
period they found that 1 ml of 0.005 N copper reagent 
is equivalent to 0.222 mg (or 0.000649 millimoles) of 
isomaltose [12]. Recently, the Northern Regional 
Research Laboratory, U. S. Department of Agri- 
culture, kindly furnished us a sample of substantially 
pure but amorphous isomaltose. A study with the 
Somogyi phosphate method revealed that the re- 
ducing power of this sample of isomaltose was only 
94 percent that of gentiobiose. A lower reducing 
power for the disaccharide having the alpha con- 
figuration is in accord with the relative reducing 
power of maltose and cellebiose. Because dextran 
and isomaltose have the alpha glucoside structure, 
whereas gentiobiose has the beta structure, iso- 
maltose would be preferred as a reference compound. 
Nevertheless, one hesitates to standardize a method 
on an amorphous sugar. For this reason, it seems 
desirable to select an arbitrary factor until pure 
crystalline isomaltose is generally available for use 
in standardization. The average of our values for 
the reduction equivalents of isomaltose and gentio- 
biose is 0.00065 millimole of sugar for each mulliliter 
of 0.005 N reagent. This is » to that obtained 
by the group at Commercial Solvents Corp. for 
isomaltose. Use of this figure gives the following 
equation for converting copper-reducing values ob- 
tained by the Somogyi method to molecular weight: 
M,=—w/(V 0.00065), where MM, is the number- 
average molecular weight of the dextran, w is the 
weight of sample in milligrams, and V is the net 
titration in milliliters of 0.005 N solution. 


2. Experimental Procedures 


2.1. Preparation of Samples 
a. Gentiobiose 


Gentiobiose methyl alcoholate was prepared by 
the method described on page 463 of reference 13. 
The alcoholate was converted to the pure anhydrous 
beta modification of sugar by dissolution in water, 
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Com parison of number-average molecular weights of dextrans obtained by various methods, using gentiobiose as a standard 


Folin-Wu ds Mean 


somenst ds 


Oy 


44, 000 
43, 200 
40, 000 
35, 600 
34, 200 


40, 800 
40, 900 
40, 900 
32, 800 
31, 600 


36, 200 
35, 300 


73, 
33, 
19, 
31, 
35, 


700 
400 
500 
200 
300 


71, 700 
, 400 
, 700 
29, 700 
5, 700 


76, 400 
20, 900 


32, 500 


34, 100 


followed by evaporation of the water and alcoho! 
under reduced pressure. The resulting sirup crys- 
tallized spontaneously. The material was recrys- 
tallized by dissolving in water, concentrating the 
solution to about 90-percent total solids, adding 5 
parts of methyl cellosolve, and seeding with 8-p- 
gentiobiose crystals. The crystals that formed were 
separated, washed with methyl cellosolve, and dried 
under reduced pressure at 60° C. At equilibrium 
la]? was found to be +9.60 (water, C=2)* in 
agreement with the reported value of + 9.6. 


b. Isomaltose 


The isomaltose was supplied by Allene Jeanes of 
the Northern Regional Research Laboratory. It 
was a carefully purified amorphous sample originally 
prepared from the octaacetate. It contained by 
analysis 41.1 percent of carbon and 6.6 percent of 
hydrogen. The carbon content corresponds to a 
product containing 97.6 percent of the anhydrous 
disaccharide. Hence the dry weight of the material 
used for the determination of reducing power was 
based on the assumption that the material contained 
2.4 percent of moist «re. 


c. Dextran samples 


The samples of clinical dextran received in pow- 
dered form were dissolved in water, filtered, and 
lyophilized (freeze dried). Dextran received in 
“saline’’ solution was purified by two precipitations 
with methanol. The resulting material was dis- 
solved in water, filtered, and lyophilized. For the 
analytical measurements lyophilized samples were 
placed in tared volumetric flasks and dried to con- 
stant weight at room temperature under less than 
0.1-mm Hg pressure in a desiccator connected with 
a dry-ice—acetone trap. After the weight of the dry 
sample had been determined, water was added and 
the solution was made to volume. For the analyses, 
suitable quantities of the solution were transferred 
with a pipette to the reaction tubes or flasks. 


2.2. Modified Scales Method 


_ The reagents and procedure were the same as those 
given on page 189 of reference 13, except that a 


* 2 grams of gentiobiose dissolved in water and made up to 100 ml. 





Reducing values of gentiobiose by the modified Scales 
method [7] 


10-min boiling time 


Net titra- 
tion 0.04 
N iodine 


Sample Sugar/ml of titration 


ml Millimoles 
5. oS 0. 00485 
5.4 00490 
11.78 00493 
14. 00497 
17.39 00503 


10-min boiling time was used in place of 6 min. 
The reduction values found for gentiobiose are given 
in table 2. 

With samples of gentiobiose the reducing power 
changed with concentration. Hence a curve was 
made by plotting the millimoles gentiobiose per 
milliliter of 0.04 N reagent with respect to the net 
titration in milliliters. Readings from this curve 
were used as factors to convert the net titrations 
obtained with dextran samples to gentiobiose equiv- 
alents. The following formula was used to convert 
the results obtained with clinical dextran to molecular 
weights: M,=w/(VX/), where w is the weight of 
sample in milligrams, f is the millimoles of gentio- 
biose equivalent to each milliliter of the titration 
read from the curve, and V is the net titration in 
milliliters of 0.04 N solution. Some results obtained 


with clinical dextran samples are given in table 3. 
TaBLe 3. Number-average molecular weights of clinical 
dextrans obtained with modified Scales method 


(Boiling time 10-min; total heating period 14 min 


J, milli- 
moles of 
| tentiobiose 
| per milli- 
liter 


V, net 
titration 
0.04 N 
solution 


Dextren sample w, weight 


number of sample Ma 


ml 
10. 37 
10.4 
6.10 


0. 00491 19, 200 
00492 19, C0) 
00486 33, 800 
00484 33, 800 
00486 30, 500 
00486 30, 400 
00483 45, 400 
00483 46, 100 
00483 5, 900 
(484 36, 600 
00486 33, 400 
OO485 34, 900 
00486 33, 400 
00485 33, 500 
00485 34, 000 
00486 3°, 100 
00484 |, 700 
00484 . 600 


SCOe COUN ee we 


2.3. Somogyi Phosphate Method 


The procedure and reagents used were essentially 
the same as those described by Somogyi [11] but are 
repeated in part, since the original article gives 
certain alternative procedures. 

The alkaline copper reagent was prepared from 
70.6 g of Na,HPO,.12H,0, 40 g of Rochelle salt, 100 
ml of N NaOH solution, and 80 ml of solution con- 
taining 8 g of CuSO,.5H,O, 180 g of Na,SO,, and 
water to make | liter. For some of the measure- 
ments a measured volume of the sugar solution was 
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heated with an equal volume of the copper-phosphate 
solution. After cooling, a precisely measured 
amount of a standard potassium iodate solution and 
an appropriate amount of potassium iodide and sul- 
furic acid was added, and after complete oxidation 
of the cuprous copper the residual free iodine was 
titrated with 0.005 N sodium thiosulfate. In other 
measurements the “combined” copper reagent was 
used. This reagent is prepared by incorporation of 
the requisite quantity of potassium iodate in the 
alkaline copper solution. Two concentrations were 
used; one contained 2 ml and the other 1 ml of N 
potassium iodate per 100 ml of the combined copper 
reagent. 

The analyses were made in 25X250-mm test tubes 
in groups of eight. Five-milliliter quantities of the 
solutions containing from 0.002 to 0.01 millimoles of 
gentiobiose, isomaltose, or dextran, or 5 ml of water 
for blank measurements, were transferred to the test 
tubes. Each sample was mixed with 5 ml of the 
copper reagent. When the “combined” reagent was 
employed, its volume was measured precisely. The 
tubes were closed with funnel-shaped stoppers and 
placed in a rack suspended in a boiling-water bath. 
After 25 min the rack and tubes were transferred to 
a cold-water bath and brought to room temperature. 
The tubes were then treated separately and succes- 
sively. In measurements using the unmodified cop- 
per reagent a precisely measured volume of 0.01 N 
KIO; was pipetted into the tube; in those made with 
the “combined” alkaline copper reagent this step 
was omitted, since the reagent already contained the 
necessary potassium iodate. In any case, the sides 
of the tube were washed down with 1 or 2 ml of 
2.5-percent potassium iodide without agitation to 
minimize oxidation by air. Immediately following 
the addition of potassium iodide, 1.5 ml of 2 N H,SO, 
was added with shaking, so that the contents of the 
tube were mixed quickly and acidified. When all 
the cuprous oxide was dissolved the excess iodine 
was titrated with 0.005 N sodium thiosulfate and a 
starch indicator. The titer of the sample sub- 
tracted from the titer of the blank containing a like 
quantity of potassium iodate is the net titration. 

his corresponds to the reducing sugar present. 

Prior to the selection of a heating period of 25 
min, the effect of boiling time on the reaction was 
studied. The measurements reported in table 4 
TaBLe 4. Effect of the length of the heating period on reduction 

of 2-mg samples of gentiohose by the Somogyi phosphate 


method 
(Combined copper-iodate solution) 


Milligram of 
gentiobiose 
per milliliter 
of net 
titration 


Net 
titration 
0.005N 
thiosulfate 


Time of 
boiling 


0. 271 
243 
220 
217 
217 
213 
211 


s<seccocecr™ 








show that the amount of reduction is substantially 
constant after 20 min 

The results obtained by the Somogvi-phosphate 
method with gentiobiose and isomaltose are given in 
tables 5 and 6. The amount of reduction was found 
to be roughly proportional to the amount of sugar, 
in the range covered by net titrations from 6 to 13 
ml. Smaller samples gave slightly lower factors and 
larger samples higher factors. For this reason the 
weights of the dextran samples were adjusted to give 
titrations of from 6 to 13 ml, and the same factor was 
used for all samples The results obtained with 
clinical dextrans, using a 25-min heating period, are 
given in table 7. The molecular weights based on the 
gentiobiose standard were obtained by the formula 


M, ; — ——— 
<0.216/342.3 V>< 0.000631 

in which V7, is the number-average molecular weight, 
w is the weight of the sample in milligrams, V is the 
net titration in milliliters of 0.005 N reagent, and 
0.000631 is the reduction equivalent of dextran 
expressed in millimoles per milliliter of reagent. 
Multiplication of the results by 0.939 gives values 
based on isomaltose, and multiplication by 0.969 
gives values based on an empirical factor that is the 
mean of the values found for isomaltose and gentio- 
biose. This empirical factor is in agreement with 
a value found by Commercial Solvents Corp. for 
isomaltose 
TaBLe 5 


Reduction equivalents by the Somogyt phosphate 


reagent determined with gentiobiose 


Net Milligram of 
titration sugar per 
0.005 N milliliter of 
solution titration 


Sample 
gentiobiose 


Separate addition of iodate 


ml 
5. 80 
ht! 
8. 47 
549 
10. 28 


¥. S55 
10. 25 
12. 09 
12.06 
12.06 


13.0 


13. 70 
13. 67 


Mean 216 


Combined copper and iodate 
solution 


From a somewhat similar study [12] the group at C.8.C. selected a heating 
period of 20 min 





| 
| 
| 
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Re duction Somogyt combined 


phosphate-iodate reagent determined with isomaltose 


TABLE 6 equivalents of the 


roppe 


min boiling time 


Net Milligram of 
te titration isomal tose 
altose - 

1.0005 N per milliliter 
solution net titration 


TABLE 7 Vumber-average molecular weights of clinica 


dertrans determined by the Somogyi phosphate method 


gentiobiose standard 
Heating tin 
litration 


0.005 N 
thiosulfate 


Weight of 


Dextrar 
sample 


sample number 


ml 
70 
&. 4 
& 8) 
mi | 
10. 95 
10. ¥2 
11. 2h 


11. 6 


NBs-C8 


Polyglucose,> P 


SOH5SRHZA 


* M.=w/() 0.000631 
» This substance, a polysaccharide prepared by chemical means, is included, 
even though it is not a true dextran 


2.4. Somogyi Carbonate Method [14] 


Reagents and procedure were the same as those 
given in reference 15, except that a 30-min heating 
time was used. 

The procedure was standardized by use of definite 
amounts of gentiobiose. Because the reduction 
equivalents, given in table 8, did not show any defi- 
nite trend with change in concentration, the average 
reduction equivalent was used rather than a separate 
factor for each titration value. That js 1 ml of 0.005 
N thiosulfate was found, on the average, to be equiv- 
alent to 0.214 mg of gentiobiose. Some results 
obtained with dextran samples are given in table 9 
The number-average molecular weights were ob- 
tained from the formula M,=w/(V 0.000625), in 
which w is the weight of the dextran sample in mil- 
ligrams, and V is the net titration in ‘milliliters of 
0.005 N thiosulfate. 





Th 
view 
tions 
proce 
of ref 
with 
el ctl 
“gal 
stand 


TABLE 8. Somogyi carbonate method 
Reduction equivalents of varying amounts of gentiobics 
(30-min heating period 


Milligram of 

gentiobiose 

per milliliter 
10.005 N 
solution 


Net titra 
tion 0.005 N 
thiosulfate 


entiobiose 


BLE Number-average weights of clinical 


rtrans determined hy the Somogyt carbonate method {1 A 


motecular 


Heating time 30 mit 


ritration 
0.005 N 
thiosulfate 


Weight of 


Dextran sample number 
’ f , sam ple 


mag 
S440 4 =. 
M4657A 
S4668 4 


3-155 
B-50074 
1 
R-2 
B-742 
NRC-2 


NRC-3 


NRC-4 


awe 


s46574A 


CSC275 


CSCZl. 


NBS-C8 


NRC-2> 


NRC-2¢ 


* M,=w/(V X0.000625). 
> Reclaimed from sample N RC-? after dissolution in water, precipitation with 
nethanol, dissolution again in water, treated with vegetable char, filtered, and 
he filtrate lyophilized 
Same as NRC-2-A, except no char treatment was used in its preparation 


2.5. Folin-Wu Method [16] 


The method of Folin and Wu was studied with the 
view that it might be useful for routine determina- 
tions in the hands of laboratory technicians. The 
procedure was essentially that described on page 198 
of reference 17. The color measurements were made 
with both Klett-Summerson and Lumetron photo- 
eleetric colorimeters. As with other microreducing 
sugar methods, it is essential that the analyst 

andardize his procedure by measurements of the 





reducing power of a series of standards covering the 
range to be used. Table 10 gives some results 
obtained with anhydrous dextrose and gentiobiose. 


Percentage of transmission of typical standards 


Jor the Folin and Wu method 


Tasie 10 


i] range 400 to 465 my 
of blank is 100 


No. 42 blue filter approximate spectr 
mission 


Gentiobiose 


Dextrose 


Milligram Milligram 
{su rin Transmis f sugar in 
, _ a SUE 
sion " 
100 milliliter ' 100 milliliters 


In the analysis of dextrans, a standard, ordinarily 
gentiobiose, was always included in each set of 
determinations. Thus, the small variation in 
reducing power from one set to another, as well as 
from day to day, is compensated for by the applica- 
tion of suitable corrections. The dextrose equivalent 
or the gentiobiose equivalent for dextran is then 
obtained by calculation, or by graphic means. 


It has been our experience that it is desirable to 
employ standards yielding a transmission reading 
close to that obtained for the dextran under analysis. 
Some values obtained by this method are included 


in table 1. 
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An Interferometer Procedure Applied to the Study of the 
Chemical Durability of Silicates, Enamels, and Metals 


Robert G. Pike and Donald Hubbard 


The chemical durability of quartz, fused silica, obsidian, flint, opal, and glass were 
determined by the interferometer method and the results interpreted in the light of the 


Donnan theory of membrane equilibria. 


the uneven distribution of migratable ions, as shown by the Ag*+Br~ titration 


A more sensitive test of the same phenomena is 


The tests 


used are sensitive enough to differentiate between the different types of silica and may be 
applicable to the study of the chemical weathering characteristics of all natural silicates. 
Results of preliminary experiments on enamels and metals show that these tests may serve 
as indicators of the corrosive resistance of these materials under a wide range of conditions 


l. Introduction 


The results reported in earlier publications have 
demonstrated the utility of the interferometer pro- 
cedure for determining the chemical durability of 
g particularly in the study of the voltage 
anomalies of the glass electrode [1, 2].!_ These results 
suggested that the same procedure could be used to 
determine the relative chemical durability of natural 
minerals and thus throw some further light on the 
probable mechanism and rate of chemical weathering 
in nature. 

Another investigation has also been described 
which shows that there is an uneven distribution of 


classes, 


| tained at a temperature of 80° C for definite periods 


of time. Using the unexposed portion of the sample 
as a reference plane, the attack on the exposed sur- 


| face is measured by determining the displacement of 


migratable ions at a glass-solution interface in accord- | 
ance with the dictates of the Donnan theory of | 


membrane equilibria [3]. 
it seemed desirable to determine if the same principles 
could be applied to some of the natural silicates. 

With these objectives in mind, quartz, fused silica, 
opal, obsidian, flint, and pitchstone were chosen for 
study. The glasses, Pyrex brand chemical ware 7740 

hereafter referred to as Pyrex) and Corning 015, on 

which many tests had previously been made, were 
also included in the study as control and comparative 
samples. 

The results obtained in the first experiments 
suggested that the interferometer procedure might 
have even wider applicability. Therefore, the, 
the method was extended to the study of porcelain 
enamels and metals. 


2. Experimental Procedure 
2.1. Interferometer Method 


When available, specimens about 4 by % in. were 


ground flat on a cast-iron lap with fine alumina pow- 
der and polished on a wax wheel with levigated alum- 
ina. These specimens were treated for chemical 
ducability by the interferometer procedure [1, 2]. 
his consists of immersing the sample one-half its 
depth in a solution of known pH covered with an oil 
‘yer to prevent evaporation of the water, and main- 


Figures in brackets indicate the literature references at the end of this paper. 


In view of these findings | 





the interference fringes with the conventional Pul- 
frich viewing apparatus, using an unfiltered helium 
lamp for illumination. Swelling of the surface is 
plotted in all figures as negative attack. To make 
the results comparable with previous publications 
[1, 2, 3, 6] the exposures were made at 80°+0.2° C 
for 6 hr to Britton-Robinson universal buffer solu- 
tions [4] (hereafter called Britton buffers) ranging 
from pH 2 to pH 11.8, with some exceptions as 
explained later in the text. 


2.2. Hygroscopicity 


The procedure for determining hygroscopicity was 
that reported in previous publications [6, 2], which 
consisted of weighing the water sorbed [5] upon ex- 
posing about 1.5 g of powdered glass (passed through 
a U.S. Standard Sieve No. 140) to the high humidity 
(approximately 98 percent) maintained by a satur- 
ated solution of CaSO,.2H,O. To insure speed in 
obtaining equilibrium and uniformity throughout the 
humidity chamber, the walls were lined with a blotter 
wick, and the atmosphere of the closed system kept 
circulating with a fan. The samples were exposed 
for 1- and 2-hr periods at 25° C and the results re- 
ported as milligrams of water sorbed per cubic centi- 
meter of sample. 


2.3. Uneven Distribution of Migratable Ions 


For the experiments on the uneven distribution of 
migratable ions between the outer solution and the 
silica rich layer of the specimens, samples of approxi- 
mately 10 g were prepared by being crushed and 

assed through a No. 50 and ceaheal on a No. 140 
J. S. Standard Sieve. They were then leached for 
various periods of time at pH 4.1 and 80° C with 
occasional agitation. After rinsing and decanting 
the samples twice, they were treated by agitation for 
5 min in 50 ml of saturated Ag(NH;), Br solution [3]. 
The samples were then separated by decantation 
from the liquid phase. The sample and liquid phases 
were acidified with pM on os H,SO, and titrated 








potentiometrically for excess Ag* remaining in the 
gelatinous silica-rich surface and for the excess Br 
ions appearing in the outer solution [3]. 


2.4. Materials 


Flint samples were cut from the black center | 


portion of a nodule from Travis County, Tex. 
Pitchstone was from Georgetown, Colo. Obsidian 
was a typical natural black glass. These three 
materials were obtained through the courtesy of 
Edward P. Henderson at the United States National 
Museum. Quartz samples were cut from a near 
optical grade crystal. Fused silica was commercial 
optical grade. Opal from 

through Wards Natural Science Establishment. 
Pyrex was chemical glassware 7740. Enameled 
metalware samples were regular commercial grade, 
furnished by the Enamel Metals Section of the 
National Bureau of Standards. Metal specimens 
were commercial samples. 


3. Results and Discussion 
3.1. Silicates 


Nevada was obtained 


The extent of surface alteration (swelling or 
attack) at various pH values has been demonstrated 
with certain glasses {1, 2, 3]. The swelling has been 
rationalized in terms of the Donnan theory of mem- 
brane equilibria [3], and the attack in terms of rate 
of solution of silicates. In the figures the amount 
of swelling is plotted as negative attack. Quartz, 
flint, obsidian, and fused silica show no detectable 
swelling or attack at pH values below 9.5. Above 
this value, the rate of attack increases to \ fringe 
at pH 11.8 for fused silica. The difference between 
the fused silica and crystalline quartz is easily 
demonstrated (fig. 1 data given in table 1). The 
fact that obsidian, a natural glass, falls between 
the two in chemical durability is perhaps surprising 
but may possibly be a result of a high-alumina con- 
tent characteristic of obsidian glass. The attack on 
flint and crystalline quartz was not detectable at 
pH 11.8 after a 6-hr exposure. However, there is a 
detectable attack after 24 hr. 


TABLE I. 


Surface alterations for exposures of 6 hr, 80° C at pH 
Material 


8.2 118 


Fringes) Fringes Pegs Fringes rage nae Fringes: mgicm* mg/cm! 
De ND ND ND ND 4.2 


ND 
ND ND 
ND 


ND 
Hie—s! He-8 | | 
ND ND | ND 
b ag*8 18 


ND 
ND 
ND 
Pitchstone | 

Corning 015 lMat+A 

4y-A ! 
e-8 34 Al MA 


Pyrex 
Opal(Si0;.nH;0) 18 





5% 
NaOH 


A sample of the crystalline quartz was also given 
an exposure of 2 weeks at 80° C in pH 4.1 and showed 
a definite swelling at the contact of the solution oi! 
surface. Longer exposure times could be used if a 
method could be found to prevent diffusion of the 
water through the oil layer. 

In order to show more clearly the relative resist- 
ance to attack in strong alkalies, the samples wer 
also treated for 6 hr in 5 percent NaOH at 80° C and 
the results listed in table (Because of the ver\ 
low resistance of opal it was treated for only 15 min 
and the results extrapolated to 6 hr.) 

The Corning 015, a glass of poor chemical dura- 
bility, is included as a means of comparison with 
previous tests on optical glasses. 

The results obtained on opal are the most inter- 
esting of the group. Swelling of nearly 1 fringe is 
obtained at the lower limit of the tests (pH 2) and 
continues at that rate to about pH 8 and from there 
to pH 11.8 the swelling rapidly changes to an attack 
of 3% fringes. It was noted during the tests that the 
swelling that took place in the solution disappeared 
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Figure 1. Relative chemical durabilities over the range pH 
to pH 11.8 after exposure for 6 hours at 80° C in Britton buffers 


Chemical durability, hygroscopicity, and interfacial uneven distribution of migratable ions for silicate materials 


Excess Ag(NHa)a ions in surface of speci- 
mens leached at pH 4.1, 80° C in— 


Water sorbed in 


1 hr 2hr 0 hr | 6 hr 24 br 72 br 


M-eq/¢ 
».6X10-* 


| M-egig M-eq¢ 


—0.4XxX10" 


M-eqi¢ 
7.9 | <*O.1x1l0O" 


1.0 


<@1 0.9 
<.5 i) 0.9 


\30.6X10+ | 39.3 
| 


166 


5.8 3 


52.3 





* ND, not 5 detectab le. 
><, less t the smallest titrateable increment. 
* 


, at 
4 D, dectable, but not measurable 


‘8; peel slightly less than 4. 


C, surface cut at the + buffer boundary. 
, “, slightly more than 4 


, 





the specimen was exposed to the air. Therefore, | posed to the atmosphere for varying lengths of time. 
\ditional tests were made, the results of which are | From this curve it appears that after about 50 hr the 
hown in figure 2 (data given in table 2), where | amount of swelling approaches a maximum. In 
itack is plotted against the number of hours after | other words, the extent of attack is equal to or greater 
emoval of the samples from the solution. It can be | than the extent of swelling. If the swelling as noted 
een that after a period of 5 hr the initial swelling | in these curves is caused by the uneven distribution 
f more than 1'4 fringes has changed to an appreciable | of migratable ions, the osmotic pressure should be 
ittack, which leads to the conclusion that both solu- | _ ; ; : , 
‘on and swelling had occurred during exposure. In | TU 3+, Piet of time on the ewelling of opal in Briton bulfer 
figure 1 it can be seen that swelling is more rapid | from the buffer and exposure in air P 
han attack for all values below pH 10.2. This is é 
shown more strikingly in figure 3 (data given in 
able 3), where surface alteration (attack, swelling) 
s plotted against time of exposure at pH 4.1. Also oh eo ohe 12 he 26 be 48 hr 
shown are the changes in the surface alterations of 
the same specimens, taken from the buffer and ex- Fringes Fringes Fringes Fringes Fringes ‘Fringes 

3 : . 1 


4 Se l 1+ 8 it 8s +S 448 


Surface alteration after exposure in buffer for 


. . , Hours of 
TABLE 2. Surjace alteration oj the two specimens of opal upon e\ posure Surface alteration after removal from buffer 
exposure to air after immersion in Britton buffer at pH 4.1 in air 
for 24 hours at 80° C 


ace alteration in 
Sample 
number 
0 hr 2 hr 6 hr 9hr 24 hr #8 hr 


Fringes Fringes | ringes Fringes Fringes Fringes 
Ilo Ss us lio A t,—A ytA byt A 
s 


pA . A i \ 


* Symbols explained in table |! * Symbols explained in table |! 
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FiGurRE 2. Surface alteration of two specimens of opal in air 
after exposure to Britton buffer solution at pH 4.1 for 24 hours 
at 80° C. 


SWELLING ATTACK 


SUPPACE ALTERATION, FRINGES 


TE, HOURS 
Ficcre 3. Typical surface alterations of opal with time of exposure in Britton buffer solution at pH 4.1 at 80° C followed by 
exposure in air 


, Surface alteration immed ately after removal from the buffer solution. 
, Surface alterations after removal from the buffer and exposure to air. 
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Figure 4. Typical reduction of swelling of opal and Corning 
015 by increased concentration of electrolyte, exposed for 6 
hours at 80° C. 





TABLE 4. 


Surface alteration of opal and Corning 015 glass in | 


sulfuric-acid solutions of various concentrations for 6 hours 


at 80° C 


Surface alteration 


HySO, 
Corning 
Opal Os 
Percent Fringes Fringes 
01 2+8* 20-8 
LO | Ihe8 219-8 
0.0 =| 1) 8 he-s 
25.0 48 Ds 
™.0 ND, 8C ND 
74.0 ND, 8C ND 
wO0 ND, SC ND 


* Symbols explained in table 1 


reduced by high concentrations of electrolytes in 


accord with the Donnan equation for osmotic pres- 


sure [7]. This effect is shown in figure 4 (data given 
in table 4), where surface alteration of opal and 


Corning 015 is plotted against concentration of | 


H,SO,. It can be seen that at 0.1 percent H,SO, a 
swelling of greater than 2 fringes is obtained, and at 
65-percent H,SO, the amount of swelling has been 
reduced until it is no longer detectable. It is real- 
ized that considering 65-percent H,SO, as a high 
concentration of electrolyte in terms of the Donnan 
theory of membrane equilibria is questionable and 
that any one of a number of other explanations could 
suffice [3]. 

The hy opicity of the samples is plotted in 

gure 5 (data given in table 1) as water sorbed in 
mg/cm? against time in hours. Note that the relative 
“samapeen of Corning 015 and opal in this figure have 
een reversed from that in figure 1, indicating that 
the surface alteration (swelling) of the opal results 
from some factor other than its hygroscopic property. 
In other words, hygroscopicity is not a reliable 
indicator of chemical durability. 

In figure 6 (data given in table 1) the excess 
Ag(NH;).* ions in milliequivalents per gram (here- 
after designated as M-cq/g) appearing in the surface 
is plotted against time of leaching at pH 4.1 at 80° C. 
The relative positions of these curves are compatible 
with those i wie alteration shown in figure 1. 
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Figure 5. Hygroscopicity (water sorbed) for quartz, fused 


SiO, obsidian, opal, and an electrode glass (Corning 015) 
after exposure in an atmosphere of 98-percent relative humidity. 
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Figure 6. Uneven distribution of migratable ions at the 
interface between the surface layer and the outer solution after 
various periods of leaching in Britton buffer at pH 4.1 at 
80° C using Ag (NH;)f and Br- as the indicator ions. 


This uneven distribution of the migratable ions is 
caused by the nonmigratable negative charged sili- 
cate ions acting as their own semipermeable mem- 
brane, and thus causing a definite osmotic pressure 
across the opal-solution interface. 

It is well recognized that chemical attack Ewe a 
considerable part in the weathering of rocks and 
minerals. These tests strongly poioat out the fact, 
not always realized, that the weaker acids and bases 
or even neutral salts may sometimes be a much more 
destructive force than concentrated solutions. These 
studies also indicate that another phenomenon, os- 
mosis, plays a very important and persistent role * 
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he weathering of numerous rocks and minerals. 
rhe interferometer procedure might we'l be used to | 
study the relative resistance of many of the rock- 
forming minerals to acids and bases over a much 
wider time and concentration range than was done 
here. 

3.2. Alkali Aggregate Reaction 


High-alkali cements are known to react with cer- 
tain aggregates, particularly opal, used in concrete 
[8]. This reaction causes a swelling and subsequent 
deterioration of the concrete structure. Hansen [9] 
has postulated that the cement paste might act as a 
semipermeable membrane, which would allow the 
alkalies to migrate through the concrete and react 
with the opal. The reaction product would not be 
capable of diffusing through the hardened cement 
paste. Thus osmotic pressures would be built up, 
which could cause disintegration of the concrete. A 
modification of this hypothesis as advanced by 
Parsons and Insley [10] suggested that a reaction 
takes place at the walls of the pores in the concrete 
to form a membrane impermeable to further migra- | Ficure 7. 


own semipermeable membrane. However, the inter- 
ferometer shows that expansion of the opal itself 
occurs in the acid and not in the alkaline range. Any 
expansion in the alkaline range must be because of 
some reaction products or to osmotic pressures set 
up through some membrane surrounding the irdivid- 
ual opal grain rather than to any Donnan membrane 
equilibria in the opal itself. Another interpretation 
of this effect has been advanced by Vivian [11], who 
believes that alkalies react with opal to form a gel 
that absorbs water, causing swelling of the gel and 
subsequent failure of the concrete. 


3.3. Porcelain Enamels 


WATER SORBED, mg/m? 


In figure 7 (data given in table 5) is plotted the 
chemical durability of three ponent porcelain 
enamels, 9 class A acid resistant, a class D acid re 
sistant, and a hot-water tank enamel of a type used 
in the manufacture of the so-called glass-lined tanks. 
(Hereafter referred to as gless-lined tank enamel.) 
{12}. The difference in durability in both the acid 
and alkaline range can be observed from these curves. 
The hygroscopicity of these enamels is shown in 
figure 8 (data given in table 6) and is seen to be | 


Ficure 8. 


—— TF) + oo 
| 


ao} 
20, 
} 


r 


° 


. 


° 


== 


. 


b- 


SURFACE ALTERATION, FRINGES 
ATTACK 


SWELLING 


be 





- o——__@—+@- + <}—-—_@ 


CLASS 0° | 


Jocassal 


1 


Li 
ss 





| 


GLASS -/LINED TANK 





, Glasslined tank; @, class A, and +, class D acid-resistent enamels 


emt : Chemical durability over the range from pH 2 to 

tion of alkali silicates. The authors have considered | 11.8 of three types of porcelain enamels, after exposure for 
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the possibility that the opel grains might act as their 6 hours at 80° C in Britton buffers. 
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| Corrung O15 








Hygroscopicity (water sorbed) of three typical 


greater than that of Pyrex but much less than that porcelain enamels compared with two commercial glasses, 
of Corning 015. | after exposure to a relative humidity of 98 percent. 


TaBLe 5. Chemical durability of three types of porcelain enamel 


Surface alteration at pH— 


Enamel 
41 | 6.0 
} 


Fringes | Fringes | Fringes 
Acid resistant, class A ND | ND 
Glass-lined tank K+8 | ND ND 


Acid resistant, class D 3 368 le+8 ND 


* Sym bols explained table | 


91 


Fringes 


8.2 9.0 11.0 


Fringes Fringes Fringes 
ND ND I MA 
ND DA 


Mio A 








Comparison of the hygros« opicities jor three types of 
porcelain enamels and two commercial glasses 


TABLE 6 


Water sorbed in 


Specimet 
b he 
Porcelain enamels n N" mgom 
Class A, acid resistant l Is 
(jlass-lined tank 21 27 
Class D, acid resistant a4 y 
(i lasses 
Corning O15 21 v1 
Pyrex 7740 s 4 


The chemical durability of the class A enamel at 
high concentration of acids is shown in figure 9 
(data given in tables 5 and 7), and the repression of 
swelling of the tank enamel at high concentration 
of electrolytes is shown in figure 10 (data given in 
table 8). The uneven distribution of the migratable 
ions for the class D enamel gave the low value of 
less than 0.2” 10~* M-eq g with no leach, and the 
very high value of 415 10~* M-eq/g after 6 hr of 
leaching at pH 4.1. This value is 46 times greater 
than for the glass-lined tank enamel and is a very 
striking example of the sensitivity of this test. 


TasLe 7 Chemical durability of two commercial porcelain 
enamels in aqueous solutions of H.s0, and CH,COOH 


Surface alteration 
Solution 


Class A, acid Glass-lined 
resistant tank 


HeSOy CHyCOOH HO, CH»COOH HesOy CH;»COOH 


Percent pH pH Fringes Fringes Fringes Fringes 
01 1 27 *ND ND 

Lo o.75 2 ND . 
10 0 ND | ‘ 
2.0 0. 62 ND 

00 1.49 1.3 ND ND ND ND 
75.0 245 0.¢ ND ND ND ND 
“oO &1 ND ND 

“0 22 ND ND 


* Symbol: explained in table 1 

> The iron backing of the glass-lined tank enamel was severely attacked at 
0.1 to 25 percent HeSO. and CH»COOH, H: was evolved, which chipped and 
shattered the enamel, thus destroying the specimen 
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Figure 10. Comparison of the uneven distribution of migratabl: 
ions at the enamel-solution interface with the chemical dura 
bility of a porcelain enamel used in glass-lined tanks, after 
exposures for 6 hours at 80° C. 


TABLE 8. Uneven distribution of Ag( NH, > at the enamel- 
solution interface after leaching in various solutions for 6 
hours at 80° C 


Enamel Treatment annie . 
pl Solution M-eqg¢ 
None *<0.2xk10" 
S11 WT, HeSOy 10 
. — 20 Britton buffer 16 
Glass-lined tank 23 CH;COOH 20 
4.1 Britton buffer ” 
11.8 do 4 
: . j None 0.2 

Chass D acid resistant { 4.1 Britton buffer 415 


* Symbols explained in table 1 


Because the composition of the enamels was not 
known in detail, it is difficult to discuss the mecha- 
nism or nature of the attack indicated by the inter- 
ferometer data. The evidence, however, that por- 
celain enamel is subject to alkaline attack correlates 
well with field tests that have been made on enamel- 
ware [13]. It appears that the nature of the attack 
on the enamels is analagous to that on commercial 
glasses. It would seem, then, that some of the 
problems of the chemical durability of different 
enamels could be advantageously studied by these 
techniques. 


3.4. Metals 


Although the interferometer has been used widely 
in science and industry, the procedure for measuring 
the chemical durability has not been applied to the 
study of the corrosive resistance of metals and alloys 
to acids and bases of all concentrations. The chem- 
eal durability of cold-rolled steel, aluminum, copper! 
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RE 11 Chemical durability of various metals and alloys 
after exposure to Britton buffers for 6 hours at 80° C. 


Chemical durability of metals exposed to Britton 
buffers for 6 hours at 80° C 


TABLE ¥Y. 


Surface alteration at pH 


2.0 41 10.2 11.0 11.8 


Fringes 
ND 


Fringes Fringes Fringes 
ROA ND ND 


Fringes 


Fringes Fringes 
* 0A ND 


Cold-rolled 
steel 

Aluminum IA ND ND ND ND. SC, BA 
black 

ND 1+A 

ND le A 

ND 


ND 
ND ND 


~O A 


ND 
4108 
ND 
ND 


2A 
1A 
ND 
ND 


ND 
Ds 
ND 
ND 


Copper 
Brass 
Silver hiotA 
Platinum ND 


WA 
4A 


Symbols explained in table 1 


brass, silver, and platinum is shown in figure 11 
data given in table 9). Because of the rapid reac- 
tion of some of the samples, exposures as short as 

min were required, and the results were extra- 
polated to 6 hr in order to put all the data on a 
comparable basis. The swelling of brass and the 
slight attack on silver might be more convincingly 
demonstrated by increasing the exposure time to 
24 hr or more. When aluminum was treated at 
pH 10.2, a hard black almost nonreflecting coating 
was formed, although no dimensional change could 
be detected by the interferometer. In fact, in many 
cases there was a color change in the treated area 
but with no detectable dimensional change. Hass 
i4] states that the oxide film on aluminum may be 
as thin as 20 A. Because the lower limit of resolu- 
tion of the interferometer is about 30 A, no swelling 
could be observed. It is also interesting to note 
that the points at which surface alteration becomes 
detectable correspond very closely to the values 
determined by Wells [15], namely, pH 3.95 and 
pH 10.94, where aluminum hydroxide becomes sol- 


| 


uble in the acid and alkaline regions. In several 
samples where a slight attack or dull matt swollen 
surface had resulted, repeated gentle polishing on 
4/0 emery paper increased the observable attack in 
the exposed area to a definite maximum. This in- 
dicates that chemical attack had progressed to a 
considerable depth below the surface. In all cases 
the maximum attack noted after the emery-paper 
treatment was reported as the chemical durability. 

The results on all metals tested show good cor- 
rosive resistance from pH 6 to pH 10. There is, 
however, a very definite attack on all the metal 
surfaces except platinum at pH values below this 
range. Because of the sensitivity of the interfer- 
ometer measurements, the method offers a means of 
studving the corrosive resistance of metals and 
allovs under very diverse conditions. 
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‘ransitionsand Phases of Polytetrafluoroethylene (Teflon) ' 
Charles E. Weir 


Transitions in polytetrafluoroethylene 
and temperature. 
were studied. 


sure 
tion 
sure 
with increasing temperature, 


(Teflon 
A high-pressure transition (5,000 atm) 
These are believed to result from three polymorphic forms. 
of the high-pressure transition decreases and that of the low-pressure transition increases 
A triple point is indicated at approximately 70° ¢ 


have been studied as a function of pres- 
and a lower pressure transi- 
The pres- 


‘, and a third 


equilibrium line is found to proceed to higher pressures with increasing temperature from 


this point. 
the triple point indeterminate 


1. Introduction 
In the course of a continuing study of pressure- 
iume-temperature relationships in high-polymeric 
aiterials measurements have been made on polytet- 
ifluoroethylene (Teflon). This report describes the 
essure-temperature dependence of the transitions 
Teflon. These transitions show a surprisingly 
omplex behavior that is unique among all the poly- 
mers studied to date. 

The room-temperature transition was discovered 
by Rigby and Bunn [1] * and independently by the 
author [2] at about the same time. The transition 
that occurs at about 20° C is a reversible first-order 
iransition producing reversible changes in the X-ray 
diffraction pattern. The changes in the X-ray pat- 
tern have been interpreted as arising from specific 
changes in the crystal structure of the Teflon [1]. 
Quinn, Roberts, and Work [3] studied the transition 
extensively by volumetric methods and observed that 
the transition was complex, the complexity consisting 
of a double break in the volume-temperature curve 
associated with some degree of hysteresis. They 
evaluated the total volume change (that is, the sum 
of the two observed discontinuities) as 1.23 percent. 
Pierce, Bryant, and Whitney [4] have rece oa studied 
ithe X-ray pattern of Teflon over a wide temperature 
range. 

Bridgman [5] and the author [6] have reported 
first-order transitions in Teflon under high pressures 
at room temperature. Bridgman reported the tran- 
sition at 6,500 kg/em* as involving a volume change 
of 2.26 percent, whereas the author found the transi- 
tion centering about 5,500 atm at 25.5° C, with 
un associated volume change of 2.35 percent. Dif- 
ferences between these data may well be due toevari- 
ations in samples. 

Finally, Furukawa [7] has measured the specific- 
heat-temperature properties of Teflon. His data 
show the double discontinuity referred to, and he 
evaluated the total latent heat of the two breaks as 
approximately 8 j/g. It is immediately apparent 
that this latent heat is not sufficiently large to fit the 
assumption that the high-pressure and room-tem- 
perature transitions are related. The 2:1 ratio of 
volume changes now tends to confirm this conclusion. 

As a result of the present studies it is concluded 

This paper is a preliminary summary of work on a project sponsored by the 


fice of Naval Researe 
? Figures in brackets refer to lite -rature references at the end of this report 
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Experimental errors involved render tests of all thermodynamic conditions at 


| that Teflon undergoes two separate transitions, one 
at high and one at lower pressures at temperatures 
| 
' 


below 70° C. The phase diagram for these transi- 
tions indicates at least three forms for the crystalline 
regions of Teflon. The evidence points to the fact 
that these are most probably polymorphic forms. 


2. Experimental Method 


The apparatus and experimental procedure have 
been described in detail previously [8] and will not be 
discussed fully here. Briefly, an experiment consists 
in forcing a leakproof piston into the bore of a heavy- 
walled cylinder containing the specimen and the man- 
ganin pressure gage immersed in a suitable pressure 
transmitting medium, and taking pressure-volume 
readings. For locating a transition a graph of pres- 
sure Versus piston position is constructed, the transi- 
tion appearing as a discontinuity in this curve. 

In these experiments the hydraulic press and 
pressure vessel were enclosed in a thermostatted air 
bath that could be maintained at elevated or reduced 
temperatures. Temperatures were measured with a 
calibrated mercury thermometer set in a well in a 
massive brass block fastened to and in intimate 
thermal contact with the pressure vessel. A maxi- 
mum variation of £0.35 deg C was noted on this 
thermometer at any temperature. However, the 
enormous heat capacity of the pressure vessel 
probably reduced the variation in temperature of 
the specimen to a considerably smaller value. 

Inasmuch as the pressure gage was subjected to 
various temperatures in these experiments, it was 
calibrated at most temperatures used. Fixed points 
for calibration were taken from the International 
Critical Tables [9] but originate from data obtained 
by Bridgman. The transitions used for calibration 
were: Water-ice VI at lower temperatures; freezing 
of chloroform and CCl, I to II at intermediate 
temperatures; and urea I to IIT at higher tempera- 
tures. Small inconsistencies noted in these calibra- 
tions were not studied, the variations being of such 
small magnitude that any error resulting was negli- 
gible in these experiments. However, subjection of 
the gage to the elevated temperatures resulted in 
limiting the measurements to a maximum tempera- 
ture of approximately 80° C. After continued 
exposure to this temperature, marked drifting of the 
resistance of the pressure gage was observed. This 
drift is probably due to the action of the confining 


95 





liquid on the shellac-treated gage coil. A study of 


the behavior of this or other materials to higher | 


temperatures would require considerable modification 
of the apparatus to permit the pressure gage to 
remain at room temperature. 

All data in these experiments were measured 
isothermally. Isobaric measurements were 
tempted, but the large heat capacity of the apparatus 
effectively excludes interpretation of this type of 
data, the actual temperature of the specimen being 
unknown. 

A single specimen of Teflon was used in all measure- 
ments. It was obtained from E. I. du Pont de 
Nemours & Co., Inc., and had a density of 2.219 
g/cm’ at 20° C. The initial volume of the specimen 
was approximately 21 cm® at 20° C, initial volumes 
at other temperatures being calculated from the 
data of Quinn, Roberts, and Work [3], who studied a 
specimen cut from the same sample. 


3. Results and Discussion 


Typical compression curves of Teflon are shown in 
figure 1. These measurements were made at the 
temperatures indicated, but the curves are shown 
displaced along the vertical or compression axis to 
reduce overlapping. The numerical compression 
data for these curves will be given in a subsequent 
report because they are not of primary interest — 

It is seen in figure 1 that two transitions exist at 
intermediate temperatures, a high-pressure transition 
involving approximately a 2-percent change in vol- 
ume and varying with temperature around 5,000 
atm and a low-pressure transition involving a 1-per- 
cent volume change at somewhat lower pressure. 
This low-pressure transition occurs even at 38° C 
and at lower temperatures, but at pressures below 
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1,000 atm which is the lower limit for these data 
At 81.0° C, however, no corresponding transitio; 
can be found. 

The following points of interest are to be noted i) 
figure 1: The rounded corners at the transition pres 
sure, combined with the existing curvature arisin; 
from pressure dependence of compressibility, rende 
evaluation of AV for the transition rather uncertain 
There is little doubt, however, that AV for the high 
pressure transition increases, whereas AV for thy 
low-pressure transition decreases with increasing 
temperature. 

Both transitions occur over a small pressure in 
terval rather than at a fixed pressure. Expecte: 


small differences in structure of crystalline regions 


in Teflon would cause this effect as well as the round 
ing of the corners at the transition. 

The behavior at the high-pressure transition with 
increasing temperature is interesting. At low tem- 
peratures this transition is “smeared out’ over « 
rather large pressure interval, as seen in figure | 
With increasing temperature it becomes noticeab)) 
sharper, and at temperatures above 60° C the tran- 


_ sition may be observed easily, since the piston may 





be withdrawn 0.040 in., with no change in pressur 
at the transition point. This type of behavior has 
been referred to many times by Bridgman [10]. 

All data have been obtained on decreasing pressure 
although there is no evidence of a region of indiffer 
ence [10] in these transitions. On increasing pres- 
sure, however, with only the low-pressure phase 
present, considerable excess pressure must be applied 
to produce initiation of the transition. 


Taste l. Experimental transition data on Teflon 
Transition 
Ito ll IL to Ul Ito ll 
P t P t / t 
atm Cc atm 4 4 atm Cc 
520 2.7 5, 550 &.0 4.540 75.0 
590 31.0 5, 650 9.2 5, 050 81.0 
1,060 38.5 5, 40 16.9 4, 040 S1L.0 
2, 130 50.5 5, 350 21.0 
2,470 57.0 5, 260 $2.6 
3, 200 4.0 5, 080 38.6 
4,000 6a 5 4, 910 3.5 
4.740 57.0 
4, 560 6.0 
4, 30 4.0 


The experimental data obtained on these transi- 
tions are shown in table 1 and figure 2. The several 
regions are denoted as corresponding to forms I, Il, 
and III, and the P,7 data experimentally obtained 
on the transitions are those given in table 1. 

The equilibrium line II to III is well defined and 
linear. The least squares equation of this line is 
found to be 

P=—17.5t+5761 (1 


for pressures in atmospheres and temperatures in 
degrees C. The negative slope requires in a transi- 
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Figure 2 Proposed phase diagram for 
on from IT to TIL, in which the volume decreases, 
hat heat be absorbed. This type of transition, 
though rare under l-atm = pressure, has been 
ibserved many times under elevated pressures [10] 
With increasing temperature AV increases rather 
owlv along the II to IIT line, and from figure 1 it 
s obvious that IIL is considerably less compressible 
han Il. This difference in compressibilities has 
ot been evaluated quantitatively. The I to Il 
juilibrium line is less well established but is cer- 
ainly nonlinear. Here AY is much smaller and 

dee rapidly with increasing temperature 

Due to the in AY with temperature, the 

point shown at 68.5° C ts subject to some uncertainty 

\t low pressures (that is, below 1,000 atm) the 
yperimental points are also subject to some uncer- 
unty, arising from the fact that the apparatus 
whaves indifferently well at pressures below 700 
im. That the course of the curve is essentially 


creases 


deci ease 


orrect at low pressures, however, is shown by the 
agreement between the latent heat of 1.54 cal/g 
calculated from the graphical slope at zero pressure 
ind the observed value of 1.91 cal/g measured by 
Furukawa [7]. The latter value contains the total 
atent heat of the double transition, as mentioned 
previously. 

No evidence was found for any transition line 
originating at the second discontinuity noted by 
Quinn, Roberts, and Work [3] at 30° C. The low 
pressures required and the small magnitude of the 
volume change involved would effectively obscure 
apy such observations with the present apparatus. 
Inasmuch as X-ray studies have apperently disclosed 








no phenomena at 30° C [1, 4], 
Teflon in the vicinity of 30° C 
study. 

In figure 2 a triple point is shown at approxi- 
mately 70°C. That such a point exists is strongly 
indicated by the abrupt change in slope experienced 
by the I] to ILL transition line in this region. For 
this reason, it is believed that a true triple point 
exists and that the phases I, Il, and TIT are truly 
polymorphic. The alternative choice that the tran- 
sitions occur in incoherent phases seems rather 
unlikely. It appears from the measured data that 
the requirement that A//=0 for a closed path en- 
closing the triple pom| does not appear to be ful- 
filled. However, because the AV velues cannot be 
accurately determined and Al for the IT to Il 
transition becomes extremely small near this point, 
the apparent discrepancy may not be real 

The diagram shown i figure 2 
more than a portion of the behavior of Teflon 
X-ray studies have shown melting of 
crystals at approximately 330° C [1]. At least one 
additional transition line must therefore appear 
at higher temperatures [11]. In this connection, it 
is to be noted that the extrapolation of the IL to TI 
transition line cuts the temperature axis at 329° C 


the behavior of 
requires further 


cappnot represent 


because 


Inasmuch as the details of this diagram cannot be 
further clarified with the present apparatus, it is 
presented with full recognition of its lack of complete- 
ness and the fact that the proposed phases should 
be verified by X-ray studies 
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Hydrodynamic Effects of Gales on Lake Erie 


Garbis H. Keulegan 


The coefficients of wind stress end the sea roughness ere derived from records of water 
levels and wind intensities relating to gales passing over Lake Erie during the past fifty vears 
The reduction of wind velocities meesured at various elevations in four cities on the lake 


shore to a single “effective wind velocity’’ over the lake is explained 
nation of the wind tides agrees well with the observations 


\ theoretice! determi- 
The coefficient of wind stress 


derived from the observations is nearly the same as that obtained by Neumann from obser- 


vations in the Gulf of Bothnia. 


Both the coefficient of wind stress and the sea roughness 


are found to decreese with increasing wind velocities 


l. Introduction 


For a wind blowing steadily with a uniform in- 
tensity over a lake, the relationship between the 
wind intensity and the frictional stress induced at 
the water surface may be found by considering the 
wind tides produced at the ends of the lake. Wind 
tides are those displacements of the water level from 
the initial undisturbed level that are due solely to 
the action of the wind. The relationship involves a 
characteristic quantity called the coefficient of wind 
stress. There will be an element of uncertainty in 
the coefficient derived if the wind stress correspond- 
ing to a wind of constant intensity is a function of 
distance, or fetch, and only the wind tides at the 
lake ends are known. In such cases the analysis 
vields average values of the wind stresses, and if it 
be true that the stress is also a function of fetch, the 
derived coefficient of wind stress will be different for 
lakes of different lengths 

Another source of uncertainty in the values of the 
coefficient of wind stress is in the derivation of the 
formula which relates the wind stress to the wind 
tides. In working out the formula, due attention 
must be given to the nonuniformity of depth of water 
in the lake and of width of the lake surface, and also 
to the magnitude of the friction at the bottom of the 
lake. The turbulence pattern in the waters of a lake 
under the action of wind is at the present time far 
from being depicted. The partition of velocities in 
the drift current at the upper layers and in the re- 
turning gravity current in the lower layers, and the 
magnitudes of the mixing lengths, are questions which 
have not been studied. 
nature prevents one from appraising accurately the 
the frictional force at the lake bottom. Besides, if 
the boundary of a lake exhibits very marked irregu- 
larities and there are also numerous islands distrib- 
uted in the regions close to the lake ends, the 
derivation of the exact wind tide formula is very 
difficult and in some cases may be impossible. These 
difficulties are present to some extent in the case of 
Lake Erie. 

A second problem in connection with the action of 
wind on lakes is the relationship between the effective 
roughness of the sea and the velocity of the wind. 
Meteorologists regard the Prandtl law of velocities 
as representing satisfactorily the state of movement 
of the wind in the layers of air above the water 


Lack of information of this | 


surface. This law requires a roughness length which 
is called the roughness of the sea. Once the coeffi- 
cient of wind stress based on the wind velocities at 
points of known distance from the surface of the 
water is known for a given lake, the sea roughness is 
readily obtained from the value of the coefficient. 
Since sea roughness is an inferred quantity and is a 
function of stress, the likelihood exists that the sea 
roughness is a function of fetch. If the determination 
of roughness is made on the basis of average stresses 
at the lake surface, the value of roughness may vary 
for a given wind velocity for lakes of different lengths. 

It is proposed in the present work to obtain the 
values of the coefficient of wind stress and of the sea 
roughness for Lake Erie. The determinations are 
brought about by considering the maximum wind 
tides observed during the severe westerly gales pass- 
ing over Lake Erie in the past half century. 


2. Sources of Original Data 


The hydrographic charts of the United States Lake 
Survey showing the surface changes of Lake Erie 
during gales constitute the chief source of material 
for the present investigation. These charts show 
the water stages at the extremities of the lake, and 
the values of the hourly movement of wind and its 
direction for a few stations on the southern coast. 
Two examples of such charts are given by Freeman 
1}. 

The stations for which charts of water stages are 
given have been changed from time to time. Buffalo 
and Amherstburg, or Gibraltar, are the localities 
near the ends of the lake for which the water stages 
are given for westerlies during the earlier years. 
For easterlies the corresponding localities are Toledo 
and Buffalo. In these charts compiete wind data 
are shown for Toledo, Cleveland, Erie, and Buffalo. 
In the charts of later years, water stages are shown 
for Toledo and Buffalo. However, the representa- 
tion of wind data has not always been complete, 
and in such instances the missing information has 
been received from the U.S. Weather Bureau. 

It was desired to handle the analysis on the basis 
of the wind tides observed at Toledo and Buffalo. 
The data for Toledo, where missing, were deduced 
from the data for Gibraltar and Amherstburg. This 


is feasible as in general a correlation exists for the 
The desired correla- 


wind tides of nearby localities. 








tion factor could be determined, for in some of the 
charts the simultaneous values of the water stages at 
Toledo and Amherstburg are given. In this de- 
termination use was made also of Blunt’s work [2]. 

In the determinations of the coefficients of wind 
stress, the observed wind tides must be corrected for 
barometric pressure changes. Data needed are the 
hourly values of the barometric pressures at Toledo 
and Buffalo for the periods of the gales considered 
For these the U. S. Weather Bureau has kindly 
allowed us access to the proper records 


3. Determination of Effective Winds 


3.1. City Wind Velocities at a Common Elevation 


The evaluation of the coefficient of wind stress will 
be made with respect to wind velocities at a height of 
25 ft above the water surface of the lake. As these 
velocities were not observed and the only available 
velocities are those measured over the cities of Toledo, 
Cleveland, Erie, and Buffalo, the lake wind velocities 
will be deduced from the city velocities. The re- 
duction to the desired values will be made in three 
steps 

As the anemometer heights during a given period 
have been different in the cities mentioned, or have 
been changed in a given city from time to time, the 
first step of the analysis is to reduce all the measured 
wind velocities to the values appropriate to a common 
elevation. The elevation selected was 165 ft above 
the ground. 

The region of air above the ground may be con- 
ceived to consist of two lavers, the boundary laver 
below and the laver of frictional influence above. 
The boundary layer is relatively thin, and the move- 
ment of air in it is controlled entirely by frictional 
stress as horizontal pressure gradient and deflecting 
forces may be neglected. In a vertical traverse the 


frictional stress is constant, and the direction of the. 


stress vector coincides with the direction of the wind. 
For this layer, Prandtl [3] gives the law of velocities 


5 log -» (1) 


where 2 is the elevation of a point at which the velo- 
city is u, u, is the shear velocity, and ¢ is the effective 
roughness. Denoting the frictional stress at the 
ground by ro and the density of air by p,, 


Us =~ To/ Pa (2) 
When the actual roughness is large, the reference 
level from which to measure the height 2 is a matter 
of importance. It may be that in cities with an 
average height of buildings of 2A the reference level 
can be taken at the distance AK from the ground. 
Then, following Prandtl, the effective roughness e, 
of a city may be taken as 


€.=2K/30. (3) 


| 





a ~ - | 


- 


0.94 






. —f-- 


uw (msec) 





u, m/sec 
n 








0.4 0.0 o4 oe 1.2 1.6 2.0 
fog z, z inmeters 


Relation of wind velocity, u. to he ight over the city of 
, fram measurements by Siotani and Yamamoto 


Figure | 
Tok yo, 


We may now consider, to test the utility of th 
rules given above, a known case of measured veloci- 
ties over a large city. In the autumn of 1946 
Siotani and Yamamoto [4] obtained the air velocities 
over the central part of Tokyo at elevations of 60 
55, 45, 35, and 25 m above the ground, using hot- 
wire anemometers on the tower at the Central 
Meteorological Observatory. The neighboring 
buildings being 10 to 20 m in height, A in this case 
Was 4.0m, 

in figure 1 are reproduced some of the velocity 
values from the original curves of these authors 
In the plot 2 is the height in meters measured from 
the reference level of 7.5 m above the ground. The 
straight lines drawn through the points meet at a 
common point on the abscissa axis, log «© 0.3 
Thus the effective city roughness is e&=0.5 m, a 
result satisfying the rule implied by eq 3. Denoting 
by uw, the value of the velocities to be read from 
the straight lines of the figure for 1, the shear 
velocity may be computed as u, (uv) /(5.75 log 
«). The shear velocities thus computed are also 
entered in figure 1. Representing the data in 
figure 1 in dimensionless form, the results are the 
three points in figure 2. The distribution of the 
points confirms the relation given by eq 1, and this 
fact points to the possibility of applying the rule to 
cities in general. 

Estimating the average height of buildings in 
American cities as 30 ft, the effective roughness ma) 
be taken to be 1 ft. Measuring the elevation of 
points 2 from a reference level 15 ft above the 
ground, denoting the position of an anemometet 
by 2», the measured velocity by u,, the standard 
height by 2 (= 150 ft), the velocity at this height by 
Uy», the roughness of the city by & (=1 ft), the 
formula to reduce the anemometer velocities to the 
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Ficture 2 Relation of wind velocities measured over the city 


of Tokyo to the Prandtl law 


velocities at standard heights is 
Uo 


“ log log (4) 
The curve in figure 3 is established through this 
relation and is used to reduce the anemometer 
velocities to the velocities corresponding to the 
standard height 2, (=150 ft). For the purpose of 
comparison a second curve is shown, evaluated with 
«2 ft. The disparity between the curves being 
small, it is seen that errors arising from the uncer- 
tainties in the exact value of the city roughness 
are not very appreciable. 
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Ficure 3. Ratio of wind velocity at a standard height of 150 | 


feet above a city, Ux, to the wind velocity indicated by an 
anemometer, tig, as a function of the height of the anemometer, 
», and the roughness of the city, e. 
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3.2. Wind Velocities Over the Lake Surface 


The second step is to infer the wind velocities over 
the lake surface from the wind velocities u.. corre- 
sponding to the standard height 2, (= 150 ft) over 
the city. To establish the desired values resort is 
to be made to the characteristic relationships in- 
volving the various quantities of the layer of fric- 
tional influence. Regarding the hyvdrodynamies of 
the laver, the analysis of Rossby [5] proves to be 
particularly fruitful. 

The quantities that will be of value here are the 
height, 77, of the boundary layer, the wind velocity, 
{’y, at the upper surface of the boundary laver, the 
gradient wind, (’,, and the angle @, between the 
direction of the gradient wind and the direction of 


the wind in the boundary laver. The Rossby 
relations are 
{ " { P ( cos o sin ®, ), (5) 
y2 
11-0.54k 0.065, (6) 


© 

sin @,, k 
/ 
and 


log { e/ fe 1.604 cot @ log sin @, + 1.441 (7) 
Here f represents the quantity f= 22 sin A, 2 being the 
angular speed of the earth’s rotation, and \, the lati- 
tude. The roughness of the surface below is «. The 
importance of the effective roughness, of the ground 
or of the sea, for the quantities pertaining to the lower 
lavers of the atmosphere can easily be seen, 

Once it is assumed that the magnitude and the 
direction of the gradient wind are the same over the 
lake water and over the coastal regions next to the 
lake, eq 5, 6, and 7 provide the means by which a 
relationship may be established connecting the ve- 
locity of the wind over the surface of the lake with the 
velocity of the wind over a city on the borders of the 
lake. 

For a second definition of terms the diagram im 
figure 4 may be consulted. Assuming a gradient 
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wind constant both in intensity and in direction 
blowing in the higher regions over the lake and over 
the city, let 4/7, denote the height of the boundary 
laver over the lake and //, that over the city. Let 
¢, denote the effective sea roughness and «& the rough- 
ness of the city. Let «, denote the wind velocity at 
the upper surface of the boundary layer over the 
lake and w, that over the city. Let u» denote the 
wind velocity over the city at the standard height 
>, (= 150 ft) and u,, that over the lake at the standard 
height 2, (=25 ft). Let Ow be the angle between the 
lireection of the wind over the city, and the direction 
of the gradient win 1, and @, the corresponding angle 
over the lake. Then, according to the Prandtl law 
of velocities given by eq |, the formula which reduces 
the wind velocities over the city to those over the 


lake is 


log ~ log ~~ 

oe My €; €, 
. . (S) 

M22 (Ue log Mh log H; 

€) €, 


the computations being made through eq 5, 6, and 
7. The results are shown in the form of curves in 
figure 5, with the ratio w,,/u. plotted as function of 
u» for various sea roughnesses. It is seen that the 
ratio u))/%» for a given wind velocity depends to an 
appreciable degree on the sea roughness. In using 
the graphs for the desired reductions of wind veloci- 
ties for a given storm, it Was necessary to assume a 
value for the sea roughness. The value selected was 
¢=0.3 em. Evaluations from the data of wind tides 
given later in this paper confirm this chosen value as 
an average. 

It will be remembered that the curves of figure 5 
are based on a city roughness of value e=1 ft 
(30.5 em). Accordingly, the figures show that if the 
sea roughness be ¢=0.3 em, the wind velocity ob- 


served at the lake surface at a height 25 ft above the / 


water will be about 1.12 times the city wind velocity 





0.9 
° 10 zo 30 40 so 60 7° 


Us, mph 


Figure 5. Relation between winds at a standard height of 25 
feet over a lake, uy, to winds at a standard height of 150 feet 
over a city, Ux, for various sea roughnesses. 
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150+15) ft above th 
ground for the higher velocities of the city wind 
It will be found, considering the data of figure 3 alse 
that if the city wind velocity be observed at a heigh: 
315 ft above the ground, it gives the wind velocit 


observed at a point 165 | 


over the lake. This is in agreement with the finding 
of Kadel [6] in an investigation of wind observations 
at Buffalo. 

The roughness of the sea being much smaller than 
the roughness of a city, there is a change in the wind 
direction in passing from a city over a lake surfac« 
The change IS (do—d,,), and the estimated values 
as deduced from the Rossby relations are shown in 
figure 6. The magnitude of change depends on th: 
wind velocity and the sea roughness. If the sea 
roughness be chosen as ¢,=0.3 em, the value of th 
average change is about (@.—,,) =7° for the higher 
wind velocities. 


3.3. Effective Lake Wind Velocity 


The third step is to determine the effective lake 
wind velocity, which is the velocity of the wind 
necessary to produce the observed wind tide under 
the condition that the wind blows with uniform in- 
tensity along the axis of the lake. When the actual 
winds at Toledo, Cleveland, Erie, and Buffalo during 
a storm are examined, it will be seen that these winds 
differ from each other in intensity and also in diree- 
tion. The need for a formula to give the effective 
wind is readily seen. 

The outlines of Lake Erie are shown in figure 7 
The direction of a wind blowing from the south is 
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Sketch map of Lake Erie showing locations of cities 


considered and notation for directions 


hosen as the reference direction from which = to 
measure wind directions. The lake axis is drawn so 
that it divides the area of the lake into two equal 
parts. Its direction is approximately a=68 In 
the cross section through the point z on the lake axis 
the velocity of the wind at an elevation of 25 ft 
ibove the water surface is u,,, and its direction is 8@, 
an angle less than 180°, Denoting the effective lake 


wind by Vo and the length of the lake axis by Z, 
, 1 (°/ 
V? I ui, cos (@—addr. (9) 


This formula is applicable when (@—a) is small and 
the resulting wind velocity is large, for then the effect 
of the earth’s rotation may be neglected and it may 


be imagined that the resultant movement of the 
water is in the direction of the lake axis. Putting 
V/ = Uy) / Ua Where %))/M%. Is determined by eq 5, 
eg 4 becomes 
1 f2 
V2 AMP?u2, cos (@—a)dr. (10) 
L Je ae 


In this transformation no provision is made for 
the circumstance that the wind directions are not 
the same over the lake and over a city on the shore 
of the lake. Theoretically, in the above, denoting 
the difference (@—a) as 8, the factor cos 8 must be 
replaced by cos (8+-y) where y is the difference 
Os2— ds) Now (B+y) 8B cos y—sin B 
sin y, Where y is a small quantity not exceeding 7°. 
\lso the larger and significant values of wind tide 
are associated with winds for which the absolute 
values of 8 are less than r/4. Therefore it is readily 
seen that the resulting error is small when y is 
neglected. 

When the exact variation of u,,—.Wuy» as a function 
sis not known, the formulation must be based on 
discrete values of wind velocities. Let Vy, Ve, Ve, 
and Vy be the wind velocities ua» at Toledo, Cleveland, 
Erie, and Buffalo. The corresponding values of M 
and @ are denoted by the same subscripts. Let the 
lake axis be divided into three segments of equal 
ength. It may be supposed that the winds pre- 
vailing over the end segments are deducible from V, 
aml Vy. The middle segment may be divided into 
vo equal parts in which the winds are deducible 


cos cos 


of 
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| Figure 8 The reduced data for wind tides, AH, and eflective 
winds, V. for storm of December 31, 1911 
from Ve and Vy. Thus instead of eq 10, we may 
write 
V?=4[.M7V 7, cos(@p—a) + M5V3 cos(6,— a) 
MM? V2 cos(@ a) 4M? v; COS(A» a)| 
Now for the higher wind velocities, in which we 


are interested, the variation of 7 is not large as 
may be seen from the curves in figure 5. Therefore 
without significant loss of accuracy, we may write 
| \"2 

172 


MV, 


LIV} cos (@7— a) + Vi, cos (0g— a) 


SV cos (@¢—a)+ 4V2 cos (6¢— a). (11) 
In this representation \, is, in a sense, the effective 
city wind velocity for the southern coast of the lake 
over the entire length between the lake extremities. 
Values of .V7 as a function of V, are shown in table 
1. The table is established on the basis of the curve 
in figure 5 associated with ¢,=0.3 em, after identi- 


fying vu» as V, and ww, as V. 


Tapie 1 Dependence of the multiplying factor M on the 
effective city wi nd velocity V 
1 Ml Vy M 
(mph mph) * 
loo 11% 1 300 1. 254 
an 1.175 1, 400 1. 257 
300 1. 198 1, wm 1a 
00 1. 205 1, 600 1. 2 
um 1.214 1, 700 1. 24 
wow 1, 223 1, 500 1 6 
700 1. 229 1, a0 1. 268 
sou) 1. 24 2, 000 1. 270 
ooo 1. 239 2, 100 l ? 
Loo 1248 2.2" 1 J 
1, 100 1. 247 2, 300 ! ‘ 
1, 200 1. 251 2, 400 1. 275 















The method of velocity reduction explained above 
was applied to all the wind data of the storms under 


consideration. The final values of effective wind for 
the storm of December 31, 1911, are shown in figure 8 
as a matter of illustration. 


4. Determination of Wind Tides 


If the nodal point of the displacements of the sur- 
face water in Lake Erie during a storm be known and 
also the level of the undisturbed water prior to the 
storm, the coefficient of wind stress may be obtained 
by considering the wind tide at either end of the lake. 
Practically, however, this method of evaluation has 
to be avoided as it is difficult to determine the nodal 
point and the level of the lake with exactness. It is 
more satisfactory to base the determination of the 
coefficient on the relative displacement of the water 
at the lake ends. This relative displacement is called 
the total wind tide, and is denoted by A//. The 
displacements are observed at Toledo and at Buffalo. 
In the case of a westerly gale if hy is the wind tide at 
Buffalo and A, is the wind tide at Toledo, the total 
wind tide is AJ =hy—hr. 

This method of total wind tides has another advan- 
tage in the case of Lake Erie which may be mentioned 
very briefly. The particular phases of the storm 
effects that are specifically useful in determining the 
coefficient of wind stress are the maxima of the wind 
tides. Now, because the lake ends differ in shape, 
the rates of rise and fall of the water levels are differ- 
ent for the two ends. For Buffalo, the top part of the 
curve of tides is generally pointed; for Toledo, flat. 
Moreover, the maxima at Buffalo tend to 
earlier than at Toledo. These differences become 
more pronounced when the maturing storm increases 
rapidly in intensity. By adopting the total wind 


occur | 


tide as the datum to be treated, the effects of these | 


peculiarities of the individual behaviors are somewhat 
reduced and thus the data become more reliable. 

In some of the charts of the westerly gales received 
from the Lake Survey the displacements of water for 
the western end of the lake are shown for Gibraltar 
or Amherstburg. In these cases the wind tide indi- 
cations shown had to be reduced to the Toledo values. 
For the reductions it was necessary first to establish 
the correlation between the displacements at these 
western localities. This was possible, since in some 


of the charts the simultaneous values of the displace- 


ments are shown. Also Blunt [2] gives values of the 
maximum displacements for these localities. Intro- 
ducing all these in figure 9 and drawing a straight 
line through the plotted points, it is seen that on the 
average 

hr=1.71h, (12) 


where Ay is the displacement of water either at 
Gibraltar or Amherstburg. This is the correlation 
to be applied. 

A final matter to be considered relative to the 
observed values of the wind tides is the correction for 
the effect of barometric pressures. Ordinarily, when 
the gale is moving over the lake area the pressures at 
the two ends of the lake are unequal. When the gale 


















































ee ee : aor 
| 
+ | ~ -_——— o-— 
| “ LL 
4 
| | , 
-\¢ 4 
7 ° 
& > v 
- 5 > 7 
< Po oe v 
: os 
a © 
Vv 
' oe 
° 
° ' 2 3 4 5 6 7 . 
hr, feet 
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is a westerly, the greater pressure is at Buffalo. The 
difference of the atmospheric pressure at Toledo and 
Buffalo every hour of the storm period may be 
expressed as a head of water to be subtracted from 
the observed total wind tide. 


5. Data for Winds and Wind Tides of Lake 
Erie Storms 


The curves of wind intensity and total wind tide 
for the storm of December 31, 1911, shown in figure 8, 
are typical of all the storms considered in the present 
study. Taking the presentation in figure 8 as an 
example, it may be seen that the entire manifestation 
of the storm may be broken into three epochs; the 
epoch of maturing storm, the epoch of relative 
steadiness, and the epoch of recession. The hydro- 
dynamical behavior of the water in the lake is 
expected to show marked differences during these 
epochs. In the initial epoch the wind must blow for 
some length of time before the response of the water 
to the action of the wind is omental. The reason 
for this condition is that any manifestation of wind 
tide must be associated with the flow of water from 
one end of the lake to the other. The action of the 
wind must establish a layer of drift current, the 
thickness of which must increase with time either 
under action of molecular viscosity, or turbulent 
viscosity, or both. As the water is being collected 
at the leeward end of the lake, a returning gravity 


, current is created which likewise increases in intensity 


with time. The second epoch represents that steady 
condition in which the transport of water through 
the body of the drift current is counterbalanced by 
the transport of water in the returning gravity current 
that is maintained by the unchanging surface 
gradient of the lake waters. The third epoch is 
associated with decreasing wind intensity. ith the 
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strength of the wind decreasing, the stress of the 
vind cannot maintain the adverse gradient of the 
surface water with the result that a surge of water in 
wave motion with a period equal to the seiche period 
of the lake is produced. In some records these 
seiches with periods of 4 hr are clearly seen during 
he time of the falling storm. The inertia effects are 
vy pronounced in the extreme epochs: in the 
termediate epoch they are considerably reduced. 

For the determination of the relation between 

ind velocity and wind stress, only the manifestation 
if water level changes in the intermediate epoch will 
« considered. Referring once again to figure 8, 
he reference time is the instant of maximum wind 
tide In each case of the storms considered the 
eference time is evaluated by taking the average of 
the times of maximum deflection of water at the two 
nds of the lake. As the total wind tide of the 
ntermediate epoch we have taken the average value 
of the wind tide over a duration of 4 hr around the 
reference time. For the corresponding wind values 
we have taken the average value over a duration of 
> hr immediately preceding the reference time. 

The wind velocities and wind tides thus obtained 
from each record are entered in table 2 with the dates 
of the storms indicated. In some of the wind tide 
charts the storms showed two peaks with an extended 
flat region between them. The total tides of such 
regions are also included in the data of table 2. and 
are identified with asterisks. 


Panie 2 Wind tides of Lake Erie during westerly gales 


No Date 1 1 sll 
mph mph ft 
l Nov. 21, 1900 WS 249 13.12 
2 Oct 0, 1405 1.4 wt 6. OS 
; Oct 0, 1906 is. 3 i471 ¥.7 
‘ Jan 0, 1407 is. 1 $12 12.4 
fi Dex 7. 1909 “1 10s Ww. 51 
"i Dec. 31, 1911 —u 1516 9.53 
7 Jan. 31,1014 4.7 1207 7. 05 
s Dex 4, 1917 45.2 Ist5 Ww. 17 
hd Dee @, 1017 3.1 lous ee 
10 Dec, WO, 1907 4.3 1176 7. #2 
il Dec. 18, 1921 5.1 wvss 12.40 
12 Dex 8, 1027 7.4 243 13. 24 
13 Dex ¥ 1927 27 7" 41 
14° Dex 4, 1027 26.1 in 44 
15° Des ¥, 1927 22. 1 487 1.73 
16 April 1, 1929 51.3 Lid 13.31 
17 Jan. 22, 1939 ws 1508 9” 
Is Sept. 25, 141 5 1277 oo 
iy Jan im? 0.4 1632 12.53 
al) Jan s, 142 19.5 79 2. 38 
1 Nov. 22, 1646 4.2 171 &. 36 
2 Mar. <5, 1947 5.7 1276 8.4 


The square of the wind velocities, V?, near the 
lake surface and the corresponding wind tides, A//, 
taken from table 2 are plotted in figure 10. Com- 
paring the position of the plotted points with respect 
to the full line drawn, it is seen that the data cor- 
responding to high wind velocities show small 
deviations from the curve, whereas those corre- 
sponding to low wind velocities show proportionately 
larger deviations, and all towards lower values of 
\/1. One inference is that the wind stress for low 
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wind velocity is relatively smaller. Unfortunately, 
the observations for the lower wind velocities are not 
sufficiently numerous to establish this point bevond 
doubt. 


6. Integration of the Wind Tide Equation for 
Lake Erie 


First the wind tide differential equation for a lake 
of arbitrary shape will be derived. Let positive z 
be drawn along the direction of the wind, this diree- 
tion coinciding with the axis of the lake. At the 
point z the wind tide is A, it being the elevation of 
the raised water surface with respect to the level of 
the undisturbed water. See figure 11. Corres- 
ponding to A at z the surface width of the water is 
B, the cross-sectional area is A, and the mean depth 
is //. The cross section is normal to the axis of the 
lake. The velocity of the wind is V, and the cor- 
responding wind stress is r,. At the bottom the 
stress induced by the gravity current is rt». The 
density of water is p and that of air is p,. Two see- 
tions Ar apart, one at s and the other at z+ Az, are 
taken. If the total pressure force on the section at 
s be denoted by P, the pressure force on the section 
at z+ Az will be P+ Apg (dh/dz)Ar. Hence, the dif- 
ference of the two pressure forces is Apg (dh/dz) Ar. 
The inertial action of flow in the drift current above 
and in the gravity current below being ignored, the 
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Ficgure 11, Notation diagram for wind tides 


17 is the mean depth in the cross section of the lake at x 
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only remaining counterbalancing forces are the frie- 
tional forces amount to B(r,+-7))Az. Thus, equat- 
ing the forees and dividing by Ar, Apg (dh/dz)=B 


(r,+-r,), and since A=B (77+h), then 


dh T, TT 


= ’ (13) 

dr pg IT +h) 
which is the wind tide equation for a lake of arbi- 
trary shape. As the volume of water in the lake is 
not changed, the condition of continuity reduces to 


I 
Bhd, 0 (14) 


where L is the length of the lake axis 

Equations 13 and 14 may be changed into dimen- 
sionless forms. The average section depth may be 
expressed in terms of the average lake depth /7), 


l (on 
a | Bids (15) 
BL Jo 


The section surface width may be expressed in terms 


of the average lake width, B, 


os] 


(16) 


4 (" Bds 


The distances along the axis may be expressed in 
terms of the lake length, L. 


» 
Putting 


Te+ To L a7 
pall, Hy’ “ 
R (18 
eae: ” 
and 
: 
[=> (19) 
eq. 13 and 14 reduce to 
dy /1, HM = (20) 
de (14-91/) 
and 
"~B 
= ndt—0 (21 
{, | dete 


It is important to remark that ¢ is the parameter of | 


stress, », the parameter of wind tide, ¢, the parameter 
of distance. 
B/B are functions of ¢. 


As ordinarily 9 is a small quantity, eq 20 may be | 


written in the form 


dy H, H, 
=f l—— . 22 
dé a Al 7”) (22) 





In a lake of arbitrary shape //,/H7 and 


| 
| 


The solution of the equation may be expressed in t| 
form 


n=oF\(t)—o°F x(t), (2 
where 
Fy (e) C+ | Mo je (24 
Jo Hf 
. ’ 6 I, . - - » 
fF (¢) Ct]. (a) Rods, (2 


the quantities C; and C, being the constants of iy 
tegration. Now, from eq 21 


a 
| B ie ceydi 0, 26 
0B 

and 
| B k(o) dé 0, 27 
Jo B 


and these serve to evaluate the constants C, and ¢ 
To carry out the determinations of the functions 
F, and F,, it is necessary to obtain the shape quanti 
ties /7/77, and B/B for Lake Erie. These quantities 
as determined from the Lake Survey Office map of 
Lake Erie are shown in figures 12 and 13. Using the 


data from the figures, and following the method of 


computation shown above, the functions F; and Ff 
are determined. These are plotted in figure 14. For 
the subsequent analysis the end values 



































Fy(0) 0.641, F.(0)=0.323, 
(28 
F,(1) 0.512, F,(1)=0.118, 
are important. 
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Mean depths of Lake Erie as function of distance 
from Toledo. 
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riaure 13 Surface widths of Lake Eric as function of distance 
: s 
> 0.6 T T T T - 
1.25 
i x/L 
p Figure 15 Comparison of computed wind tides for Lake Erie 
f : with observations 
: Oct. 29, 1892 (Blunt Oct. 14, 1843 (Blunt). J, Sept. 25, 141. ©@, 
Jan. 2. 1942 Jan. 3. 1942 
! 
The quantity on the right hand, computed with 
reference to the values in figure 14, is shown in figure 
15 by the full curve The plotted points are from 
observations of different storms. The points repre- 
sented by open circles and open squares are taken 
from data given by Blunt [2]. Apparently there 
appears to be a fairly good agreement between 
theoretical and observed values 
7. Evaluation of the Coefficient of Wind 
Stress 
Representing the end values of 7 as no and n, it is 
= seen from eq 21 that mn =oF)(0)—oe*F2(0), and n, 
t ‘ ‘ oF (1) —o* F2(1). Hence (m)— An=o|F\(1)— 
: as ; ; F\(0)|—oe|F,.1) — F,(0)], and introducing the values 
IGURE 14 alues of the functions F, and F. for Lake Erve in eq 28, An=1.153 o +-0.205 0%, or o=0.867 An 
0.178 07. Treating the second term in the right hand 
member as a small quantity, 
At this point, some observed values of the lake- 
water displacement may be compared with the theo- o =0.86749—0.134An° 
retically determined curve of the analysis above. 
lhe most suitable way to represent the curve is to | Introducing values of ¢ and » from eq 17 and eq 18, 
consider the ratio 2h/(h;—ho), where h is the surface | and remembering that A//=h,—/o, we have 
displacement at z, and hy and h, are the extreme 
values, as a function of s/L. Since o is a small tat+to L 0.267 All] 0.134 (= . (20 
' quantity, F, may be neglected and thus, poll, H, aoe ae . HI, ): oe 
9 oF (+ rar : 
2h 2F'(¢) (29) rhis is the relation that relates the wind stress to the 


hi—h, F\A)—F (0) total wind tide for Lake Erie 
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Expressing the shear at the lake bottom as a 
fraction of the wind stress, 


Tr NTs, (31) 
the relation in eq 30 may be written as 


ait ) HM, pglly. (32) 


l = ae ) 
T; es [ 0.86: HW, 0.134 ( Hi, L 


Introducing the definition of the coefficient of wind 
stress, 


o xp. ™, (33) 


where p, is the density of air, and V the velocity of 


the wind, we have from eq 32 
0.867 MT) AIT p gil, 
we [ 0.16 H,| Hop. V2’ (34) 


which is the proper formula to evaluate the coefficient 
of wind stress from the total wind tide observed at 
Lake Erie. In the application of the formula the 
depth of water //, is taken as 58 [t., the ratio of the 
density of water to the density of air, p/p, as 1,250. 
The wind velocities V are to be expressed in feet per 
second and the wind tides A// in feet. It may be 
supposed that the frictional stress of the lake bottom 
is about one-tenth of the wind stress. Accordingly, 
nis 0.10. 

Applying the above formula to the set of wind 
intensities and corresponding wind tides appearing 
in table 1, a set of values of the wind stress coefficient 
are obtained which are entered in table 3. The same 


TasBLe 3 (‘ve ficre nt of wind astreas deduced from Lake Erie 
gales 
No ! xl No ! xxl 
mph m ph 
1 weS 214 12 7.4 24 
2 $1.4 287 13 27.2 2.35 
} a I 2.78 4 1 215 
‘ “1 217 15 22.1 1. 5a 
) “1 273 16 51 2. 8s 
“ ay 265 17 iS 2. 4 
7 “4.7 2 76 Is 6. 7 2m 
s 43.2 2.2 iv “4 4. 22 
” Mt 1.78 2 lv 2 
w M4 2.71 21 “4.2 2% 
* > 5 > > 


data are plotted in figure 16 as full circles. Plotted 
in the same figure as open circles are the wind stress 
coefficients which Neumann obtained in an exami- 
nation of the Pelman observations for the Gulf of 
Bothnia [7]. Neumann's original values were reduced 
to correspond to winds 25 ft above the surface of the 
water. 

An examination of figure 16 reveals that in the 
region of the higher wind velocities the coefficient of 
wind stress is practically independent of the velocity 
of the wind. Furthermore, the coefficient of wind stress 
for the two bodies of water are of like value. The 
average value of the individual determinations is 
0.00254 for Lake Erie, and 0.00236 for the Gulf of 
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Comparison of coefficients of wind stress in Lake 
Erie and the Gulf of Bothnia. 


@. Lake Erie 


Fieure 16, 
. Gulf of Bothnia 


Bothnia. Now, the data from the Gulf of Bothnia 
for low wind velocities indicate that for this region 
the coefficient of wind stress is not constant but 
decreases with increasing wind. At the moment no 
similar data are at our disposal for Lake Erie and 
efforts are being made to remedy this gap in our 
investigation. 


8. Evaluation of Sea Roughness 


We shall adopt for the determination of sea 
roughness the procedure used by Neumann. Writing 
eq 33 in the form (V/V,)?= x7', substituting in 
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Fieure 17. Variation of sea roughness, €;, with wind ve- 


locity, V. 
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eq 1, and changing ¢ to ¢,, there is obtained 


l " 


5.75 log -—. (335) 
€, 
‘Xx . 


which relates the sea roughness €, to the elevation 

for which the wind is determined and the co- 
efficient of wind stress x. In the present case, since 
y Is determined from wind velocities prevailing at 
a height of 25 ft, 2,,—25 ft, or 762 em. As mentioned 
above, the average value of the coefficient of wind 
stress from the Lake Erie data is 0.00254. This 
vields for the sea roughness the value €,;=0.27 em 
As will be remembered, in the computations required 
to reduce the city wind velocities to the lake velocities 
it was necessary to assume the value of the sea 
roughness beforehand. The assumed value was 
« —0.3 em, and this selection is now seen to be 
satisfactory 

The full line drawn through the data points of 
figure 16 may now be used to investigate the variation 
of sea roughness with wind velocity. Using values 
from the curve and the formula above, eq 35, the 
results of the computations are given in figure 17 
in the form of a curve. With increasing wind 
velocities the magnitude of the roughness decreases 
When the higher wind velocities are attained, the 








roughness approaches a constant value, practically. 

We refrain at the present moment from discussing 
the theoretical significance of the facts that both the 
coefficient of wind stress and the sea roughness 
decrease with increasing wind velocities 
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Twisted Square Plate Method and Other Methods 


for Determining the Shear Stress-Strain 
Relation of Flat Sheet *' 


Walter Ramberg and James A. Miller 


A method is presented for determining the stress-strain relation in shear for isotropic 


flat sheet 


stress-strain curve 


The method requires measurements of deflection or of extreme fiber bending strain 
in the center portion of a twisted square plate 


The 


octahedral stress derived from the 


in shear for a twisted plate of aluminum alloy agreed within about 5 


percent with the octahedral stress derived from tensile and compressive tests of the same 


material 
point. Unfortunately, 
size of the specimen 


the use of an annulus of constant thickness 


l. Introduction 


There has long been an interest in the stress-strain 
relations in shear of materials in the form of sheet. 
The interest has become more pronounced as working 
stresses and strains have increased and as it has 
become apparent that the mechanical properties of 
thin sheet could not be deduced directly from 
mechanical tests of heavier sections. The most satis- 
factory method, perhaps the only really satisfactory 
one, for determining the stress-strain curve of thin 
material in shear requires the material to be in the 
form of tubing with circular section. As long as the 
tube mat rial is isotropic, its stress-strain curve in 
shear can be computed conveniently from the torque- 
twist curve of the tube by following the procedure 
given in [{1].? 

There is no dearth of suggestions for methods of 
determining the stress-strain curve in shear of flat 
sheet. Figure 1 indicates 10 suggested methods for 
stressing thin sheet in shear that have come to the 
authors’ attention. The first two methods deal with 
tubular specimens. 

Method (a) is nothing more than the bending of a 
strip of the sheet into a tube of large diameter, which 
is then tested like the conventional long tube for 


which a satisfactory theory exists [1]. The tube 
would have to have a large diameter so that the 


strains introduced in bending the strip would be 
small compared to the strains to be measured. 
Difficulties might be experienced with deviations 
from the theoretical uniform shearing stress caused 
by variations in clamping forces along the edges of 
the specimen. Difficulties might be experienced also 
with the discontinuity at the joint and with the 
transmission of sufficient torque to obtain stresses 
in the plastic range. 

Method (b) was suggested to one of the authors 


his paper is based on work sponsored by the National Advisory Committee 
Aeronauties and the Mechanics Branch, Office of Naval Research, U. 8 
partment of the Navy 
Che major portion of this paper was presented before the First | 
ress of Applied Mechanies in Chicago on June 11 to 16, 
shed in the Proceedings of that Congress under the tithe 
ress- Strain Curve in Shear by Twisting Square Plate 
gures in brackets indicate the literature references at the end of this paper. 
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The difference was much greater for specimens of mild steel with a definite vield 
the method is difficult to perform on thin sheet because of the small 
An examination of alternate methods indicates particular promise for 


by J. M. Frankland of Chance Vought Aircraft 
Four strips of the sheets are assembled into a tube 
of square section. The strips will be subjected to a 
uniform shearing stress by equal and opposite torques 
applied to the ends, if the torsional rigidity of the 
joints along the four corners is negligible compared 
with the torsional rigidity of the tube as a whole. 
It is probable that this corner effect is not negligible 
for tubes that are sufficiently thick-walled to trans- 
mit, without buckling, torques that extend the shear- 
ing stresses into the plastic range. 

Several methods have been suggested for stressing 
a specimen of sheet in shear by forces acting in the 
plane of the sheet. 

Method (ec) is the well-known square-picture-frame 
method, in which a square specimen is held between 
four pairs of rigid clamping bars with pin joints at 
the corners. The specimen will be placed in a state 
of uniform shear by equal and opposite forces acting 
at opposite corners of the frame, provided that the 
sheet is perfectly flat, the clamping bars are infinitely 
rigid, the edge effect due to deviations from ideal 
clamping does not extend into the middle portion 
of the sheet, and the pin joints are frictionless. 

Method (d) tries to get around the difficulties with 
the edge effect in method (c) by making the speci- 
men long so that the shearing stress is practically 
uniform along the middle portion of the specimen. 
However, an end effect will remain because the shear- 
ing stress reduces to zero at the corners. This end 
effect may be evaluated by making tests of the same 
material in picture frames of different lengths. The 
procedure would be costly and time consuming. 

Method (e) was used by Coker and Filon [2] and 
Bollenrath [3] to stress thin sheet in shear. The 
loads are applied symmetrically to subject each one 
of the two “bays” of the specimen to a state of stress 
in which the shearing stress along the center line of 
each bay varies from zero at the ends to a nearly 
constant value near the middle. As in method (d), 
the end effect may be eliminated by testing speci- 
mens of different lengths. 

Method (f), suggested by W. R. Osgood, subjects 
a single strip to a state of stress which approximates 
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Fieure | Methods suggested for stressing thin sheet in shear. 


that in each bay of the specimen for method (e). 

Method (g) avoids end effects by testing the sheet 
in the form of an “endless strip” or annulus. Equal 
and opposite torques are applied at the inner and 
outer edges of the annulus. The method has the 
advantage that, except near the periphery and the 
inner circle, the shearing stresses would be known 
accurately, for isotropic sheet, in the plastic as well 
as in the elastic range. There are some disadvan- 
tages, however. The shearing stress varies inversely 
as the square of the radius of the circle along which 
it acts. Hence there will be a nonlinear strain varia- 
tion along the gage length of any finite gage mounted 
at 45° to the tangential direction. As a result, the 
average strain along the gage will differ from the 
strain at the midpoint of the gage, unless the annulus 
is large enough to make the nonlinear variations 
negligible. If the annulus is large enough, it may be 
difficult to apply sufficient torque to stress the ma- 
terial into the plastic range. The testing machine 
would require heads that are rotatable boat their 
axis only, to prevent bending of the specimen by 
rotation about other axes. The method would give 
only an average stress-strain curve for anisotropic 
sheet. 


Method (h) eliminates some of the difficulties o 


method (g). The annulus is machined into a dis! 
that varies in thickness inversely as the square o 
the radius. The material in the annulus will be sul 

jected in that case to a uniform shearing stress | 

equal and opposite torques applied at the inner an 
outer edge. The disk must be sufficiently thick nea 
the periphery to develop strains in the plastic rang 
before buckling occurs and the variation in thicknes 
with radius must be small enough at the radius whe: 
the strains are measured to assure a practically un 

form shearing stress throughout the thickness. A 
practical difficulty is that of machining the specime: 
accurately to the desired thickness variation with 
out changing the properties of the material by th 
machining operation. 

The methods so far mentianed subject the sheet 
to a shearing stress that is substantially uniforn 
throughout the thickness of the sheet. In the r 
maining two methods the shearing stress is varied 
in the thickness direction from one extreme to an 
opposite extreme. 

In method (i) a state of shear is produced in a 
square plate by twisting the plate into an anticlastic 
surface by equal and opposite tranverse forces ap- 
plied to adjacent corners. The state of shear wil! 
be uniform for planes at a fixed distance from th: 
median plane which will undergo zero shearing stress 
The upper and lower faces of the plate will be under 
equal and opposite shearing stress. Objections might 
be raised against this method in so far as it neglects 
deviations from uniform shear due to various causes, 
such as having large transverse stresses near the 
corners and no shearing stresses along the edges 
These objections will be discussed in detail in the 
body of the paper. Experimental difficulties may 
arise in applying the corner loads to the specimen 
without interfering with strain and deflection meas- 
urements on the faces of the plate. 

Method (j), suggested by A. H. Stang, uses a 
much simpler scheme for producing a state of stress 
similar to that in method (i). A specimen in the 
shape of a long strip is placed in a state of shear by 
applving equal and opposite torques to the ends of 
the strip. Shearing strains can be measured with 
strain gages or they can be obtained from the twist 
per unit length of the strip. However, if the strip is 
made long enough to eliminate end effects, it is 
likely to deflect laterally. This will superimpose 
bending stresses on the torsional stresses. The same 
objections might be raised in this case as in case (i 
with regard to edge effects. The shearing stress at 
the extreme fiber will be practically constant in the 
center portion, far from the edges of the strip, but 
it will decrease rapidly to zero as the edges are 
approached. The value of the extreme fiber shearing 
stress near the center of the strip will depend on 
this edge effect. It can be computed exactly in the 
elastic range, but no exact solution is known for the 
plastic range. 

Consideration of the above 10 methods suggests 
that at least three of them, the square-picture-frame 
method (c), the annulus of constant thickness method 
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¢), and the twisted-square-plate method (i), should 
be tried out experimentally. 

An experimental study of method (¢), the square 
picture frame, was conducted by Wm. R. Osgood 
about 5 vears ago. Method (c¢) was selected because 
it appeared to have particular promise and because 
it had been used in the aircraft industry to estimate 
shear properties of sheet A fixture had been de- 
veloped at the El Segundo Plant of the Douglas 
Aircraft Co, to make such tests. A duplicate of this 
fixture was obtained, and tests were made at the 
National Bureau of Standards for the National 
Advisory Committee for Aeronautics on several 
15™4.50.063 in. 245-T aluminum alloy plates 
The results were disappointing. They were reported 
to the NACA, but were not published. The strains 
were fear from uniform even in the center portion of 
the plate farthest from the edges. The extreme 
fiber strains in this optimum region deviated up to 
10 percent from the average strain. The deviations 
persisted even after an attempt had been made to 
improve the conditions of clamping along the edges 
by remachining parts of the clamping fixture 
“Two vears after completion of the work for the 
NACA, it was decided to try method (i), the twisted 
square plate, since Nadai had shown in 1915 [4] 
that this method provided a simple means for 
determining the shear modulus of sheet material 
More recently the test has been extended to ortho- 
tropic materials at the Forest Products Laboratory 
15] to determine the shear moduli of wood 


2. Nomenclature 


The following nomenclature is used in the paper 
a=half-span for deflection measurements 
b—width of specimen. 
k—modulus of elasticity mn 

compression. 
G—modulus of elasticity in shear 
h=thickness of specimen. 
m=moment per unit length. 
M= bending moment. 

M,=twisting moment. 

P load at each corner of plate. 

?—redius of curvature along diagonals of 
plate 

displacements parallel to r-, y- and -- 
axes, respectively. 

r, y= axes parallel to diagonals of plate, also 
planes normal to the corresponding 
axes. 

axes parallel to edges of plate, also 
planes normal to the corresponding 
axes, 

axis parallel to thickness of plate, also 
plane normal to 7-axis. 

y=shearing strain. 

y,— octahedral strain. 

¢=normal strain. 

v= Poisson's ratio. 

normal stress. 
7=shearing stress. 

r,— octahedral stress. 


tension or 


113 





3. Principle 
3.1. Elastic case 


An elasiic square plate will be in a state of pure 
shear, in which the shearing stress is proportional to 
the distance from the neutral surface of the plate 
when the plate is bent into an anticlastic surface 
with equal and opposite principal curvatures. The 
bending of the plate into such an anticlastic surface 
produces equal and opposite extreme fiberbending 
stresses along the two principal axes. Nadai has 
shown that this remarkably simple state of stress 
may be realized by applving equal and opposite 
tranverse forces to adjacent corners of the plate. 

Nadai’s proof |4] makes use of Saint-Venant’s 
principle of the equivalence of equipollent loads 
Nadai notes that for a load P at each corner, each 
edge of the plate is loaded by equal and opposite 
forces P/2 acting at two corners (fig. 2). The load 
along the edge will remain unchanged upon applying 
balancing forees + P/2 at frequent intervals along 
the edge, as indicated in fig. 2. Any element of 
length of edge Ax’ between two adjacent equal and 
opposite forces P/2 will then be subjected to a twist- 
ing moment (/?/2)Az’. According to Saint-Venant’s 
principle, this edge load is equivalent to any edge 
load which results in a constant twisting moment 
per unit length 


b 
te 
te 


uniformly distributed along the edge of the plate 
The stress distribution, in regions of the plate suffi- 
ciently far from the edge to satisfy Saint-Venant’s 
principle, will then be that corresponding to uniform 
twisting moments + m,.,.,—Mm, in the z’, y’ planes, 
normal to the #’, y’ axes, respectively. The twisting 
moments + m,), are associated with principal bending 
moments m,, m, on the y, 2 planes normal to the 
diagonals of the plate. Equilibrium of moments 
acting on the edges of a triangular element of the 
plate (fig. 3) bounded by planes parallel to the z, y, 
and y’ planes shows that 


m, my, Miytyt P/2. (2) 
The plate is bent into an anticlastic surface by equal 
and opposite bending moments per unit length //2 


acting on the y, # planes normal to the diagonals of 
the plate. The surface will have its planes of prin- 
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Fieure 2 Twisting of square plate by forces on corners 
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Ficure 3. Edge moments acting on element of twisted square 
plate. 


cipal curvature normal to the diagonals of the plate. 
If R is the radius of curvature along the diagonals, 
the deflection relative to a plane tangent to the center 
of the plate is given by 
r y* 
w “e (3) 
2R 


The distribution of strain and of stress throughout 
the plate may be computed on the assumption of the 
Kirchhoff plate theory that plane sections remain 


plane. 
linearly with distance from the neutral fiber 


22 22 
ég™ h €; €, ho €2) (4) 


where €,, €: are the principal extreme fiber strains; 
z, the distance from neutral fiber; and A, the thick- 
ness of plate. The maximum shearing strain on the 
2 plane is 
22 

Y Yrty’ €z— €, (5) 
where y; is the shearing strain at the extreme fiber 
yi=€:—e. Because of the antisymmetry of the 
deformation €,;=—e, so that 


71 2é). (6) 


The shearing stress corresponding to the strain ¥ is 
given by 

T Gy (7) 
nside the elastic range. The twisting moment per 
unit length is given by 


A/2 
Metyt f _r2dz. (8) 


—h/2 


Substituting (1), (7), (5), and (6) in this equation 
gives 

rd. ie 
== €,Gh’. (9) 


3 3 





The normal strains ¢,, ¢, must then increase | 





From which 


G (10) 


2h? «, 


Hence the shear modulus can be computed from the 
slope P/e, of a plot of load against extreme fiber strain 
by substituting P/e, in (10)° 

The shear modulus may be determined also from 
the deflection w, of the center of the plate relative to 
a line connecting points a distance a from the center 
along one of the diagonals, figure 2. If plane sections 
are assumed to remain plane, the curvature 1/R along 
the diagonal z is given by the familiar formula 


l Pw €; 
R- dz? h/2’ vad) 


where w=normal deflection. Substituting (11) in 
(9) gives 
ee 36°P 


2 9 *») 
h 2h°G (12) 


Wy 


so that 
3a’? P (13 
, 3) 
2h  w, 
which gives the desired relation between shear mod- 
ulus and center deflection under a given load. 


3.2. Extension into plastic range 


Equations (1) to (6) apply beyond the elastic range 
as long as plane sections remain plane. Equation 
(8) for the relation between twisting moment and 
shearing stress applies beyond the elastic range also. 
We may replace z in the integral of that equation by 
y, using eq (5) so that eq (8) becomes 


P h? ‘ 


1 
Mm, => | rydy. (14) 
“Yi Jo 


zy’ » 
The change in load dP corresponding to a small 
change in extreme fiber shearing strain ¢@y, is obtained 
by differentiating this equation in respect to y, with 
the result 


, » 2 lek 
dt 2h | redeem T1171 t (—2P+hAh?r,), 
Yi 71 


dy, 7 vi Jo 
(15) 
where 7, is the extreme fiber shearing stress. Solv- 
ing this equation for 7; gives 
1 dP : 
Ti rd ke i +2P P (16) 
: ay. 


The extreme fiber shearing strain can be mente 
simply from measured values of extreme fiber bend- 
ing strain « or from measured values of deflection 
w, by making use of eq (6) and (12) 


_ 2hw, 


a (17) 


Y= 2€ 
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lhe stress-strain curve in shear can be derived, there- 
fore, from measurements of P(e) or P(w,) by using 
\7) to compute the shearing strain y,, plotting 
P(y,), and computing 7, from the ordinates and 
ope of P(y,) in accordance with eq (16). 

It is interesting to note that eq (16) is of the same 


orm as the equation 

2 dM,, 

bh: [« de, +2} 
which is derived by Nadai [6] to compute the stress- 
strain curve o,(e,) from the bending moment-strain 
urve M(e,) for a rectangular beam of width, 6, and 
height, A, in pure bending; « is the extreme fiber 
bending strain produced by the bending moment, 
\/. The identity in form of eq (16) and (18) might 
be expected because the square plate is bent plasti- 
cally into an anticlastic surface by equal and oppo- 
site bending moments per unit length about axes at 
ight angles to each other 


a; (18) 


4. Design of Specimen 


In designing the specimen, consideration must be 
viven to practical deviations from the theoretical be- 
havior because the theory neglects edge effects, mem- 
brane stresses, and deformation by transverse shear- 
ing forces. These effects must be considered in some 
detail in order to dimension the specimen so that the 
deviations from ideal behavior can be either neglected 
or taken into account. 

The edge effect will be considered first. In the 
theory, the concentrated forces at the corners are 
replaced by twisting moments uniformly distributed 
along each edge, which are applied by horizontal 
shearing forces rather than the vertical forces shown 
in figure 2. It follows that the theoretical solution 
will hold only for regions of the plate at a sufficient 
distance from the edges, to make Saint-Venant’s 
principle applicable. Dimensional reasoning indi- 
cates that the “sufficient distance’ should be pro- 
portional to the thickness, A. The factor of propor- 


tionality may be estimated roughly by examining | 


cases for which solutions have been obtained from 
the theory of elasticity. 

The plane stress solution for a bar of constant 
width A compressed by a concentrated load at the 
center of the free ends is given in [7]. This shows 
that the maximum compressive which be- 
comes infinite directly under the load, differs less 
than 3 percent from the average stress at a distance 
h from the end. Seewald’s plane stress analysis of 
the local stresses near a concentrated transverse load 
acting at the center of a slender beam of rectangular 
section with simple supports at the ends [8] shows 
that these stresses become negligible compared to 
the extreme fiber bending stresses derived from 
simple beam theory at a distance from the load 
vreater than the depth of beam. 

This suggests that the errors resulting from the 
application of Saint-Venant’s principle will be small 
n regions of the plate at a distance from the edge 


stress 
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which is greater than h. This condition is easily 
satisfied for a plate that is thin in relation to its 
lateral dimensions. Unfortunately, it will be shown 
below that a very thin plate is not acceptable because 
it will develop appreciable membrane stresses under 
the loads required to produce yielding in shear. 
The membrane stresses are negligible only if the 
plate has a width noi exceeding a few times the 
thickness. In such a plate, it is probable that the 
deviations from the assumed stress distribution at 
the edges will have an appreciable effect upon the 
magnitude of the stresses in the central portion of 
the specimen. 

A correction for the edge effect in the twisted 
square plate was suggested to the authors by their 
associate, Samuel Levy. He pointed out that the 
twisting of an elastic square plate can be considered 
as a special case of torsion of a rectangular plate 
considered by Kelvin and Tait [9]. The elastic 
stress distribution at crosssections of a long rectan- 
gular plate far from the ends is given in [10]. The 
extreme fiber shearing stress is plotted in figure 4 
for the case 6/h=10 which was chosen, after a pre- 
liminary analysis, as a practical compromise for the 
twisted plate. Kelvin and Tait [11] showed that 
the shearing stresses T.-, parallel to the plane of the 
strip resist only one-half of the external twisting 
moment .f4,=Pb transmiited by each transverse 
section of the strip. The other half is resisted by 
the transverse shearing stresses r,., acting on the 
section. 

A redistribution of stresses will take place as the 
strip is twisted beyond the elastic range. The 
nature of this redistribution becomes clear by con- 
sidering the extreme case of stress in a twisted strip 
of perfectly plastic material when all the material 
has yielded. This case was solved by Nadai [12] 
with the sand heap analogy. The resultant stress 
distributions are shown dotted in figures 4 and 5. 
Integration of the stresses for the perfectly plastic 
state shows that, in the plastic as well as the elastic 
range, \],/2 is resisted by the shearing stresses r,., 
parallel to the surface of the strip and M,/2 by the 
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shearing stresses r,.., normal to the surface of the 
strip. 

The shear stress distribution of figures 4 and 5 
and the condition that one-half of the twisting 
moment .M/, is resisted by 7,.,, and the other by 1, 
show that the stresses across transverse sections of 
the strip can be approximated by application of 
equal and opposite concentrated forces near the 
ends of the section and a uniformly distributed 
twisting moment along the entire section. 

This distribution of equivalent external forces 
coincides with that shown in figure 2 as acting on the 
edge y’=/2 of the square twisted plate provided 
that the corner forces P are moved in from the 
corners by a small distance to make 72 act at the 
same point as the resultant of the transverse shear- 
ing stresses r,-. which are shown in figure 5 for 
points along the z’-axis. It may be concluded that 
the distribution of forces acting on two edges of the 
square plate in figure 2 nearly balances the shearing 
stresses acting across any transverse section of a 
long strip of width, 6, and thickness, 4, which is 
twisted by a torque \f/=Pb applied to the ends. 
The square plate can be regarded, therefore, as a 
portion of such a long twisted strip with a distribu- 
tion of shearing stress across the section y’ =0 which 
approximates that in the transverse sections of the 
strip. Consequently, the extreme fiber shearing- 
stress distributions indicated in figure 4 would also 
apply to ihe square plate. The shearing stress in 
the central portion would be greater than that com- 
puted from eq (16) because of the above-mentioned 
edge effect. In the elastic case, this increase would 
be about 7 percent for 6/h= 10, according to figure 4. 

In general, we find from [10] that the extreme 
fiber shearing stress in the central portion of the 
longer sides of a rectangular section can be expressed 
by 


T 
Ti, * (19) 


where &, is a factor depending on 6/h. For values of 
b/h>5, the relationship reduces to 


Tie l 9 
. 20) 
tr, 1—0.63h/b \ 


In the cope case of a perfectly plastic material 
which has yielded throughout the section, we have 


from the sand-heap analogy [12] 


: l 
7, 1—0.333h/6 


For 6/h=10, we obtain from (20) and (21) 


rT). 1.0677, elastic cause) 


T) 1.0347, plastic case ) 


The actual correction for edge effect bevond the 
elastic range should lie somewhere between the limits 
given by eq (20) and (21). It should be closer to 
(20) as long as the plastic vielding is localized in 
the outer fibers of the plate. 

The effect of membrane stresses will be considered 
next. The stiffening effect of membrane stresses is 
noticeable experimentally during twisting of a thin 
square plate as the transverse displacements w(zr,y 
become comparable to the thickness of the plate 
Figure 6, taken from Figure 11 of [13] shows the 
effect for a 7.27.20.126-in. plate of 24S-T 
aluminum alloy. The two straight lines in the 
figure are the strains ¢, along the diagonals computed 
from eq (9) for the linear elastic case. It is seen that 
the measured load P agrees closely with the caleu- 
lated values at very small strains, but that it deviates 
increasingly with increasing strains. The deviation 
is close to 30 percent for a strain ¢;=0.0008. 

The effect is analyzed in the appendix for the 


elastic range on the assumption that the transverse 


deflection of the plate continues to be given by 
eq (3). Equation (67) in the appendix shows that 
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Fiaure 6. Extreme fiber bending strains at center and along 
diagonals of 7.27.2 0.126 in. plate of aluminum alloy 
248-T twisted by transverse forces P applied at the corners 


116 

















i relative error in extreme fiber bending strain «, 
used by the stiffening effect, is equal to 


Ae, b\' 


: 0.03666 ( h ) et. (23) 


fhe error goes up with the square of the extreme 
ber strain e; and with the fourth power of the width 
hickness ratio, b/h. It may be made small by 
naking 6/h small. 

Lastly, consideration will be given to the transverse 
hearing stress. The effect of transverse shearing 
stress will be negligible as long as this stress is very 
small compared to the extreme fiber shearing stress 
wrallel to the plane of the plate. The average 
ransverse shearing stress 7, on a section of a cylinder 
it a radius, r, from one of the corners of the plate 
fig. 2) Is 

pP 
° (24) 
mrh 2 

lhe extreme fiber shearing stress in the elastic range 
s from (7), (6), and (9 
r= Gy 2Ge,=3P/h’. (25) 
The ratio of the average transverse shearing stress 
to the extreme fiber shearing stress in the elastic 

range becomes 


T 2h h 
: 0.2 (26) 
T) 77 7 
This requires h/r to be small, that is, the effect 


becomes negligible as the plate is made very thin. 

Equations (20), (23), and (26) were used to esti- 
mate the deviations from ideal behavior caused by 
edge effect, membrane strain, and transverse shear 
for square plates with width/thickness ratios b/h=5, 
10, 20. The results, rounded off to the nearest per- 
cent, are given in table 1] 


TABLE | 
hh 
1 ” 
Fidge effeet (r,,/1 1 018 0. 07 0.08 
Membrane stiffening effect: Se « 00 a2 3s 
Transverse shear effect r,/t, 2 10 Os 


In table 1 (7,./7;)—1 is taken as the correction 
according to eq (20) for the elastic case. The ex- 
treme fiber shearing strain y,=2e, is assumed as 0.016 
to place it well bevond the elastic range. The mem- 
brane stiffening effect was computed from eq (23) 
for the elastic range, since the stiffening will be rela- 
tively insensitive to the yielding that takes place 
near the extreme fibers of the twisted plate. The 


iverage transverse shearing stress was computed for 
1 radius r=6/5 from the corner leaving the inner 
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portion of the plate 0.2<r/6<0.7 available for 
strain and displacement measurements. 

Table 1 indicates that a plate with b/h=10 
provides a fair compromise between the three differ- 
ent disturbing effects that were considered. 


5. Tests 


5.1. Test procedure 


Specimens were made from aluminum alloy 75S—T6 
sheet nominally 0.125 in. thick for which the tensile 
and compressive stress-strain curves had been ob- 
tained both parallel and perpendicular to the direction 
of rolling [14]. This material was selected because of 
its approximately isotropic properties. The speci- 
mens were taken from near the center of the sheet, 
where there was little variation in thickness. The 
specimen selected for trying out the procedure was 
0.1324 in. thick and 1.25 in. square. This gave a 
width/thickness ratio b/h=9.44, which approximates 
the value b/h=10 suggested in the previous section, 
The specimen was large enough to provide space for 
mounting a deflectometer and for attaching strain 
gages at several locations. It had additional mate- 
rial at each corner, as shown in figure 7, to reduce the 
concentration of transverse shearing stress at the 
corners. A conical seat was machined at each corner 
of the square, a pair on one side at the ends of a 
diagonal and a pair on the other side at the ends of 
the other diagonal. 

The specimen 3S, figure 8, was loaded through 
1 /8-in. steel balls in the conical seats at the corners 
with a beam-and-poise screw-tvpe testing machine 
of 50,000-lb capacity used in the 5,000-Ib range. In 
this manner of loading, the load acted along vertical 
lines through the centers of the balls, which were 
almost a thickness of the sheet from the neutral 
surface. This resulted in the moment becoming 
larger instead of less as the specimen twisted. This 
effect is opposite in sign to the stiffening effect of 
membrane stresses previously mentioned. Two cor- 
ners of the specimen rested upon balls. One of these 
was seated in a block of steel resting on the upper 
head of the testing machine, and the other was seated 
in the block A, figure 8, which was free to move on 
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Figure 7. Location of wire strain gages on plate. 














Fiavure & Plate under load with deflectometer in position for 


measuring deflections along principal diagonals 


rollers in the direction of the diagonal The other 
two corners were loaded by means of a hanger con- 
nected to the loading head of the machine by a rod 
having a ball-and-socket joint at each end. The 
upper cross pieces, B, in which the balls were seated 
are shown in figure 8. They were connected to a 
common piece below the head by steel tapes,C. The 
ball joints at each end of the rod and the pin joints 
at each end of the tapes had the effect of equalizing 
the load on the corners of the plate 

SR-4 type A-7 strain gages were mounted on the 
specimen in pairs on opposite sides of the specimen, 
The gages were connected with an SR-4 model K 
portable strain indicator through a multiple-pole 
mereury-contact’ switch. A gage cemented to a 
piece of the sheet material was used as a temperature 
compensating gage. Since each portion of the spec- 
imen was subject to both tensile and compressive 
strains at right angles to each other, the gage-factor 
setting was taken as the axial sensitivity of the gage 
computed from the gage factor provided by the 
manufacturer and the data on axial and transverse 
sensitivities for gages of this type [15]. Strain 
readings were estimated to 0.1 division (0.000001). 
A correction was applied to compensate for the strain 
sensitive wire being farther from the neutral surface 
of the specimen than the surface of the sheet; each 
strain was multiplied by the ratio of the thickness of 
the specimen to the sum of the thicknesses of the 
specimen and the gage. 

A deflectometer using two Tuckerman optical 
strain gages with 0.4-in. lozenges, figure 8, was used 
to measure deflections relative to the center of the 
plate. The Tuckerman strain gages were mounted 
so that their knife edges were connected to a pair of 


legs, D, placed along one diagonal of the plate, while 
their lozenges were connected to a pair of shorter 
legs, E, placed along the other diagonal. One of the 
gages was cemented to the cross piece holding the 
two longer legs. The knife edge of the other gage 
was mounted on a ‘-in. tube which was inserted in 
the cross piece midway between the two legs. The 
lozenges were actuated by a ‘j-in. rod, which was 
screwed into the cross piece holding the two shorte: 
legs and which was reduced at the lower end to pass 
through the tube with ample clearance. The loz- 
enges were set near their lower limits of travel so that 
the vertical component of the clamping force would 
assist in maintaining contact of the two shorter legs 
with the specimen. The deflectometer was mounted 
on the specimen symmetrically with respect to the 
center with the point of one of each pair of legs in a 
small punch mark on each diagonal and the other in 
a scratch along the diagonal. The span was about 
0.7 in. The deflectometer measured the sum of the 
deflections at the center with respect to each pair of 
points or twice such a deflection, since the loading 
was equal and opposite at each pair of corners. With 
this deflectometer, deflections could be read to the 
nearest 0.000002 in 


5.2. Test Results 
a. Elastic Range 


Wire strain gages were cemented on the specimen 
at locations | and 3, and 2 and 4, respectively oppo- 
site them on the other side; figure 7. The speci- 
men and deflectometer were set up as shown in 
figure 8, and the specimen was loaded in the elastic 
range. The value of the shear modulus is from (6 
(12) 

(97 


yi 2e, Zw h 


where 7), is the shearing stress after correction for 


edge effect. Within the elastic range, from (20) and 


(25), 
3P . 
T\ , ~¢ (<8 
hl 0.63h 6b) 
or with b/A=9.44, 
T 3.21P/h (29 


Substitution of r,. in (27) gives the shear modulus 
in terms of either measured extreme fiber strain «, 
or measured center deflection w,. Table 2 shows the 
values thus obtained from the initial linear portion 
of plots of deflectometer readings and wire strain-gage 

readings as a function of load P. 
The last value of @ in table 2 was computed from 
G a (30 


2(1+-v) 


| by substituting for the Young's modulus an average 
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Shear-modulus values measured at central portion 
of square plate 


Obtained fron a 
thin 
Deflectometer 3.82 10 
Strain gage | t.S7 
Strain gage 2 +82 
Strain gage 4 tM 
Strain gage 4 LAA 
Equation (30 ww 


value = 10.37 10° lb/in? obtained from tensile and 
compressive tests of the aluminum alloy 755-T6 
sheet. and for Poisson's ratio »=0.33, the nominal 
value for this alloy 

Table 2 shows that the experimental values of 
shear modulus agreed within 1.5 percent. The 
average for the wire gages was about 0.7 percent 
ereater than the value obtained from the deflection 
data and about 1.2 percent less than the value com- 
puted from eq (30 A lower value from the deflee- 
tion data would be expected because of the effect of 
transverse shear. 

The deflectometer was removed, and gages 5 and 7, 
figure 7, and gages 6 and 8, respectively opposite 
them, were cemented to the specimen. The loading 
was repeated in the elastic range. The results are 


given in table 3. 


larie 3 Nhea modulus values measured in outer 


portion 


oj square plate 


Obtained 


from gages 


The results show greater strains at locations 5 and 
(i than at any other location. This indicates a devi- 
ation from the uniform state of stress, which becomes 
serious in the corner portions of the plate 


b. Plastic range 


The specimen was loaded until a shear strain of 
about 2 percent was reached at the extreme fiber. 
Deflections were measured with the deflectometer 
and strain with the four central wire gages. The 
shear stress was partially corrected for edge effect 
by substitution of r,, eq (16), in eq (20). Equation (20) 
was used bevond the elastic range to avoid discon- 
tinuity at the end of the elastic range. It should be 
an adequate approximation in the important region 
near the knee of the stress-strain curve because only 
a small proportion of the material would be in the 
plastic range. 

The stress-strain data are plotted in figure 9. 
The values of G determined both from the strain 
data and deflection data were a little over 2 percent 
higher than in the first test, the average value for 
the four wire gages being alntost 1 percent above the 
value computed from (30). No explanation could be 
found for this small deviation. 
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c. Comparison with stress-strain curves under axial! load 


It is interesting to compare the stress-strain curve 
in shear with that under axial load. Such a compar- 
ison can be made by deriving the octahedral stress- 
strain curve in each case [16]. The octahedral stress- 
strain curves for both states of stress should coincide 
if the maximum shear strain y, on the octahedral 
plane (the plane defined as having a normal which 
makes equal angles with the three axes of principal 
stress) is a single-valued function of the maximum 
on the octahedral plane, regardless 
and proportion of the principal 


shearing stress r 
of the direction 
stresses. 

The octahedral stress-strain curve 7,(y,) 
puted from an average axial stress-strain curve o(e) 
by substitution in the following formulas D6): 


Was coln- 


“oot 2(« >) (31) 
G _ ik 

The axial stress o was obtained from tensile and 
compressive stress-strain curves of the material 
parallel to and perpendicular to the direction of 
rolling [14] as the average of the tensile and com- 
pressive stresses corresponding to a given axial 
strain e. The value of G used, 3.93 10°lb/in, was 
the average determined from the strain data in this 
test. The value of / used, 10.37% 10° lb/in’, was 
the average obtained from tensile and compressive 
tests on the sheet. The stress-strain curve in shear 
determined from the strain data, figure 9, was con- 
verted similarly into an octahedral stress-strain 
curve by substitution in 


2 ; Toy 2 - 7) 29 
ye "ery 3 | ait ) (a2) 
The resulting curves of octahedral stress versus 
total strain y, and versus permanent strain y,—7r,/@ 
are plotted in figure 10. It will be noted that the 
two stress-strain curves are in good agreement ex- 
cept for deviations up to 5 percent near the knee and 
at shearing strains exceeding 0.009. 

In this connection, mention should be made of 
the much larger discrepancies found in exploratory 
twisting tests of a 0.125 0.625 0.625-in. plate of 
SAE 1025 steel. These tests were made to extend 


| the above comparison to a material with a pro- 
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Fiat re 10. Octahedral stress-strain curve, aluminum alloy 
758 T6 sheet 


nounced vield point. The shear modulus was found 
to be in good agreement with the value computed 
from Young’s modulus and Poisson’s ratio. How- 
ever, the octahedral stress corresponding to the y ield 
point obtained from the twisting test was 38 percent 
greater than that obtained from the average of the | 
vield points of tensile and compressive specimens of 
the material. Twisting tests of long strips of 
width/thickness ratios equal to 5 and 10 gave 
octahedral stresses corresponding to the vield point 
respectively 17 and 26 percent greater than for the 
average of the tensile and compressive tests. A 
twisting test of a long round specimen 0.59 in. in 
diameter, on the other hand, showed a corresponding 
octahedral stress 15 percent lower. 


6. Conclusion 


The stress-strain curve in shear up to strains of 
the order of 0.01 can be obtained for isotropic flat 
sheet from measurements of deflection or of extreme 
fiber bending strain in the center portion of a twisted 
square plate, provided that a correction is made 
for the edge effect. The shear moduli obtained from 
the twisting test of a 0.13% 1.25 1.25-in. square 
plate of 755-T aluminum alloy agreed within 2 
percent with the value calculated from the average 
of the tensile and compressive moduli and an accepted 
value of Poisson's ratio. The octahedral stress de- 
rived from the stress-strain curve in shear for the 
twisted plate agreed within about 3 percent with 
the octahedral stress derived from tensile and com- 
pressive tests of the same material except near the 
knee, where the systematic deviations were as high 
as 5 percent. Much greater deviations were found 
in exploratory tests of mild steel. 

Unfortunately, the method is difficult to perform 
on sheet material because the specimen must be 
small in size. Alternative methods should be inves- 
tigated. An examination of other methods indicates 
that the annulus of constant thickness method can 
be used for determining the stress-strain curve in 
shear, even in the presence of a large radial variation | 
in shear stress. Preparations are under way for | 
carrying out such tests at the National Bureau of | 
Standards. 





7. Appendix 


7.1. Membrane stresses in elastic range 


The membrane stresses in a twisted square plate 
can be estimated from the equations for the state 
of strain in a flat plate with large deflections. These 
are given by Love [17] as 
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It is assumed that the plate of figure 2 is bent into 
an anticlastic surface so that the transverse deflection 
is given by 


r?-—y* ; 
—y (34) 
2R 


w 


| just as in the case of small deflections, eq (3). 


The curvature is from (11) 


1 2e, ez 
R h- (35) 
The magnitude of the membrane strain depends on 
the displacements u, ¢ parallel to the z, y axes as 
well as on w. 

The displacement u in the s-direction can be 
approximated by a polynomial of the form 


‘ on 
u= > e,,2'y’, (36) 
‘J 


where 7, j are positive integers. The coefficients 
c,, of many of the terms in the series must be equal 
to zero in order to satisfy the following conditions 
of symmetry and antisymmetry in the plate, figure 2: 


ur, yjy=—usr,—y) 
(37) 


~~ 


ur, y) -u(—r,Yy) 

These conditions are satisfied if 7, j are confined to 
=135... 

(38) 
j=0,2,4... 


The displacement r at the point (y, z) must be equal 
to the displacement u at the point (z, y), so that 


e= Deyy'r’. (39) 
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Up to and including terms of the fifth power, the 
series (36) and (39) become 

U= Col + Cry? + Cyt? + Cygry' + Cyr? + Cyt” ) 
>. (40) 
C=C + Crt Y + Ceo’ + Crt ty + Cot? + Coy? 
To satisfy the conditions that the membrane stresses 


Try be zero along the edge 
(41) 


of the plate, and that the strain energy stored in 
the form of membrane stress must be a minimum, 
the six unknowns ¢j, ¢), and Cs) were 
determined as follows. 

Equilibrium of forces on an element of the plate 
bounded by the edge and by planes y=const., 
r=const. shows that the membrane stresses along 
the edge are related to the membrane 
by 


c Mie 


Ci4, C32, 


stresses 


Ce, Gn, Fu 


(42) 


The stresses ¢,, o), Tt, must be adjusted to make 
‘ | — 9 . | 7er The »mbr. . 
tT, and tyy mm (42) equal to zero. 1 membrane 
stresses are related to the membrane strains ¢, «&,. 
yx by Hooke’s Law 


BE \ 
Ts: " 
~ y) 


where /, vy are Young’s modulus and Poisson’s ratio, 
respectively. 

The membrane strains may be computed from 
the displacements, using eq (33). Carrying out the 
computations and using the following notation for 
convenience 


e__.* Y. b 
ie ee a 
b (44) 
Co Cho; t Cc; Lb ( 3eab?+ 9 Re 
¢,=¢,,b* ¢ Cy 6': «¢ Ca,b' 
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in (43) and (43), in (42) and using 


Substituting (45 
the notation 


(46) 


gives, for the normal and shearing stresses acting 
on the edges, 


0, oly (Cy 


+. (47) 


(4¢,+-2e, 


E i 
r; (¢ 


+ *) 
2(1-+-v) m 


4 


+ (5e,—¢.(& 


n*)|(¢3 Cs 





The shearing stresses r,,,, will be zero for all points 
£, n along the edge of the plate only if the expression 
in the square brackets is equal to zero. This re- 
quires that 

Cy (48) 
The normal stress ¢,- along the edge can be expressed 
as a polynomial of the fourth order in & by noting 
from (41) that, for points along the edge of the plate, 


| 
” 3 (49 
\2 
Collecting terms by powers of & and substituting 
(48) leads to 
2(1—v) 5 N 
K C; 2¢,-+ Cst5 C.+ = | 4¢,— 20¢,-4 
) 2 ¥2 
a(4e,4+- 10¢,+-¢,— p*)| + Fl4e3 + 30¢;-4 
3¢4+ a( —4c,—300¢,—3e,4+4- p*)|4 
\ 2&°| — 20¢,-- 60,4 al(4e,+ 40¢€;)| 4 
£*(20c,+ 6c¢,+ a(—4ce,— 40c;)|=0 (50) 


The stress ¢, will be zero for all points & along the 
edge only if each one of the five brackets on the 
right is equal to zero. Examination of these brackets 
shows that only three of them are independent of 
each other. One may solve the resulting three 
equations for any three of the coefficients ¢o, ¢s, 
Cs, Cg in terms of the remaining coefficient. 

Solving in terms of ¢> and assuming for Poisson's 
ratio the typical value »=1/3 and hence from (46) 


a= 1/2, one obtains from (50) 


2 a 2 
3¢o + p/8; Cs 5 Co- a0 °? (=O. (51) 


Substitution of (51) and (48) in (45) shows that the 








Strains reduce to 


Yay = [4es— p? + 200,(8 4 
It remains to compute the constant ¢) by minimizing 
the membrane strain energy stored in the plate. 
The membrane strain energy per unit volume, 
Vom More, + oye, + Tyr) 18 With (43) and (52), assum- 
ing vy=13 

, S yp : : Rs 

Ve . (Se ¥ (53) 

16 g 


The elastic membrane strain energy \’,, stored in the 
plate is obtained by integrating (53) over the volume 
of the plate 


Vohd yds 


noe fe fe * Vidnds. (4) 


Substitution of (52) in (53), of the resulting expres- 
sion in (54), and integration leads to 


\". = bh? [ 268 n : Col s+ I (c2 + 10¢0es)4 


The membrane energy V,,, is minimized by computing 


OV’, 
OC» 


0 (56) 


and substituting (51) to obtain an equation having 
¢o as the only unknown. Solving for ey gives 


3 


146 ps 0.0205 p*. (57) 
a0 


€o 


The membrane strains and stresses throughout the 
plate are then obtained by substituting (57) in (51) 
and the resulting constants in eq (47) and (52) and 
the latter then in (43). 

The stresses can be plotted in dimensionless form 
by multiplying them by A’/b*y,r;, where b/A is the 
width/thickness ratio for the plate and y,,r; are the 
extreme fiber shearing strain and stress due to 
bending. These are from (6), (7), and (11) 


h. al Se 
att! el | (98) 


The results of such computations are plotted in fig- 
ure 11 as solid lines. The plots show that the mem- 


brane stress at the center of the plate consists of a 
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uniform compressive stress 


b* se 
Oo, C, 0.082 he Vi Ti- (O90) 


The membrane stresses are largest at the center of 
the edge of the plate. The membrane normal 
stresses at this point are 


b? : 
fa Vets CO, 0. (60) 
52 


Cy 0.136 
A check on this solution seemed advisable since it 
did not satisfy the conditions of equilibrium at all 
points of the plate; it would be necessary to apply a 
system of body forces per unit volume with compo- 
nents .Y, )’ in the vs, y directions in order to maintain 
equilibrium in all plate elements. The equilibrium 
conditions would be satisfied more nearly if ¢ is 
chosen to minimize the average value of \?+ ¥? for 
the entire plate instead of minimizing the strain 
energy. This leads to 


C O.OLSTp*, (61) 


which is about 10 percent below the previous value. 
The membrane stresses corresponding to this value 
of ¢9 are shown as dotted lines in figure 11. They 
differ less than 10 percent from the stresses previously 
obtained. 

A second check on the adequacy of the theory was 
obtained by computing from it the nonlinear relation 
between load 7? and extreme fiber bending strain ¢, 
and comparing it with the experimental values given 
in figure 6. The relation between 7’ and e, for the 
elastic range was obtained by equating the incre- 
ments of work done by the corner loads dW) and the 
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Fieure 11, Elastic membrane stresses in twisted square plate. 























he increment dV), stored in the form of bending 


rains and the increment dV), stored in the form of 


embrane strains 
J d¥ 1-8... (62) 


‘ach one of the three terms in this equation can be 
egarded as a function of the curvature ratio p=6 BR. 

can be shown easily that the work done by the 
orner forces ? corresponding to an increment of 
irvature ratio dp is 


dW, Phdp (633) 


lhe increase in bending strain energy can be derived 
om the bending strain, eq (4) as 


dV, a od p (4 


lhe increment in membrane strain energy is obtained 

-~- » | -—-_ - ry’ ° 
from (55) after substituting (57) and (51). This 
leads to 


d\ we 0.001146 Ehb*p dp (65 


Substituting (63) to (65) in (62), solving for 7’, and 
replacing p by 
h b 


» 2 
p y €\: (Ob) 
R h 


vives the desired relation for ?(e,), 


pe } i 
P = afi - 0.03666 (7) } (67) 


This reduces to eq (9) for small values of €; with 
13. Substitution of A=10.4%10° Ibinw, A 
0.126 in., b=7.2 in. for the plate of figure 6 leads to 
the solid line shown in that figure. The observed 
points come close to this curve, thereby providing a 
second check on the adequacy of this approximate 
theory for computing membrane stresses in the 

elastic range. 





esulting Increase in strain energy dV" consisting of 
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It is shown that the stress-s train relation in shear of isotropic thin sheet can be deter- 


mined from a test of an annular specimen 


The annulus must be clamped uniformly along 


the inner and the outer edge, and the relative twist of two circles concentric with the edges 


must be measured 


The shear strain at the inner circle can be computed from the shear 
strain at the outer circle and the slope of the torque-twist curve 


Tests were made on 0.032- 


inch aluminum alloy 75S—T6 sheet well bevond the elastic range to a shearing stress of 48,000 


Ib/in2 


EPON adhesive VI was used to bond the sheet to a clamping ring at the outer edge of 


the annulus and between the sheet and the end of a cylinder through which the torque was 


applied. 
1. Introduction 


relation in shear of thin sheet 
tal is of practical importance for estimating the 
strength of shear webs. It is of theoretical impor- 
anee in studying the plastic behavior of metals 
inder combined stress 

Many methods have been suggested for determin- 
ng the stress-strain curve in shear of thin sheet [1], 
and at least two of them have been tried. Good 
esults have been obtained so far only from twisting 
tests of square plates of aluminum alloy [1]. Unfor- 
tunately the twisting test of a square plate makes 
severe demands on the skill and patience of the test 
ngineer. Furthermore, the specimen becomes too 
small for precision in the measurement of strain and 
n the application of the twisting loads for sheet less 
than O.1 ia. thick. Hence it seemed desirable to try 
out some of the other methods that have been pro- 
One of the most promising among these is 
the twisting test of an annulus of constant thickness. 
This method was accordingly investigated at the 
National Bureau of Standards as part of a program 
supported by the Office of Naval Research, United 
States Department of the Navy. 


| he stress-strain 


posed. 


2. Theory of Twisting Test of Annulus of 
Constant Thickness 


Consider an annulus of constant thickness A, 
twisted by a torque .V/,, which applies a constant 
shearing foree per unit length along the inner edge 

a, figure 1, and along the outer edge r=6. If the 
material is isotropic, the shearing stress at any radius 

is given by 


me - 


2arh 


The shearing strain at any point can be derived 
from the radial displacement u and the cireumferen- 
tial displacement ¢ by substitution in the equation 
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for plane strain in polar coordinates r, @ |2, p. 66] 


or. Ou r 
Or rod 13 


In an isotropic annulus the displacements u, 7 caused 
by twisting are independent of @, so that (2) reduces 
to 
dr a d r 
y=tt_tard(*) 7 
"dr or drX\y 


Integration of this expression leads to a formula for 
the circumferential displacement 6 of a point on an 
outer cirele with radius r, relative to a platform 
resting on an inner cirele with radius 7), figure 1. We 
note from figure | that this displacement is 


Dividing y by r in (3) and integrating gives 


7” i ] * 
| di (.) 
Jr? ; / 
so that 
5 "2 4 
c y 
dr (b 
/ Im 7 


From this an equation for the shearing strain y; at 
the circle 7 r, can be obtained as follows. We can 
replace ras the independent variable in the integral 


by r, by taking the logarithms of the two sides in 
eq (1 
log r=log — —? log r; (7) 
Jah 


and then differentiating for a given plate under given 
twisting moment (\/,, A constant): 


dr dr 
2 . (S) 
T r 




















Fieurr | Diagrammatic sketch of annulus 


radius of inner gage circle; ro, radins of outer 


gage circle 


a, inner radius: 6, outer radius; r 
° r ) rj . 
Inserting (8) in (6) gives 


(Q) 


where r,, tz are the shearing stresses at radii 7, 72, 


respectively. Differentiating with respect to 17, 
gives [3, p. 29] 
dé ry *r ad ¥ 
dM, 2 [ I. au A r )dr+ 
dr, V2 dr, ¥1 
dM,\ 1, ~aM| r )} (10) 


Since the ratio of y to r does not depend on the para- 
meter AM,, d(y rT) dM, 0, and consequently the 
integral on the right-hand side is equal to zero. The 
derivatives dro/dM, and dr,/dM, are obtained di- 
rectly from (1), noting that r is equal to r;, r, respec- 





tively. This leads to 
db ry YY ) 1) 
AM, Amrh\ nr ner 


This equation can be solved for the shearing strain 
y; at the inner cirele 


dé 
‘dM, 


(12) 


” 
n=nt M 


Equation (12) makes possible a step-by-step con- 
struction of the stress-strain curve in shear, provided 
that the measurements are started within the elastic 
range. The shear modulus G may be determined 
from a single measurement of 6 in that range by | 
noting that inside the elastic range 


y vy '/ : 
fo % G (13) 
Inserting in (11) and solving for G: 
rs l l dé 
G on 
, jek (A 2) dM, (14) 
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3. Procedure for Twisting Test . 
3.1. Dimensions of Specimen 


The dimensions of the specimen must be chosen to 
satisfy several opposing requirements. 

The required external torque must be readily at- 
tainable. If M, is the maximum torque that can 
be applied and +, the shearing stress at the inner 
edge r=a of the annulus, we have from (1) 


M, 
\ Qar,h 


a (15) 


The torque must be applied uniformly to the 
edges of the specimen. This makes it inadvisable 
to use clamps consisting of heavy clamping rings 
held together by individual bolts. In such a clamp- 
ing arrangement one would expect higher shearing 
stresses to be transmitted near the bolts, where the 
clamping pressures are high, than at locations rela- 


tively far from the bolts. It would be better to cut f 
a circular disk with a radius well in excess of 6 from 
the material and then apply the torque to one or { 


both faces of this disk through adhesive lavers, as 
indicated in figure 2. The torque \/, is then limited 
by the shear stress 7, transmitted through the adhe- 
sive at the inner edge r=a of the annulus. 

The magnitude of +, can be computed from the 
torque M, on the assumption that each adhesive 
laver behaves like an elastic lamination of constant 
thickness connecting two rigid bodies that are ro- 
tated relative to each other. The shearing stress is 
then proportional to the distance r from the axis of 
twist. If there is only one layer, the torque is 
equal to 


*a a r - 
M, r2arrdr—2r rr, dr = a®*r,. (16) 
Jo J a - 


The torque \/, produces a shearing stress r, at the 
inner radius of the annulus given by (1) 


M, ay 
Tv; »2° ‘ 
‘ 2Qrarh 
With (16) this leads to the condition 
' 7 
2Qra*hr, < 5 a°r., (18 
or 
ash >4r,/r,. (19 
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__f—!—4 specimen, |) , 
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\.¢dhesive~ | } j A 
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Fieure 2. Cross-sectional view of annulus and grips. 


The sheet in the annulus must be stable against 
buckling into a pattern of shear wrinkles up to the 
highest torque applied. The authors have not found 
a solution for the buckling stress of an annulus. 
However, an approximation to this buckling stress 
is obtained from Southwell and Skan’s solution for 
the infinitely long elastic strip with clamped edges 
4]; with »=0.3 this reduces to 


o1 | Ghe 
oe ae 


The modulus of elasticity in shear, G, must be re- 
placed by a reduced modulus G, where G@ will be 
somewhere between G and the tangent modulus at 

since buckling must not take place until the 
material at r=r,; has been well into the 
plastic range. This leads to the condition 


stressed 


Gh? 
21.1 ’ 
(b—ay 


(21) 
where 7 denotes the shearing stress at the mean 
radius ¢=(a+-6)/2, of the annulus. Condition (21 
sets an upper limit to the width (b—a). 

A lower limit is set by the condition that (6—a) 
must be sufficient to accomodate the device for 
measuring relative displacement. 

Another condition that must be satisfied is that 
the tensile strain in a radial direction is negligible. 
From figure 3 we see that this radial strain is in 
first approximation 


y(b—a)y?+(r 
(bh—a) 


or with 


(24) 


This must be small compared to the normal strain 
of the order YA along lines at 45° to the radial 








' 
TTTT/11117 


Ficure 3. 
to relative circumferential displacement of inner and outer 
edges of the annulus 


Diagramatic sketch for estimating radial strain due 
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direction, that is, 


From (24) 


y » (26) 
Which is small compared to one in the range in 
which we are interested. 

Let us apply these equations to the specific case 
of high-strength aluminum-alloy sheet with a vield 
strength in shear, 7,, of the order of 45,000 Ib in 
(tensile yield strength about 78,000 Ib/in) to be 
tested in a torsion machine with a capacity of 27,000 
lb-in. Let us assume that a bonding material with 
a shear strength 7,=4,500 Ib/in2 can be used for 
applving the torque to the specimen, and that a 
displacement measuring device with a width of 1 in 
can be constructed to measure 6, so that b—a>1 in 
Let us assume also that the stress at the inner edge 
exceeds the vield strength in shear by 20 percent 
54,000 Ib/in2 Substituting these quan- 
tities in (17) (19) gives a<0.282/yh and 
a>48h, respectively, or, 484 <a<0.282/yh. We 
obtain an upper limit for the thickness h, 4 <0.0326 
in., by going to the case in which @ is equal to its 
lower and its upper limit. The corresponding radius 
of the inner clamp is a=0.282/,0.0326=1.56 in 
The outer and the mean radius become b=a+1 
2.56 in; c=2.06 in. The shearing stress at the 
mean radius is from (1) 7 =—7,(a?/e?) =31,000 Ib/in- 

This calculation indicates that the annulus must 
have a thickness of less than 0.0326 in. for high- 
strength aluminum-alloy sheet with a vield strength 
in shear of the order of 45,000 Ib/in 

A lower limit for the thickness may be obtained 
by solving (21) for hk. Since much of the material 
in the annulus will be in the elastic range, it seems 
reasonable to assume 3,000,000 Ib/in? (about three- 
fourths the modulus of elasticity in shear) as an 
value of G. Then (21) 


Te= 1.27; 


and 


average from 


h>(b-a)y 7/21.1G=0.0221 in. 


We conclude that the thickness should be between 
0.022 and 0.033 in. for this alloy tested under the 
assumed conditions. 


3.2. Method of Applying Torque 


Advice regarding the bonding medium for trans- 
mitting torque to the specimen was obtained from 
experts in the field, and it was decided to try cements 
of the Epoxy type upon recommendation of N. J 
DeLollis of the National Bureau of Standards and 
D. M. Alstadt of the Lord Manufacturing Co., 
Erie, Pa. 

The shear strength of the cement was determined 
from torsion tests of specimens made by joining in a 
butt joint two pieces of aluminum-alloy rod = 1.5 
in. in diameter. Similar specimens were made from 
aluminum-alloy tubing 1.5 in. in diameter and 0.094 
in. thick. The surfaces were cleaned with a 
greasing solvent and Deoxidine. The 


de- 
ling 
bonding 








medium was ePON adhesive VI with 6 parts by 
weight of accelerator A to 100 parts of the adhesive 
As soon as the adhesive was applied to one piece, 
the two pieces were pressed together until the thick- 
ness of adhesive was from 0.005 to 0.01 in. The 
pieces were held in position by clamps while the 
adhesive was cured for 45 min. at from 200° to 220 
fF. The excess adhesive was removed from the 
outer surface of each specimen. 

The maximum shear strength +r,, computed from 
these tests by substituting the torque at failure in 
16) for the rod and in an analogous equation for the 
tube, was between 6,200 and 6,800 lb/in’ for the 
tubular specimens and between 6,000 and 6,700 
lb/in.? for the evlindrical specimens 

Attempts to develop shear strength of this order 
in transmitting torque to a specimen of 0.032-in 
aluminum-alloy sheet bonded to a 3-in. aluminum- 
allov evlinder were unsuccessful. This may have 
been due in part to stress concentration in the ad- 
hesive at the inner edge of the annulus. This 
stress concentration was relieved somewhat by 
machining a circumferential notch in the evlinder 
0.01 in. from the bond and by bonding disks (A and 
B, fig. 4) of 0.020-in. aluminum-alloy sheet to the 
specimen 
The radii of A and B were 0.06 and 0.03 in., re- 
spectively, greater than the radius of evlinder C 
The evlinder, the disks, the specimen 5, and the 
ring R were held coaxially by the pin P and the jig 
J while the adhesive was being pressed out to a 
thickness of about 0.01 in. and during the curing 
of the adhesive The jig was also used for locating 
the bolt holes in the ring and for holding the ring 
while it was being machined. The bond between 
the evlinder and disk B did not fail until the value of 
r, reached 4,670 lb/in2? No difficulty was experi- 
enced with the bond between the specimen, which 
had a diameter of 6.5 in., and ring R, the inside 
diameter of which was 5 in 

Two torsion-testing machines were available for 
the tests, a Riehle testing machine of 60,000-Ib-in 
capacity and an Amsler testing machine of 150-kg-m 
13,000-lb-in.) capacity. In the Riehle machine the 
torque was transmitted to the weighing beam through 
a series of knife edges and linkages, which provided a 
very flexible mounting for one head of the machine 
The specimen was first set up in this machine 
because of its greater capacity but was moved to the 
Amsler machine when it was found that good aline- 





Fiaure 4. Expanded cross-sectional view of specimen S, disks 
A and B, cylinder C and ring R, positioned by jig J, and 
centered with pin P for bonding 





Fieure 5. Vodi fied Amsler torsion machine. 

ment could not be maintained. The alinement of 
the weighing head of the Amsler machine is main- 
tained by ball bearings. The 13,000-lb-in. capacity, 
measured by rotation of the pendulum, was insuffi- 
cient for the specimen. To increase the capacity, 
the weighing system was changed to the simple 
lever and seale type, utilizing the pendulum only 
for attaching the arm A, figure 5. The effective 
length of the lever was 48 in., and the seale was the 
spring balance B, having a capacity of 600 lb. Dur- 
ing the test the pendulum was maintained in the 
zero position so that its moment would not affect 
the torque measurement. The specimen assembly 
was mounted in the testing machine by bolting 
the ring to the plate H, figure 6, attached to the 
dynamometer head of the machine. The holes in the 
plate had been located by the jig J, shown in figure 
4 The evlinder was centered in the loading head 
of the machine 


3.3. Method of Measuring Deformation 


The device for measuring deformation of the 
specimen is shown in position in figure 7. The 





| 
| 


Ficure 6. Specimen ready for test in modified Amsler torsion 
machine. 
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FIGURE 7 Gage for measuring relative circumferential 


di s plac ement, 


surpose of the device is to measure the relative 
otation between two circles concentric with the 
ixis of the cylinder C having radii r; and 7, of 2.375 
and 1.625 in., respectively. These circles are defined 
by conical points, such as P, 120° apart, mounted 
m the rings O and I, respectively. The relative 
rotation is measured by measuring the relative 
ingential motion between the two rings. To 
nerease the sensitivity and reduce any errors due 
o bending of the specimen, the motion between the 
ends of the arms B and F is measured. Due account 
s taken of the radius at which the motion is meas- 
ired. Each arm B is attached rigidly to the outer 
ng O and provides physical support to the body 
fa Tuckerman optical strain gage T. Each arm F, 
attached to the inner ring I, is flexible in the direc- 
tion of the axis of the specimen and rigid in the 
direction of rotation. The rotating lozenge L of 
each gage is held against the outer end of an arm F 
by a pair of hairpin springs, one attached to each 
end of a rod through a hole in F and to a knob on 
each side of the gage. Guides G were rigidly at- 
tached to each arm B to limit the motion of the end 
of arm F to a plane perpendicular to the axis of the 
gage lozenge. The circumferential displacement 6 of 
the cirele with radius r; with respect to the circle 
with radius 72 is given by 


where e is the average displacement measured by 
ihe gages, and r, is the average of the radii to the 
points of contact of gage lozenges with the arms F. 


4. Results 


A disk of 0.032-in. 75S-T6 aluminum-alloy sheet 
as prepared and tested in shear as described above. 
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Stress-strain curve in shear for a sheet of }45S—T6 
aluminum alloy 0.0382 in. thick 


Ficgure 8 


The test portion of the disk consisted of an annulus 
with inner and outer radii of 1.56 and 2.5 in., respee- 
tively. The test was carried to a shearing stress of 
48,000 |b/in and a shearing strain of more than 
0.05 at the inner gage ring before the adhesive at 
the inner edge of the annulus broke at a shear stress 
in the adhesive, as computed from (16), of 4,670 
lb/in. A second annulus reached a shearing stress 
of over 45,000 Ib/in.®. The data for both annuli up 
to a strain of about 0.019 are plotted in figure 8, 
together with a stress-strain curve in shear computed 
from an average stress-strain curve obtained from 
tests of tensile and compressive specimens from the 
same sheet. This curve was computed, in accord- 
ance with octahedral shear theory [5] from 


wd 
Fiyo 


and y= 3 ( € ep ie ): 


where o is the average of the stresses corresponding 
to the strain ¢ obtained in tensile and compression 
tests of specimens in the direction of rolling and 
across the direction of rolling of the sheet, and £ is 
the average of the moduli of elasticity in tension and 
compression. The shear data are in good agreement 
with the derived curve. The modulus of elasticity 
in shear for each specimen differed by less than | 
percent from that computed from / and the nominal 
value of Poisson’s ratio (0.33) for this alloy. 


5. Conclusions 


The twisting test of an annulus appears feasible 
for determining the stress-strain relation in shear of 
isotropic sheet metal in thicknesses below those which 
can be tested by twisting a square plate of the ma- 
terial. Special precautions must be taken to obtain 
a good bond for transmitting the shearing stresses to 
the specimen and to have good alinement between 
the two heads of the torsion-testing machine. 

Using this method, the stress-strain curve in shear 
was determined for 0.032-in. aluminum alloy 75S-T6 
sheet up to a strain of 0.05, which is far beyond the 
elastic range. .The stress-strain curve thus obtained 
was in good agreement with one derived from the 
stress-strain curve of the same material in tension 


| and compression, on the assumption that the octa- 
| hedral shear stress in the metal was a function of the 


octahedral shear strain only. 
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